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Preface 


The  n-hit  simulation  tool  has  been  developed  with  the 
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thanks  to  my  two  advisors,  Dr  Peter  Maybeck  and  Captain  J B 
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Thank  you  for  your  patience,  understanding,  support,  and  your 
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In  the  development  of  estimation  and  control  systems, 
the  algorithms  developed  are,  in  general,  very  sensitive  to 
the  wordlength  of  the  onboard  computer.  The  representation 
of  the  desired  design  can  be  greatly  affected  by  the  quanti- 
zation resulting  from  truncation  or  round  off.  The  result 
can  be  even  more  severe  than  niimerical  imprecision:  numer- 
ical instability  can  be  generated  that  render  algorithms  to- 
tally useless. 

An  importcint  part  of  the  algorithm  development  is  the  - 
determina'tiron  of  ” the  appropriate  wordlength.  The  proper 

K 

wordlength  can  be  determined  by  running  the  algorithm  on  a 
simulated  n-bit  machine  for  several  values  of  n.  The  re- 
suits  can  then  be  evaluated  against  the  performance  speci- 
fication to  determine  the  wordlength  and  precision  require- 
ments . 

The  objective  of  this  thesis  was  to  process  algorithms 
written  in  FORTRAN  on  the ''Control  Data  Corporation  (CDCX”* 
6600/CYBER  74  computer  systems  such  that  the  results  ob- 
tained were  similar  to  those  obtainable  on  an  n-bit  machine. 
The  problem  of  n-bit  simulation  was  addressed  from  two  lev- 
els: the  assembly  language  level  and  the  FORTRAN  language 
level.  Each  level  involved  designing's  preprocessor^hich  , 
would  modify  the  algorithm's  code  so  that  the  numerical  ef- 
fects of  n-bit  wordlength  could  be  realized . ..  The  assembly 
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language  (COMPASS)  level  approach  failed,  however,  an  n-bit 
simulation  tool  was  successfully  developed  and  implemented 
at  the  FORTRAN  level,  ^Several  user  options  were  incorporated 
into  the  n-bit  simulation  tool  to  enable  the  user  to  simulate 
the  characteristics  of  various  computer  types. 
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I,  Introduction 


In  the  development  of  estimation  and  control  systems, 
the  algorithms  developed  are,  in  general,  very  sensitive  to 
the  wordlength  of  the  onboard  computer.  The  representation 
of  the  desired  design  can  be  greatly  affected  by  the  quanti- 
zation resulting  from  truncation  or  round  off.  The  result 
C6in  be  even  more  severe  than  numerical  imprecision:  numer- 
ical instabilities  can  be  generated  that  render  algorithms 
totally  uselt 

An  important  part  of  the  algorithm  development  is  the 
determination  of  the  appropriate  wordlength.  The  proper 
wordlength  can  be  determined  by  running  the  algorithm  on  a 
simulated  n-bit  machine  for  several  values  of  n.  The  re- 
sults can  then  be  evaluated  against  the  performance  speci- 
fication to  determine  the  wordlength  and  precision  require- 
ments . 

The  objective  of  this  thesis  is  to  process  algorithms 
written  in  FORTRAN  on  the  Control  Data  Corporation  (CDC) 
6600/CYBER  74  computer  systems  such  that  the  results  obtain- 
ed are  similar  to  those  obtainable  on  an  n-bit  machine. 

Then  the  results  would  reflect  the  degree  of  accuracy  that 
would  have  resulted  from  actually  processing  the  algorithm 
on  a computer  with  that  n-bit  wordlength.  For  the  purposes 
of  this  report,  n-bit  simulation  will  refer  to  the  effect 
n-bit  wordlength  would  have  upon  the  numerical  accuracy  of  an 
algorithm  being  processed  on  an  n-bit  wordlength  machine. 
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This  type  of  tool  would  be  useful  in  the  initial  stages 
of  the  selection  process  for  an  onboard  computer  that  would 
process  the  desired  algorithms.  It  is  essential  to  know  the 
full  impact  of  the  wordlength  of  a computer  upon  the  algo- 
rithms before  selecting  the  computer. 

For  a hypothetical  exaunple,  suppose  the  Air  Force  just 
purchased  several  mini-computers  having  16-bit  wordlengths. 
It  was  thereafter  discovered  that  these  computers  did  not 
furnish  the  degree  of  accuracy  required  by  the  algorithms. 
Considerable  time  may  be  expended  rewriting  the  algorithms 
to  process  properly  on  the  purchased  computers,  if  that  is 
at  all  possible.  In  fact,  it  may  be  the  case  that  the  min- 
imum wordlength  which  could  meet  the  requirements  would  be 
24-bits  or  that  a 16-bit  wordlength  computer  with  double  pre 
cision  capabilities  would  be  required.  Application  of  the 
n-bit  simulation  tool  before  purchase  could  have  saved  the 
government  significant  time  and  money. 

Similarly,  at  the  other  end  of  the  scale,  this  software 
tool  could  save  government  expenditures  for  computers  with 
much  longer  wordlengths  than  the  algorithms  could  ever  re- 
quire. Computer  costs  go  up  rapidly  with  longer  wordlengths 
Tailoring  of  the  computer  wordlength  to  the  algorithm  re- 
quirements could  therefore  eliminate  some  unnecessary  ex- 
penditures. 

Through  the  utilization  of  the  n-bit  simulation  tool, 
an  algorithm  can  be  processed  with  several  different  word- 
lengths.  The  results  can  be  compared  to  results  generated 
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by  the  full  60-bit  wordlength  (which  is  longer  them  any  for- 
seeable  onboard  computer  wordlength)  of  the  GDC  6600/CYBER  74 
computer  systems.  The  degree  of  accuracy  lost  due  to  each 
n-bit  wordlength  cam  then  be  determined. 


Example  of  Its  Application 

Figure  1 cam  be  used  to  demonstrate  a typical  applica- 
tion of  the  n-bit  simulation  tool.  The  Kalman  Filter  repre- 
sents the  algorithms  which  would  be  processed  by  the  onboard 


computer.  However,  this  procedure  did  not  correctly  reflect 
the  degree  of  accuracy  the  onboard  computer  (having  a shorter 
wordlength)  would  achieve.  The  original  algorithm  (prograimed 
in  fortran)  would  be  modified  by  the  n-bit  simulation  tool 
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Figure  1 N-bit  Simulator  Application 


3 


r 


so  that  its  results  would  show  the  exact  amount  of  signifi- 
cance that  an  n-bit  wordlength  computer  would  have  produced . 

It  may  be  executed  several  more  times  with  varying  wordlengths 
until  sufficient  information  has  been  gained. 

Requirements 

The  requirement  for  this  thesis  is  to  develop  an  n-bit 
simulation  tool  which  when  applied  to  a FORTRAN  programed 
algorithm  will  produce  the  near  exact  arithmetic  accuracy  as 
a computer  of  the  specified  n-bit  wordlength.  Five  options 
were  considered  necessary  to  provide  the  user  with  a versa- 
tile n-bit  simulator  which  could  be  used  to  simulate  the 
characteristics  of  various  computer  types: 

1)  Allow  the  user  to  specify  n,  the  number  of  bits  per 
wordlength, 

2)  Permit  specification  of  floating  point  word  char- 
acteristics: 

a)  The  number  of  bits  to  be  used  for  the  exponent 
emd  memtissa 

b)  Position  of  the  implied  decimal  point  with  re- 
spect to  the  mantissa,  either  to  its  left  (mak- 
ing the  mantissa  a fractional  representation) 
or  to  its  right  (making  the  mantissa  a fixed 
point  representation), 

3)  Give  the  user  an  option  to  have  arithmetic  opera- 
tions performed  with  either  truncation  or  rounding 
effects. 
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4)  Allow  the  user  to  specify  whether  arithmetic  calcu- 
lations are  performed  in  single,  double,  or  triple 
wordlength  precision  with  the  results  being  stored 
as  single  precision  quantities.  This  option  essen- 
tially simulates  the  effects  of  a computer  having  a 

single,  double,  or  triple  wordlength  accumulator 
but  only  single  precision  storage  variables,  and 

5)  Allow  the  user  to  specify  different  portions  of  an 
algorithm  for  different  wordlengths. 

Another  characteristic  of  the  n-bit  simulation  tool  is 
that  when  an  overflow  is  detected  for  a given  n-bit  word- 
length,  the  maximum  value  that  can  be  represented  with  the 
n-bit  word  should  be  substituted  for  the  overflowed  word  and 
an  overflow  message  printed.  Overflows  are  treated  different- 
ly depending  upon  the  computer  type,  but  it  was  felt  maximum 
number  replacement  for  overflows  would  better  demonstrate  the 
numerical  effects  of  shorter  wordlength  upon  an  algorithm 
than  simply  truncating  the  most  significant  bits  of  the  value 
that  overflowed.  When  a word's  value  exceeds  the  maximum 
positive  number  the  word  could  contain,  the  maximum  positive 
number  would  be  substituted  into  the  word.  If  the  word's 
value  exceeds  the  maximum  negative  number,  the  maximum  nega- 
tive number  would  be  substituted  into  the  word. 

Some  general  goals  of  this  thesis  were  to  provide  the 
user  with  a n-bit  simulation  tool  that  would  be  reliable,  be 
easy  to  use,  be  maintainable,  provide  an  option  to  communi- 
cate overflow  error  messages  when  they  occur,  and  to  employee 
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structured  design  techniques  in  programing  so  the  resulting 
design  can  attain  certain  desirable  characteristics. 


Assumptions 

Some  fundamental  assumptions  upon  which  this  thesis  was 
based  are  listed  below; 

1)  Applying  n-bit  simulation  effects  to  all  arithmetic 
expressions  is  sufficient  to  simulate  an  entire  pro- 
gram effectively  in  n-bit  precision.  The  major  pur- 
pose of  this  thesis  is  to  provide  a tool  which  can 
be  used  to  evaluate  the  effects  that  varying  word- 
lengths  have  upon  arithmetic  accuracy  of  an  algo- 
rithm. 

2)  The  algorithms  which  are  to  be  n-bit  simulated,  are 
programable  in  American  National  Standards  Institute 
(ANSI)  FORTRAN  £uid  are  executable  on  the  CDC  6600 
or  CDC  CYBER  74  computer  systems.  Potential  users 
currently  use  the  FORTRAN  language  and  execute  pro- 
grams on  these  computer  systems  regularly. 

3)  The  maximum  fixed  point  wordlength  desired  to  be 

simulated  is  48  bits.  This  is  the  largest  fixed 
point  representation  the  CDC  system  can  handle.  In 
addition,  the  largest  exponent  and  mantissa  that 
would  be  required  are  11  and  48  respectively.  These 
are  also  maximums  for  CDC  single  precision  floating 
point  representations.  Within  the  forseeable  fu-  < 

ture,  onboard  computers  would  not  exceed  any  of 
these  maximums. 
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4)  Additional  core  storage  and  execution  time  required 
for  n-bit  simulation  will  not  be  a limiting  factor 
in  the  design  of  the  n-bit  simulation  tool.  This  is 
partially  true  because  this  n-bit  simulation  tool  is 
expected  to  be  used  only  several  times  each  year  by 
any  projected  user. 

5)  The  quantities  to  be  n-bit  simulated  are  single  pre- 
cision numeric  quantities  (i.e.  no  logical,  alphanu- 
meric, or  double  precision  variables  and  quantities 
are  used  within  the  program  being  simulated).  Again, 
this  is  because  this  thesis  is  addressing  the  prob- 
lem of  numerical  precision  amd  the  stability  of 
mathematical  algorithms. 

6)  All  floating  point  values  to  be  n-bit  simulated  will 
have  three  characteristics  in  common  with  CDC  float- 
ing point  representations: 

a)  Mantissa  will  be  normalized  (i.e.  the  most  sig- 
nificant bit  of  the  mantissa  is  always  a one, 
except  in  the  case  of  zero). 

b)  Exponents  will  represent  powers  of  two,  which 
would  be  multiplied  with  the  mantissa  to  obtain 
the  floating  point  value  of  the  word,  for  in- 
stance, not  by  powers  of  16  for  hexadecimal  ma- 
chines. 

c)  Computers  being  simulated  represent  negative 
numbers  with  1*s  complement  representation, 
where  a negative  number  is  represented  by  a 1 
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in  the  sign  bit  and  the  1’s  complement  of  its 
magnitude  in  the  other  bit  positions.  Floating 
point  exponents  for  2*s  complement  representa- 
tion have  a positive  exponent  range  of  one  less 

••1  2R 

than  the  negative  exponent  (for  instance,  2 
to  2 whereas  1*s  complement  would  range 

from  2"''^'^  to 

General  Approach 

This  thesis  study  addressed  the  n-bit  simulation  problem 
at  two  levels.  Each  level  involved  designing  a preprocessor 
which  would  modify  the  algorithm's  code  so  that  the  numerical 
effects  of  the  n-bit  wordlength  could  be  realized.  The  first 
approach  attempted  to  accomplish  n-bit  simulation  through 
modification  of  the  COMPASS  assembly  language  representation 
produced  from  a FORTRAN  compilation  of  the  algorithm.  The 
COMPASS  code  would  be  modified  by  an  n-bit  simulation  pre- 
processor so  that  every  arithmetic  operation  would  be  replaced 
by  COMPASS  code  sequence  which,  when  executed,  would  generate 
the  exact  accuracy  attainable  on  a computer  of  the  n-bit 
wordlength  specified  by  the  user.  This  approach,  however, 
was  unsuccessful  in  meeting  the  requirements. 

The  second  level  approach  was  able  to  perform  n-bit 
wordlength  affects  upon  a FORTRAN  algorithm  by  modification 
of  the  FORTRAN  code.  This  modification  of  the  FORTRAN  code 
requires  a preprocessor  which  analyses  the  arithmetic  state- 
ments within  a program  and  replaces  them  with  subroutine 
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calls.  The  subroutines  perform  the  arithmetic  operations 
while  simulating  n-bit  wordlength  affects.  This  approach 
proved  successful  in  solving  the  n-bit  simulation  problem 
and  is  discussed  in  detail  in  Chapter  3. 

Chapter  Synopsis 

Chapter  2 of  this  thesis  report  contains  a summary  of 
the  analysis  and  design  accomplished  at  the  COM! ASS  level. 

A successful  prototype  of  the  n-bit  simulation  tool  could 
not  be  achieved  at  the  COMIASS  level.  The  reasons  for  the 
failure  are  elaborated  in  Chapter  2.  A working  prototype 
of  the  n-bit  simulator  was  achieved  by  way  of  the  FORTRAN 
level  approach.  This  approach  and  the  resulting  system  de- 
sign are  discussed  in  Chapter  3.  Chapter  4 briefly  summa- 
rizes the  effects  that  the  Structured  Design  techniques  had 
upon  the  system  design.  The  next  chapter  summarizes  the 
testing  performed  to  test,  debug,  and  validate  the  n-bit 
simulation  system.  Finally,  Chapter  6 draws  some  conclu- 
sions and  makes  recommendations  concerning  the  n-bit  simu- 
lation tool. 
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II.  COMPASS  Analysis  And  Design 


The  first  approach  to  solving  the  n-bit  simulation  prob- 
lem was  at  the  COMPASS  assembly  language  level.  This  ap- 
proach attempted  to  accomplish  the  numerical  effects  of  n-bit 
wordlength  upon  an  algorithm  through  modification  of  COMPASS 
code.  The  first  section  of  this  chapter  explains  why  the 
COMPASS  level  approach  was  undertaken  prior  to  the  FORTRAN 
level  approach.  The  following  sections  describe  an  ainalysis 
of  the  COMPASS  code  (produced  by  FORTRAN  compilation)  with 
respect  to  the  n-bit  simulation  problem.  Finally,  the  ap- 
proach taken  to  solve  the  n-bit  simulation  problem  and  the 
major  problems  encountered  which  resulted  in  the  failure  of 
the  COMPASS  level  approach  are  described. 

The  COMPASS  level  approach  was  anticipated  to  be  less 
complicated  than  the  FORTRAN  level  approach  to  the  n-bit  sim- 
ulation problem  for  a number  of  reasons.  In  a single  FORTRAN 
language"  statement  several  types  of  arithmetic  operations  and 
several  data  types  may  be  mixed  together  in  many  different 
ways.  Each  data  type  would  have  to  be  handled  differently, 
and  the  precedence  by  which  the  arithmetic  operations  are 
supposed  to  be  performed  would  require  further  analysis  and 
handling.  Also,  arithmetic  expressions  may  occur  in  almost 
any  type  of  FORTRAN  statement,  which  could  seriously  compli- 
cate n-bit  simulation  performance. 

On  the  other  hand,  at  the  CONTASS  level,  each  arithmet- 
ic operation  could  be  accomplished  by  a single  instruction 
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involving  only  equivalent  data  types.  The  arithmetic  preced- 
ence order  of  execution  and  conversion  of  data  types  has  al- 
ready been  accomplished  by  the  FORTRAN  compiler.  Furthermore, 
in  contrast  to  many  other  assembly  languages,  COMPASS  does 
not  include  base  addressing  and  relative  addressing.  This 
greatly  simplifies  the  insertion  of  additional  COMPASS  in- 
structions into  the  original  code,  which  would  enable  each 
arithmetic  instruction  to  be  performed  while  simulating  the 
n-bit  wordlength  effects. 

N-bit  Simulation  Accomplishment 

N-bit  simulation  for  an  entire  COMPASS  coded  algorithm 
can  be  achieved  by  simulating  the  effects  of  n-bit  wordlength 
for  each  arithmetic  instruction.  A typical  COMPASS  fixed 
point  arithmetic  instruction  is  1X7  X3+X4.  This  instruction 
adds  the  contents  of  register  X3  to  the  contents  of  register 
X4  and  stores  their  sum  into  register  X7.  To  simulate  the 
effects  of  n-bit  wordlength,  the  contents  of  registers  X3 
and  X4  would  be  modified  to  contain  the  equivalent  n-bit 
wordlength  values  (assuming  they  have  not  already).  Then  . 
the  two  registers  would  be  added  together  and  the  resulting 
contents  of  the  X7  register  would  be  similarly  converted  to 
the  n-bit  equivalent. 

Fixed  Point  N-bit  Simulation 

N-bit  simulation  is  accomplished  differently  for  fixed 
point  and  floating  point  number  representations.  The  range 
of  fixed  point  values  which  could  be  represented  is  limited 
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by  the  number  of  bits  in  the  computer  word.  An  n-bit  fixed 
point  word  would  retain  only  the  right-most  (least  signifi- 
cant) bits  of  the  CDC  60-bit  word.  For  example,  the  1 6-bit 
wordlength  equivalent  of  the  60-bit  value  shown  in  (a)  of 
Figure  2 is  shown  by  the  shaded  area  in  (a)  and  in  (b). 

N-bit  wordlength  representation  is  realized  for  fixed 
point  quamtities  by  comparing  the  60-bit  value  to  the  maxi- 
mum positive  Sind  maximum  negative  values  reprej^en table  with 
n-bits,  shown  respectively  in  (c)  and  (d)  for  16-bit  words. 

For  the  purposes  of  this  thesis,  whenever  a fixed  point  val- 
ue exceeds  one  of  these  raaximums,  an  overfow  condition  has 
occurred  and  the  overflow  value  is  replaced  with  the  maximum 
value  that  was  exceeded.  In  Figure  2,  (e)  contains  a nega- 
tive value  with  greater  thain  16  bits  of  significance.  The 
n-bit  simulator  would  react  by  replacing  it  with  the  value 
shown  in  (d) . 


Figure  2 Fixed  Point  N-bit  Simulation 
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Ploatinp:  Point  K-bit  Simulation 


To  n-bit  simulate  a floating  point  word  type,  the  cor- 
rect number  of  significant  bits  for  the  mantissa  must  be 
maintained  and  the  exponent  should  not  be  allowed  to  exceed 
the  equivalent  exponent  range  representable  in  an  n-bit 
floating  point  word.  CDC  floating  point  data  types  have  nor- 
malized mantissa.  Therefore,  in  order  to  maintain  the  spec- 
ified number  of  significant  mantissa  bits,  only  the  number  of 
bits  specified  for  it  would  be  retained  immediately  to  the 
right  of  the  least  significant  exponent  bit  (bit  #^2) . This 
can  be  accomplished  by  right  shifting  (with  sign  extension) 
the  contents  of  the  word  followed  by  a circular  left  shift. 
Right  shifting  with  sign  extension  displaces  the  contents  of 
a word  X positions  to  the  right,  where  X is  the  number  of 
bits  computed  from  subtracting  the  specified  number  of  man- 
tissa bits  from  48.  Forty-eight  is  the  length  of  the  CDC 
60-bit  floating  point  word  mantissa.  Those  bits  shifted  off 
the  right-end  are  lost.  All  positions  vacated  on  the  left 
by  the  right  shift  are  replaced  with  the  value  of  the  sign 
bit.  Figure  3 (a)  shows  an  example  of  a 60-bit  floating 
point  word  value  and  (b)  shows  the  resulting  contents  of  (a) 
after  a right-shift  of  40-bits  (for  an  8 bit  mantissa).  In 
a circular  left-shift,  bits  shifted  off  the  left  end  of  the 
word  replace  those  from  the  right  end.  In  this  manner,  the 
significant  bits  that  had  been  shifted  off  the  right  end  from 
the  original  word  are  replaced  with  zeroes  for  positive  num- 
bers and  ones  for  negative  numbers.  Therefore,  for  this 
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Floating  point  data  word  format 
12  12  13 


60 


8 bits  of 
original  mantissaJ 


(c)  After  circular  left  shift 


qooi 1 1 1 01 ooopi 01 01 1 1 0000000000000000000000000000000000000000, 


8 bit  mantissa 


Figure  3 Floating  Point  N-bit  Simulation 


example,  the  only  significance  for  the  mantissa  is  in  the 
left-most  8 bits  of  the  CDC  floating  point  mantissa.  Fig- 
ure 3 (c)  shows  the  resulting  contents  of  (b)  after  a circu- 
lar left  shift  of  40  bits.  Using  this  technique,  the  value 
shown  in  (a)  can  be  modified  to  contain  the  same  number  of 
significant  mantissa  bits  that  an  n-bit  floating  point  word 
with  an  8 bit  mantissa  would  contain,  this  equivalent  repre- 
sentation is  shown  in  (c). 

Some  significance  will  almost  always  be  lost  for  deci- 
mal number  representations,  since  decimal  numbers  can  not 
usually  be  exactly  represented  with  binary  storage.  For  ex- 
ample, the  60-bit  octal  representation  of  3.33  is  shown  in 
Figure  4 (a).  The  equivalent  representation  for  an  n-bit 
floating  point  word  with  a 15  bit  mantissa  is  shown  in  (b). 
The  right  most  33  bit,s  of  significance  have  been  lost  and 
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(a)  1721652525237366553 

(b)  172I652525OOOOOOOOO 


Figure  4 N-bit  Decimal  Representation 

the  value  left  is  equal  to  3.32971;  a significance  of  .00029 
has  been  lost  in  the  15-bit  mantissa  representation. 

Overflow  of  the  exponent  could  be  detected  by  comparing 
the  raaxiraum/rainimum  values  possible  for  a floating  point  word 
with  the  specified  bit  lengths  for  the  exponent  and  mantis- 
sa. Again  overflow  detection  would  result  in  the  maximum/ 
minimum  value  being  stored  over  the  old  value.  In  addition 
to  the  case  where  the  exponent  gets  too  large,  the  case  of 
exponent  underflow  must  also  be  addressed.  This  could  occur 
when  a value  is  computed  which  is  so  small  that  its  most  sig- 
nificant digit  is  not  within  the  range  of  the  specified  ex- 
ponent (although  it  is  within  the  CDC's  large  range).  For 
example,  5 bits  may  be  specified  for  the  n-bit  floating  point 
exponent  and  6 bits  for  the  mantissa.  This  exponent  could 
represent  an  exponent  range  of  ll6.  Therefore,  the  maximum 
values  representable  with  that  floating  point  word  would  be 
163  (the  maximum  mantissa  value)  times  2*^^  if  the  implied 
decimal  point  is  to  the  right  of  the  least  significant  man- 
tissa bit. 

If  the  implied  decimal  point  is  to  the  left  of  the  main- 
tissa,  the  maximum  would  be  .984375  times  2"*’^^  (shown  in  Fig- 
ure 5) . The  minimum  significance  representable  would  be  t'^2 
times  2“  and  1.5  times  2 respectively  for  the  two 
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different  decimal  point  positions.  Underflow  for  the  12-bit 

floating  point  word  would  occur  if  a 60-bit  division  resulted 

-24 

in  a quotient  less  than  2 . This  value  would  not  be  with- 

in the  range  representable  by  a word  with  a 5 bit  exponent. 
Zero  would  replace  the  value  that  overflowed  before  the  pro- 
gram executed  its  next  operator. 

011111*111111  = .984375  * 2*^^^ 

011111111111^  = 63  * 2'^”’^ 

000000^100000  = .5  * 2“''^ 

000000100000^=  32  * 2"^^ 

Figure  5 N-bit  Floating  Point  Value  Range  for  n Equals  12 

Now  that  the  strategy  has  been  defined  for  simulating 
the  effects  of  n-bit  wordlength  for  arithmetic  instructions 
and  the  fixed  or  floating  point  word  values  involved,  the 
following  section  will  describe  the  arithmetic  instructions 
which  the  n-bit  simulation  would  be  performed  upon. 

The  Arithmetic  Instructions 

Many  COMPASS  code  listings,  resulting  from  the  three 
FORTRAN  optimizing  compilers,  were  examined  to  investigate 
the  handling  of  FORTRAN  assignment  statements  and  arithmetic 
expressions.  The  arithmetic  operations  addition,  subtrac- 
tion, multiplication,  and  division  were  found  to  be  handled 
by  the  COMPASS  instructions  in  Table  I,  where  Xi,  Xj,  and 
Xk  may  be  any  combination  of  the  eight  X registers  (XI  - X8) . 
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COMPASS  Instruction 


Operation  type 


fixed  pt 
floating 
fixed  pt 
floating 
fixed  pt 
floating 
fixed  pt 
floating 


pt  add 
subtract 
pt  subtract 
multiply 
pt  multiply 
divide 
pt  divide 


Octal  code  for 
operation  type 


♦fixed  point  values  in  registers  Xj  and  Xk  are 
converted  to  their  floating  point  equivalents 
before  the  division  is  performed.  The  quo- 
tient Xi  is  then  converted  back  to  a fixed 
point  representation. 


Table  I Fixed  And  Floating  Point  Arithmetic  Instructions 


Each  fixed  and  floating  point  single  precision  arithmetic 
operation  generated  by  the  FORTRAN  compiler  was  performed  by 
a unique  instruction  with  the  exception  of  division.  As  not- 
ed in  Table  I,  fixed  point  division  is  performed  by  a sequence 
of  instructions.  The  divisor  and  dividend  are  converted  to 
their  floating  point  equivalents  before  a floating  point  di- 
vide is  performed.  The  quotient  then  is  converted  back  to  its 
fixed  point  equivalent  representation.  This  sequence  of  in- 
structions for  the  fixed  point  quotient  X1=X2/X5  is  shown  be- 


low; 


PX2  B0,X2  [Pack  X23  ^ „ . . vo  o ^ + 

PX3  B0,X3  [pack  X3J  I ^ 

NXO  B0,X3  formalize  X3]J  floating  point 

FX1  X2/X0  divide] 

TY1  nv’yi  Converting  result  back 

LX1  B7,X1  (Shift  to  integer?  ^ 

position]  / ^ 


Since  fixed  point  and  floating  point  divisions  share  a 
common  instruction  (FXi  Xj/Xk)  and  each  type  requires  a dif- 
ferent method  to  simulate  the  effects  of  n-bit  wordlength,  a 
choice  would  have  to  be  made  not  to  allow  one  or  the  other  of 
the  division  types  to  occur  at  the  FORTRAN  program  level. 
Floating  point  was  expected  to  be  the  more  common  usage  for 
algorithms,  so  fixed  point  divisions  would  be  restricted  from 
n-bit  simulated  algorithms.  If  the  fixed  point  division  se- 
quence always  appeared  as  consecutive  instructions,  fixed 
and  floating  point  division  could  be  differentiated;  but  the 
FORTRAN  optimizing  compilers  order  the  COWIASS  instructions 
for  the  most  efficient  execution,  and  consequently  the  COM- 
PASS code  generated  for  a sequence  of  FORTRAN  statements  may 

bd"  completely  intermixed. 

* 

One-  other  potential  problem  could  occur  for  the  fixed 
point  multiplication  instruction  DXi  Xj*Xk  shown  in  Table  I. 
DXi  Xj*Xk  is  actually  the  COMPASS  double  precision  multiply 
instruction.  However,  this  conflict  of  usage  would  not  oc- 
cur since  the  algorithms  for  the  n-bit  simulator  would  not 
require  double  precision  FORTRAN  variables  (Assumption  5 
states  this  in  Chapter  I) . 

Now  that  the  COMIASS  instructions,  which  perform  the 
arithmetic  computations  have  been  described,  two  methods  of 
incorporating  modifications  for  each  arithmetic  instruction 
will  be  described  in  the  following  section. 
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I Two  Approaches  to  Incorporating  N-bit  Simulation 

[ Two  fundamental  approaches  to  implementation  of  the  n- 

i bit  simulation  modifications  were  considered.  The  first  ap- 

I proach  would  replace  each  arithmetic  instruction  with  subrou- 

tine calls.  The  second  approach  would  substitute,  in-line, 
the  necessary  instructions  to  perform  the  n-bit  simulation 
effects . 

The  subroutine  approach  would  require  an  argument  list 
to  accompany  each  subroutine  call  indicating  the  registers 
involved  with  that  instruction.  This  argument  list  would 
key  the  handling  subroutine  to  manipulate  the  designated  reg- 
isters. Uniquely  identifying  each  instruction/register  com- 
bination with  a separate  subroutine  would  require  512  dif- 
ferent subroutines  per  arithmetic  instruction  type  (three 
registers  are  used  for  each  arithmetic  instruction  and  each 
register  may  be  any  one  of  eight  possible  X registers,  mak- 
ing 8^  or  512  possible  combinations).  Clearly  that  would 
have  been  too  many  routines.  So  an  argument  list  would  be 
necessary  indicating  the  registers  involved,  and  a technique 
would  have  to  be  devised  in  each  subroutine  to  handle  them 
correctly. 

In  contrast,  the  in-line  code  approach  would  be  simpler 
and  save  subroutine  linkage  time,  but  would  require  more  core 
space.  If  branching  would  be  necessary  within  the  substi- 
tuted code  sequence,  some  method  would  have  to  be  devised  to 
provide  unique  labels  within  each  substituted  sequence.  For 
, the  initial  COMPASS  level  investigation,  the  in-line  approach 

was  chosen. 

1' 

) 

i 
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Methods  for  Isolating  Arithmetic  Instructions 

In  both  cases  previously  stated,  the  first  obstacle  was 
trying  to  isolate  arithmetic  operation  instructions  of  Table 
I from  the  rest  of  the  COMPASS  code.  Three  techniques  were 
considered . 

The  first  technique  involved  exaimining  the  raw  COMPASS 
object  code  (machine  language).  This  would  be  a complicated 
procedure  due  to  problems  involved  distinguishing  data  from 
instructions.  Also  labels,  already  assigned  to  fixed  loca- 
tions, are  transparent  from  their  associated  object  code. 
Injecting  additional  COMIASS  object  code  to  perform  n-bit 
simulation  would  lead  to  misalignment  of  labels  with  respect 
to  their  associated  object  code  locations.  When  examining 
the  object  code,  a table  containing  all  possible  COMPASS  in- 
struction codes  and  the  associated  instruction  length  would 
have  to  be  maintained.  In  addition,  one  must  account  for 
no-operation  instructions  (NOOIS) . NOOPS  are  used  to  fill 
in  remaining  portions  of  60-bit  words  which  are  not  com- 
pletely filled  by  instructions  due  to  labels  or  other  con- 
ditions. It  would  be  a cumbersome  task  indeed.  Consequent- 
ly, this  technique  was  quickly  ruled  out. 

The  second  technique  involved  examining  the  COMPASS 
source  code  generated  from  the  FORTRAN  compilation  of  the 
algorithm.  Distinguishing  the  arithmetic  operation  instruc- 
tions would  require  a table  lookup  each  COMPASS  code  mne- 
monic. When  matches  occurred,  a sequence  of  instructions 
would  be  inserted  to  perform  the  n-bit  wordlength  effects 
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for  that  instruction  type.  A significant  advantage  of  this 
method  is  that  the  COMPASS  assembler  would  do  the  work  of 
realigning  all  COMPASS  object  code  and  label  pointers.  COM- 
PASS source  code  would  be  accessed  by  using  the  E option  on 
the  FTN  control  card.  The  sequence  of  control  cards  in  Fig- 
ure 6 could  be  used  to  access  the  COMPASS  source  code  gener- 
ated from  a FORTRAN  compilation.  The  execution  of  NBIT  (in 
Figure  6)  would  modify  the  original  COMPASS  source  code  (the 
COMP  file)  for  n-bit  simulation  producing  the  revised  n-bit 
simulated  code  in  file  C0MP2.  C0MP2  would  then  be  assembled 
and  executed. 

FTN,  E = COMP. 

NBIT,  COMP,,,  C0MP2. 

REWIND,  C0MP2. 

COMPASS,  I = C0MP2,  S = FTNMAC,  B » BIN. 

LDSET  (LIB  a FORTRAN/SYSIO) . 

BIN. 

— - — 1 — i — 

Figure  6 N-bit  Simulation  Control  Card  Sequence 

The  third  method  for  isolating  and  handling  the  arithmet- 
ic instructions  would  take  advantage  of  some  sophisticated 
features  of  the  COMPASS  assembler  called  operation  definitions 
(OPDEFs)  and  IP  conditionals  (IFCs).  Both  are  pseudo  instruc- 
tion types.  The  OPDEP  consists  of  a sequence  of  COMPASS 
source  code  that  is  saved  and  assembled  when  an  instruction 
type  within  a given  routine  matches  the  syntax  specification. 
It  usually  includes  parameters  to  be  substituted  for  formal 
parameters,  so  that  the  object  code  generated  can  vary  with 
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each  assembly  of  the  definition  and  whatever  combination  of 
registers  that  are  involved  in  the  instruction  (REP  3:  5-27). 
Using  this  feature,  the  assembler  would  do  the  work  by  re- 
placing the  arithmetic  instructions  with  the  n-bit  simulation 
sequence . 

Figure  7 shows  a COMPASS  routine  in  which  ein  OPDEF  was 
inserted,  then  assembled  with  the  rest  of  the  COMPASS  code. 
This  particular  OPDEF  syntactically  defines  the  general  in- 
struction type  IXi  Xj+Xk.  PI,  P2,  and  P3  used  within  the 
OPDEF  would  represent  the  registers  Xi,  Xj,  and  Xk  respec- 
tively. In  the  sequence  of  original  COMASS  code,  the 
1X7  XO+X4  instruction  would  match  this  syntax  and  is  re- 
placed by  this  sequence  when  it  is  assembled  as  shown  in  the 
after  assembly  code  sequence.  After  assembly,  the  object 
code  (shown  on  the  left)  is  produced  by  the  COMPASS  assem- 
bler upon  detection  of  the  1X7  X0+X4  instruction.  It  can 
be  seen  after  the  assembly  that  the  OPDEF  for  the  fixed 
point  add  instruction  replaced  the  original  one. 

Two  additional  pseudo  instructions  were  also  used  in 
this  example:  SET  (REF  3:  4-36)  and  CPOP  (REF  3:  6-7). 

CPOP  was  used  to  disguise  the  fixed  point  addition  instruc- 
tion from  its  own  OPDEF.  Otherwise,  if  am  instruction  were 
used  within  its  own  OPDEF  definition  it  would  result  in  the 
assembler  getting  caught  in  an  infinite  loop  of  replacing 
an  instruction  with  itself,  which  triggers  the  process  a- 
gain,  and  so  on. 
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loca-  Code  Generated 

Label  Opera- 

Variable 

tion 

tion 

(Before  assembly)  ^ 

KBITS 

IDENT 

SET 

TEST 

16 

f 

IXX+X 

OPDEF 

PI ,P2,P3 

LX.P2 

59-NBITS 

Inserted  into  \ 

LX.P3 

AX.P2 

59-NBITS 

59-NBITS 

original  / 

COMPASS  code  A 

XAA+A 

AX.P3 

CPOI 

59-KBITS 

0,360B,132B 

/ 

XA,P1 

A.P1+A.P3 

/ 

LX. PI 

59-KBITS 

i 

AX  .PI 
ENDM 

59-KBITS 

SA5 

AO 

SX7 

116100B 

Original  COM-  1 

1X2 

X3-X4 

PASS  code  from  < 

1X7 

X0+X4 

FORTRAN  compi-  / 

DX7 

X4*X4 

lation  V 

END 

(After  assembly) 

IDENT 

TEST 

KBITS 

SET 

16 

IXX+X 

OPDEF 

PI  ,P2,P3 

LX.P2 

59-NBITS 

LX.P3 

59-NEITS 

AX.P2 

59-KBITS 

AX.P3 

59-NBITS 

XAA+A 

CPOP 

0,360B,132B 

XA.P1 

A.P1+A.P3 

LX. PI 

59-KBITS 

AX  .PI 
ENDM 

59-NBITS 

0 54500 

SA5 

AO 

7170116100 

SX7 

116100B 

36234 

1X2 

X3-X4 

1X7 

X0+X4  ^ 

1 20053 

LX.O 

59 -KB  ITS  ^ 

) 

20453 

LX. 4 

59-KBITS  i 

Substituted 

21053 

AX.O 

59-NBITS  / 

by 

21453 

AX. 4 

59-NBITS  > 

assembler 

XAA+A 

CPOP 

0,360E,132B  \ 

for 

2 36704 

XA.7 

A.0+A.4 

1X7  X0+X4 

20753 

LX. 7 

59-NBITS  i 

21753 

AX. 7 

59-NBITS  y 

42744 

DX7 

X4*X4 

END 

Figure  7 OIDEF  Application 
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The  other  pseudo  instruction  mentioned  eeirlier,  IFC, 
could  be  used  to  incorporate  the  user  options,  for  example, 
truncation  or  rounding  effects.  Each  effect  would  require 
somewhat  different  code.  Both  code  sequences  could  be  in- 
corporated into  the  arithmetic  instruction  OIDEF,  but  only 
the  one  whose  option  was  on  would  be  invoked  as  part  of  the 
OPDEF,  upon  assembly. 

This  package  of  COMPASS  pseudo  instructions  could  great- 
ly simplify  the  n-bit  simulation  set  up.  However,  a means 
would  have  to  be  devised  to  get  the  OPDEFs  inserted  within 
each  routine,  since  they  are  effective  only  for  the  routine 
in  which  they  were  defined . If  brauiching  were  required  with- 
in an  OPDEF  to  perform  n-bit  simulation  effects,  a special 
element  called  *L  (REF  3:  2-9)  and  the  SET  pseudo  could  be 
used  within  each  OPDEF  to  offer  each  invocation  of  the  OPDEF 
a unique  label.  Figure  8 shows  an  example  of  its  usage  with- 
in an  OPDEF.  The  SET  *L  + 2 instruction  would  put  the  cor- 
rect location  into  the  symbol  table  for  N (shown  at  location 
15).  When  ZR  X.P1,N  referenced  it,  the  value  (a  location) 
currently  in  the  symbol  table  would  be  put  in  the  object  code. 
The  SA.P1  P,1  instruction  (where  the  label  is  actually  needed) 
would  be  positioned  a known  number  of  locations  down  from  the 


Location 

Label 

Instruction 

15 

N 

SET  *L  + 2 

ZR  X.P1,N  (wants  to  jump  to 

SA.P1  P,1) 

16 

(other  instructions) 

17 

SA.P1  P,1  [heeds  the  label] 

Figure  8 Example  Of  Labeling  Within  An  OPDEF 
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SET  instruction.  Since  each  Oi^DEF  invocation  would  occur  at 
a different  location,  each  reference  would  be  unique,  since 
the  value  of  N would  cheuige  with  the  location  of  each  invo- 
cation. 

Problems  Encountered  With  The  COMPASS  Level  Approach 

After  further  analysis  of  the  COMPASS  level  approach, 
three  problems  were  uncovered.  These  problems  will  be  dis- 
cussed in  the  following  subsections. 

Multiples  Of  Two 

Upon  further  analysis,  it  was  discovered  that  when  the 
FORTRAN  optimizing  compiler  1 or  2 was  used,  fixed  point  mul- 
tiplications and  divisions  by  multiples  of  two  were  not  con- 
verted into  arithmetic  multiply  or  divide  instructions,  but 
were  accomplished  by  shifting  the  register  contents  left  or 
right  the  number  of  bits  required  to  obtain  the  proper  re- 
sult. However,  the  zero  level  compiler  did  convert  such  oc- 
currences into  the  expected  arithmetic  instructions.  There- 
fore, a user  would  either  have  to  compile  the  algorithm  using 
the  level  0 optimizer  which  would  result  in  less  efficient 
execution  or  be  restricted  from  using  multiples  of  two  in  ar- 
ithmetic equations,  since  those  arithmetic  operations  could 
not  be  n-bit  simulated.  The  user  could  still  use  a multiple 
of  two  in  an  equation  for  optimization  level  1,  but  not  as 
a constant  within  the  equation.  It  could  overcome  the  op- 
timizers intelligence  by  assigning  that  multiple  of  two  (con- 
stant) to  a variable  just  prior  to  the  expression  as  shown 
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below.  This  procedure  would  force  the  compiler  to  generate 
the  COMPASS  arithmetic  instruction. 


(desired  equation:  IRSLT  = 11/8  * 12) 

I TEMP  = 8 

IRSLT  = II  / ITEFJ  * 12 

The  problem  with  multiples  of  two  could  be  overcome  by  put- 
ting an  added  restriction  upon  the  user.  However,  the  prob- 
lems described  in  the  following  subsections  are  of  a more 
serious  nature  and  could  not  be  resolved. 

No  Double  Precision  Arithmetic  Accuracy 

Another  problem  that  became  evident  quickly  was  that  no 
FORTRAN  arithmetic  expression  could  be  evaluated  in  double 
or  triple  precision,  then  stored  as  a single  precision  value 
This  is  because  at  the  COMPASS  level  the  last  arithmetic  op- 
eration of  a compiled  FORTRAN  equation  is  indistinguishable 
from  the  first.  There  is  no  way  at  the  COMPASS  level  to  de- 
tect the  first  operation  generated  from  a FORTRAN  arithmetic 
statement.  The  COMPASS  code  generated  from  individual  FOR- 
TRAN statements  is  transparent  at  the  COMI-ASS  level  and  in 
many  cases  the  code  from  several  FORTRAN  statements  are  in- 
termixed together  so  the  computer  can  execute  with  greater 
efficiency.  Thus  all  operations  would  be  limited  to  single 
precision  auid  that  user  option  specified  in  the  requirements 
would  not  be  available. 


that  extra  registers  would  be  needed,  in  addition  to  the  ones 
involved  in  the  arithmetic  instruction,  in  order  to  perform 
the  required  n-bit  simulation  effects.  These  extra  registers 
would  be  required  for  storing  the  maximum  and  minimum  values 
of  the  n-bit  fixed  and  floating  point  words.  To  determine 
overflow  conditions  these  values  would  have  to  be  compared 
to  the  contents  of  the  registers  involved  in  the  arithmetic 
instruction.  In  order  to  make  the  comparisons,  the  values 
must  be  stored  in  a 60-bit  register.  However,  due  to  the 
nature  of  the  COMPASS  level  approach,  nothing  can  be  presumed 
to  be  known  about  any  register  before  that  instruction  or 
after.  At  the  time  the  arithmetic  instruction  is  executed, 
any  or  all  of  the  16  X and  B registers  may  contain  values 
which  were  set  up  by  the  FORTRAN  compiler  to  be  used  later. 
Modification  of  a register's  contents  could  be  disasterous. 

Since  the  status  of  all  registers  outside  the  ones  in- 
volved in  the  arithmetic  instruction  is  not  known,  all  reg- 
isters not  involved  in  the  arithmetic  operation  instruction 
must  have  their  values  restored  to  what  they  were  Just  prior 
to  the  operation.  Either  a register  saving/restoring  tech- 
nique would  have  to  be  employed  or  else  some  register  with 
known  or  fixed  values  would  have  to  be  found  which  could  be 
used  as  an  extra  register  then  restored  afterwards. 

Several  possibilities  were  investigated  and  some  tech- 
niques were  devised,  but  none  which  were  reliable  10096  of 
the  time.  The  different  avenues  taken  are  described  more 
fully  in  Appendix  A.  Literature  searches  in  CDC  manuals 
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could  not  produce  any  amount  of  detail  about  the  register 
handling  conventions  of  the  FORTRAN  compilers.  Nor  could  any 
register  saving/restoring  convention  be  found.  Apparently, 
the  compilers  do  not  save  all  registers  before  jumping  to  ex- 
ternal subroutines,  as  many  computer  systems  do  (e.g.  IBM). 

Instead  it  saves  only  those  registers  it  knows  through  its 
intelligence  that  it  will  need  later. 

Without  extra  registers,  no  rounding  effects  could  be 
performed  and  no  overflows  could  be  detected.  It  logically 
follows  that  without  overflow  detection,  no  overflow  mes- 
sages could  be  printed  nor  would  there  be  a way  or  means  to 
substitute  meiximum  values  for  values  that  overflow.  The  on- 
ly action  that  could  be  performed  would  be  to  limit  fixed 
point  values  amd  floating  point  mantissas  to  the  specified 
number  of  bits  of  precision  (through  left  and  right  register 
shifts,  to  accomplish  the  truncation). 

Conclusions 

The  more  the  COMtASS  level  approach  was  explored,  the 
more  limitations  and  shortcomings  were  uncovered.  Single 
assignment  statements  could  not  be  detected,  double  and  tri- 
ple precision  operations  could  not  be  performed  as  options, 
and  no  completely  reliable  scheme  could  be  devised  to  gener- 
ate extra  registers  required  for  overflow  detection  and  han- 
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dling.  Therefore,  the  COMPASS  level  approach  was  aborted. 

The  n-bit  simulation  problem  was  next  addressed  at  the  FOR- 
TRAN level  as  discussed  in  the  following  chapter. 
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III.  FORTRAN  Analysis  And  Design 


After  failing  to  establish  an  operational  n-bit  simula- 
tion at  the  COMPASS  level,  a FORTRAN  level  approach  to  the 
n-bit  simulation  problem  was  endeavored.  The  FORTRAN  level 
approach  attempted  to  accomplish  n-bit  simulation  by  modify- 
ing the  user's  FORTRAN-programed  algorithms.  This  chapter 
will  describe  an  einalysis  of  the  FORTRAN  language  with  re- 
spect to  the  n-bit  simulation  problem.  Also  it  will  describe 
the  features  of  the  FORTRAN  algorithm  preprocessor  and  n-bit 
simulation  subroutines  that  were  designed  to  accomplish  the 
n-bit  simulation. 

FORTRAN  Problem  Analysis 

A basic  assumption  for  this  thesis  was  that  n-bit  word- 
length  effects  could  essentially  be  simulated  by  modifying 
only  the  arithmetic  expressions  within  a program.  This  is 
because  the  n-bit  simulation  tool  was  devised  to  enable  a 
user  to  evaluate  the  numerical  stability  and  precision  of 
algorithms  as  a function  of  wordlength.  In  Standard  ANSI 
FORTRAN,  arithmetic  expressions  may  appear  in  almost  all 
statement  types.  N-bit  simulation  of  arithmetic  expressions 
for  all  possible  statement  types  would  complicate  the  FOR- 
TRAN level  approach  to  a considerable  extent.  Therefore,  it 
was  decided  that  the  initial  FORTRAN  n-bit  simulation  tool 
developed  would  limit  n-bit  wordlength  effects  to  only  those 
expressions  occurring  with  FORTRAN  assignment  statements. 
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This  would  be  the  minimum  requirement  to  n-bit  simulate  an 
entire  algorithm  and  would  still  not  unduly  inhibit  the  user 
from  programing  in  his  normal  manner. 

To  simulate  the  effects  of  n-bit  wordlength  for  an  ar- 
ithmetic statement,  each  arithmetic  operation  and  its  asso- 
ciated operands  would,  first,  have  to  be  isolated  from  the 
rest  of  the  statement.  Then  each  operand  would  have  to  be 
modified  so  that  it  represents  the  n-bit  wordlength  equiva- 
lent (if  it  has  not  already).  Next,  after  performing  the 
operation,  the  result  must  be  similarly  modified.  If  the 
operation  would  have  resulted  in  an  overflow  condition  for 
the  n-bit  computer,  an  overflow  message  should  be  printed  to 
the  user  and  the  overflow  value  should  be  replaced  with  the 
maximum  value  representable  with  an  n-bit  word  (as  specified 
in  the  n-bit  simulation  requirements  of  Chapter  I). 

To  correctly  evaluate  an  arithmetic  statement,  the  or- 
der of  precedence  in  which  the  FORTRAN  arithmetic  operations 
are  normally  performed  would  have  to  be  preserved.  Also  the 
variable  types  (fixed  and  floating  point)  would  have  to  be 
determined  then  handled  appropriately  and  labels  accompany- 
ing arithmetic  statements  would  have  to  be  handled  so  that 
the  execution  sequences  of  the  algorithm  would  not  be  af- 
fected . 

Figure  9 shows  the  Backus  Normal  Form  (BNF)  for  all 
possible  syntactic  representations  of  an  arithmetic  assign- 
ment statement.  Any  arithmetic  assignment  statement  may 
consist  of  unaries,  integer  (fixed  point)  and  real  (floating 
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<arithmetic  statement> 


<identifier>  = <expression> 


<identif ier> 

* * B 

<letter>  | <identifier>  <letter> 
<identifier>  <digit>  | <function> 

<expression> 

• • M 

<sign>  <terni>  | <expression>  csign'> 
<terra> 

<term> 

cfactor>  1 <term>  <raul-op>  <factor> 

<factor> 

• • V 

• • •“ 

<f2>  1 <factor>  <exponentiation> 

<f2> 

Cf2> 

• • V 

<identifier>  | <literal>  ) 

( <expression>  ) 

<function> 

<identifier>  ( ^arguments)  ) 

<arguments> 

<expression>  ( <arguments>  , 
<expression> 

<literal> 

: : = 

<12  > 1 <real  #> 

<real 

: ; = 

<12>  . 1 <I2>  . <I2>  I . <I2> 

<12  > 

: : = 

<digit>  1 <I2>  <digit> 

<letter> 

• • ss 

A|B|c|DjE|FjG|H| l| J|  k|  L|M|N|o|P( 
Q|R| S|T|U1V| W|X|Yl Z 

<digit> 

: : = 

0|  i|2|3|4|5|6|7|8|9 

<sign> 

• • S 

<sign'>  j <null> 

<sign'> 

2 S s 

<raul“Op> 

•1/ 

<exponentiation> 

2 2s 

** 

Figure  9 Backus  Normal  Form  For  Arithmetic  Statement 
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point)  variables,  integer  and  real  literals,  arithmetic  op- 
erations, parentheses,  array  items,  and  function  subroutine 
calls.  Some  examples  of  various  arithmetic  assignment  state- 
ments are  shown  in  Figure  10.  Boolean  IF  statements  may  al- 
so contain  arithmetic  assignment  statements.  The  arithmetic 
statement  would  be  executed  if  the  logical  or  relational  ex- 
pression in  the  IF  statement  predicate  is  true.  The  basic 
form  of  a Boolean  IF  statement  with  an  arithmetic  statement 
is: 

IF  (Logical  or  relational  expression)  <arithmetic  stmt> 


Label 

Col. 

field 

7 

A = B*C^^(L0NG1-5.324/(-24)) 

1 

B = FUNC(4*J,--K,F2-LAST)+NUMB1 

C = -W**5 

25 

D = ARRAY1 (5, I)*5. 014327 

E = +(o.5**3**TWO)/DIVDER 

4 

IF(A.EQ.(10*R))F  = 2000+45.+.3*REALNUM 

RESULT  = FINAL 

Figure  10  Examples  Of  Legal  Arithmetic  Statements 

Evaluation  Of  ^ Arithmetic  Statement  Through  Reverse  Polish 
Notation 

The  task  of  simulating  the  effects  of  an  n-bit  wordlength 
machine  requires  that  the  result  of  a computation  for  an  ar- 
ithmetic statement  be  modified  to  reflect  its  n-bit  wordlength 
equivalent.  This  necessitates  the  isolation  of  each  arith- 
metic operation  and  the  operands  the  operation  is  acting  upon. 
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Therefore,  each  FORTRAN  arithmetic  statement  must  be  broken 
down  into  a series  of  arithmetic  statements,  each  performing 
a single  arithmetic  operation.  The  final  result  from  the  se- 
ries of  single  computation  statements  would  be  algebraically 
equal  to  the  result  computed  from  the  whole,  complex  arith- 
metic statement  from  which  they  were  derived. 

An  example  of  breaking  a complex  statement  into  a series 
of  single  computation  statements  is  shown  in  Figure  11.  E- 
quation  (a)  of  Figure  11  could  be  accomplished  just  as  well 
through  the  sequence  of  single  operation  arithmetic  state- 
ments shown  in  (b).  However,  the  intermediate  variables  used 
in  (b)  must  agree  with  the  dominant  data  type  (mode)  of  each 
computation.  That  is,  when  an  expression  contains  both  fixed 
8uid  floating  point  quantities,  the  final  result  would  have 
to  be  of  the  floating  point  data  type  since  floating  point 
is  dominant  over  fixed  point  in  FORTRAN. 


(a)  RESULT  = A*B*(D( I)-F/4 . 5) +2 .**EXrNT 


Precedence 
Order  of 
Evaluation 


(b)  TEMPI  = A*B 
TEMP2  = F/4.5 
TEMP  3 = D(I)-TEMP2 
TEMP4  - TEMP1*TEMP3 
TEMP 5 = 2.**EXPNT 
RESULT  - TEMI4+TEMP5 


Figure  11  Equivalent  Computation  Of  A Complex 
Arithmetic  Equation 
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When  a FORTRAN  compiler  analyzes  FORTRAN  arithmetic 


statements,  it  determines  the  proper  sequence  in  which  the 
computations  should  be  performed  and  the  necessary  mode  con- 
versions needed  for  the  differing  variable  data  types.  Then 
the  compiler  generates  the  assembly  language  instructions 
which  will  do  the  job;  forming  intermediate  results  between 
computations  until  the  final  assignment  is  made. 

An  approach  similar  to  one  a FORTRAN  compiler  might  use 
for  evaluating  an  arithmetic  statement  is  the  Early  Operator 
Reverse  Polish  Notation  (REF  5:  55).  Polish  notation  is  a 
mathematical  notation  that  provides  for  the  representation 
of  complex  expressions  in  a non-ambiguous  manner,  without  re- 
lying on  hierarchial  delimiters  such  as  parenthesis.  The 
basic  difference  between  standard  mathematical  notation  and 
Polish  notation  is  the  relative  position  of  the  operator  to 
its  operands. 

Figure  12  shows  an  example  of  a FORTRAN  arithmetic 
statement  and  its  equivalent  Reverse  Polish  Notation  repre- 
sentation. To  evaluate  the  Polish  string,  it  is  scanned 
left  to  right  until  an  operator  is  found.  Then  that  opera- 
tor is  used  on  the  preceding  two  operands-replacing  the  two 


Standard  mathematical  notation  for  a typical  arith.  equation: 

A = B+C-D*(E+F) 

Reverse  Polish  Notation  equivalent  of  equation  above; 

A E F+*-=  ^ 

Order  of  Ua— J 2 E+F 

Evaluation  i ^ i 3 D*(E+F) 

'I'  . . 4 (B+C)-(D*(E+F)) 

* — ^ ' 5 A=((B+C-(D*(E+F))) 


Figure  12  Example  Of  Reverse  i-olish  Notation 


by  the  result  auid  continuing  the  scan.  Each  arithmetic  oper- 
ator has  a hieraurchial  precedence  value.  Table  II  shows  the 
precedence  relationship  among  the  arithmetic  operations  of 
the  FORTRAN  language.  Operations  with  higher  precedence  in 
the  hierarchy  are  performed  first  in  the  evaluation  of  am  ar- 
ithmetic equation.  Operations  of  equal  precedence  are  eval- 
uated left  to  right.  The  BNF  for  FORTRAN  arithmetic  state- 
ments in  Figure  9 also  reflects  this  precedence  of  the  oper- 
ators . 


Arithmetic  Precedence 

operation  rank 


+ 


Higher 

precedence 


* 

/ 

- (unary) 
+ (unary) 

** 


1 

1 

2 

2 

3 

3 

4 


Table  II  Operator  Precedence  Order  For  Evaluating 
Arithmetic  Equations 

Through  the  use  of  Polish  Notation  techniques,  the  proc- 
ess of  generating  a sequence  of  single  operator  statements 
from  a single  complex  arithmetic  statement  is  greatly  simpli- 
fied. Once  an  arithmetic  statement  is  in  its  equivalent  Pol- 
ish Notation  form,  the  Polish  string  can  be  scanned  left  to 
right  for  operators.  For  each  operator  a function  subroutine 
call  is  generated  with  the  operands  involved  passed  along  as 
arguments.  The  subroutine  called  would  perform  the  computa- 
tion while  simulating  the  effects  of  n-bit  wordlength. 
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Additional  FORTRAN  Language  Cons iderat ions/1 roblems 

In  the  analysis  of  the  FORTRAN  language  with  respect  to 
the  n-bit  simulation  problem,  several  obstacles  and  potential 
problems  had  to  be  dealt  with  to  handle  all  cases  of  arith- 
metic statements.  Besides  the  variables  amd  arithmetic  op- 
erations, an  arithmetic  statement  may  involve  function  sub- 
program references,  arrays,  labels,  and  IP  statements.  An 
additional  consideration  is  to  provide  the  user  a means  of 
tracing  down  overflow  occurrences.  These  cases  will  be  dis- 
cussed in  the  following  subsections. 

Function  References.  Function  subprogram  references 
with  a FORTRAN  arithmetic  statement  present  a potential  prob- 
lem since  the  functions  they  reference  may  not  be  coded  in 
FORTRAN.  To  be  n-bit  simulated  completely,  all  code  executed 
must  be  written  in  the  FORTRAN  language.  System  functions 
such  as  SIN,  COS,  and  ATAN  (external  or  FORTRAN  library  func- 
tions (REP  4,  1-8))  are  not  coded  in  the  FORTRAN  source  lan- 
guage and  therefore  can  not  be  altered  by  the  FORTRAN  level 
approach  to  simulate  n-bit  wcrdlength  effects.  A conversion 
of  all  intrinsic  and  external  functions  into  their  FORTRAN 
language  equivalents  would  remedy  this  problem.  However, 
without  a FORTRAN  version,  the  most  effective  action  that 
can  be  done  to  approximate  n-bit  wordlength  affects  on  an 
intrinsic  function  would  be  to  n-bit  simulate  each  argument 
of  the  function  and  then  n-bit  simulate  the  function  value 
returned . 
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Arrays . N-bit  simulation  for  each  argument  (subscript) 
of  a FORTRAN  array  item  would  not  really  serve  any  purpose. 
FORTRAN  array  arguments  must  always  be  fixed  point  values. 
They  are  usually  relatively  small  fixed  point  values  and 
should  be  well  within  the  range  of  an  n-bit  wordlength  com- 
puter's fixed  point  words,  in  presumably  every  case.  Con- 
sequently, it  was  considered  unnecessary  overhead  to  n-bit 
simulate  expressions  and  single  arguments  used  for  array 
subscripts.  Treating  the  entire  array  item  the  same  as  a 
single  variable  was  found  sufficient. 

Since  arrays  and  function  subprogram  references  may 
have  the  same  syntactical  representation  in  an  arithmetic 
statement  but  would  be  n-bit  simulated  differently,  they 
must  be  differentiated  from  each  other.  To  distinguish  ar- 
rays from  functions,  each  program  unit  is  examined  to  de- 
termine the  variables  declared  as  arrays.  Arrays  may  be 
declared  in  either  DIMENSION,  INTEGER,  REAL,  or  COMMON  dec- 
laration statements.  So  for  each  variable  in  question,  all 
array  names  declared  within  that  program  unit  are  compared 
to  that  variable  name.  If  the  variable  in  question  is  not 
an  array,  then  the  variable  defaults  to  being  treated  as  a 
function  subprogram  reference,  and  n-bit  simulated  according- 
ly. 

Labels.  Another  potential  problem  area  for  FORTRAN  is 
where  a label  is  attached  to  an  arithmetic  statement  which 
is  the  last  statement  to  be  executed  in  the  DO  loop.  A la- 
bel in  FORTRAN  can  be  used  for  two  different  purposes. 
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These  usages  of  a label  will  be  described  before  defining 
the  problem  in  detail. 

A statement  label  uniquely  identifies  a statement  so  it 
can  be  referenced  by  amother  statement.  labels  occur  in  col- 
umns 1-5  of  a FORTRAN  statement.  The  actual  statement  ap- 


pears in  columns  7-72. 


Statement 

number 

Label 

Statement 

1 

C = X+2.5 

* 2 

GO  TO  38 

3 

14 

A = B 

4 

38 

C = B*D/5 

* Execution 
a branch 

of  statements  1 and  2 would  be 
statement  4 (through  the  use  of 

followed  by 
label  38) 

Figure  13  Normal  Use  Of  A Statement  Label 


Generally,  a label  marks  a statement  so  that  the  state- 
ment can  be  branced  ^ for  execution,  as  shown  in  Figure  13. 
However,  there  is  an  exception  in  the  case  of  DO  loops.  A 
DO  loop ‘statement,  such  as  the  one  shown  in  Figure  14,  is 
used  to  establish  a sequence  of  statements  to  be  executed  re- 
peatedly for  a specified  number  of  times.  In  the  case  of 
Figure  14,  a label  is  associated  with  the  last  state  it  of 
a DO  loop  sequence.  This  label  is  called  a DO  loop  termina- 
tion statement  label.  Upon  execution  of  its  associated 
statement,  program  execution  returns  to  the  first  statement 
of  the  DO  loop  sequence  to  execute  the  DO  loop  sequence  a- 
gain.  This  looping  would  continue  until  something  causes 
the  program  to  branch  out  of  the  loop.  The  main  difference 
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between  the  DO  loop  termination  label  and  a normal  label  is 
that  a termination  label  is  used  to  indicate  the  last  state- 
ment of  a sequence  to  be  executed  and  a normal  label  is  used 
to  indicate  the  first  of  a sequence. 


I 

t 

f 


I Figure  14  Example  Of  DO  Loop  Termination  Label 

[ 

I • 

j The  problem  with  respect  to  this  effort  arises  when  a 

I complex  arithmetic  statement  is  broken  up  into  a sequence  of 

I single  operation  arithmetic  statements;  especially  if  the 

I ' 

FORTRAN  program  is  auialyzed  and  translated  on  a statement  by 

o^atement  basis  (i.e.  each  statement  is  analyzed  without 
kiiowledge  of  the  context  of  its  use  with  regard  to  other 
; statements  of  the  program) , Each  label  type  must  be  handled 

i differently.  Figure  15  shows  an  example  of  wha^  would  hap- 

I pen  if  a termination  label  were  treated  the  same  as  s normal 

label  in  the  treuislation  of  a complex  arithmetic  statement 
to  its  single  operation  series  equivalent.  In  the  example, 
it  is  associated  with  am  arithmetic  statement.  Normal  label 
handling  would  distort  the  real  meaning  of  the  DO  loop.  In 
the  translated  code  of  Figure  15,  part  of  the  computation  is 

! 


Statement 

number 

Label 

Statement 

1 

DO  70  L = 1,20 

2 

A(L)  = A(L)+5 

* 3 

70 

C(L)  = D(L)+A(L) 

* Label  70  indicates  last  statement  executed  in  the 
DO  loop.  After  execution  of  statement  3,  L would 
be  incremented  and  the  loop  executed  once  again, 
until  L exceeds  the  value  20. 

Figure  14  Example  Of  DO  Loop  Termination  Label 
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now  positioned  outside  the  DO  loop  sequence.  Thus,  the  re- 
sults would  not  be  the  same. 


Before  translation 

After  translation 

DO  5 I = 1,7 

DO  5 I = 1,7 

5 SUM  = A(I)+SUM*2 

5 T1  = SUM*2 

SUM  = A(I)+T1  < — outside 

of  loop 
now 

Figure  15  Incorrect  Translation  Of  A DO  Loop 


To  assure  a consistent  meaning  of  labels  for  n-bit  sim- 
ulation modifications,  the  user  will  always  be  required  to 
assign  DO  loop  termination  labels  to  CONTINUE  statements.  A 
CONTINUE  statement  can  always  be  used  as  the  last  statement 
in  a DO  loop  sequence  without  changing  the  loop's  meaning. 

A CONTINUE  statement  will  never  be  modified  for  n-bit  simu- 
lation since  only  arithmetic  statements  are  changed.  Figure 
16  shows  how  a DO  loop  termination  label  could  be  moved  to  a 
CONTINUE  statement  without  changing  the  performance  of  the 
loop.  Figure  16  also  shows  the  single  operation  translation 
of  that  loop. 


Before  translation 

After  translation 

DO  5 I = 1,7 

DO  5 I = 1,7 

SUM  = A(I)+SUM*2 

T1  = SUM* 2 

5 CONTINUE 

SUM  = A(l)+T1 

5 CONTINUE 

Figure  16  Use  Of  CONTINUE  Statement  In  A DO  Loop 
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IF  Statements . Another  slight  complication  of  the  FOR- 
TRAN approach  is  translating  an  arithmetic  statement  contained 
as  part  of  an  IF  statement.  Figure  17  (a)  shows  such  an  IF 
statement.  A procedure  had  to  be  developed  so  that  arithmet- 
ic statements  contained  as  part  of  IF  statements  could  be 
translated  into  a series  of  single  operation  statements  with- 
out distorting  the  original  program's  meaning.  Figure  17 
(b)  restates  (a)  in  such  a manner  that  the  arithmetic  state- 
ment is  now  separated  from  -^e  IF  statement  and  could  be 
treuislated  into  a sequence  of  single  operation  statements, 
shown  in  (c),  without  affecting  the  original  meaning.  Two 
new  labels  and  two  GO  TO  statements  were  added  to  the  orig- 
inal sequence  in  (a)  to  preserve  the  correct  order  of  exe- 
cution. The  execution  paths  for  (a),  (b),  and  (c)  are  iden- 
tical as  shown  in  (d)  . 

Overflow  Tracing . Another  FORTRAN  consideration  was 
assisting  the  user  in  tracing  down  the  instances  where  an 
overflow  of  the  n-bit  wordlength  occurred  during  the  n-bit 
simulation  of  an  algorithm.  It  would  be  helpful  to  the  user 
to  indicate  the  name  of  the  routine  and  the  line  number  with- 
in the  routine  in  which  the  overflow  condition  was  detected. 
This  would  be  in  terms  of  the  original  FORTRAN  coded  algo- 
rithm emd  not  the  n-bit  simulation  modified  code.  In  addi- 
tion, specifying  the  operation  executed  at  the  time  of  the 
overflow  might  also  be  helpful. 

To  assemble  the  necessary  data  for  these  overflow  mes- 
sages, the  beginning  of  each  new  routine  must  be  detected. 
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Figure  17  Arithmetic  IF  Statement  Handling 


w 


New  program  units  (routines)  are  initiated  by  one  of  the  fol- 
lowing FORTRAN  keywords;  PROGRAK,  SUBROUTINE,  FUNCTION,  REAL 
FUNCTION,  or  INTEGER  FUNCTION.  All  other  routine  initiation 
keywords,  such  as  LOGICAL  FUNCTION  or  DOUBLE  PRECISION  are 
not  allowed  in  n-bit  simulation  programs,  by  one  of  the  the- 
sis assumptions. 

Once  the  beginning  of  a new  routine  is  detected  the  rou- 
tine name  is  extracted,  the  line  counter  for  that  routine  is 
initialized  to  one  and  incremented  with  each  succeeding  state- 
ment. When  an  arithmetic  statement  is  translated  into  a se- 
ries of  subroutine  calls,  the  routine  name  and  line  number 
containing  the  arithmetic  statement  are  included  as  arguments. 
Therefore,  when  an  n-bit  wordlength  overflow  is  detected, 
this  information  along  with  the  operation  being  performed  at 
the  time  would  be  printed.  Then,  the  user  could  easily  trace 
down  the  statement  which  was  the  source  of  the  overflow  con- 
dition. 

The  analysis  of  the  FORTRAN  language  with  respect  to  n- 
bit  simulation  now  has  been  described.  The  following  section 
describes  the  system  developed  which  incorporated  the  results 
of  the  FORTRAN  analysis  aind  enables  a FORTRAN  algorithm  to 
simulate  n-bit  wordlength  effects. 

N-bit  Simulation  Tool  Design 

The  design  of  the  n-bit  simulation  tool  was  divided  in- 
to two  parts.  Part  one  preprocesses  the  FORTRAN  algorithms 
before  they  are  executed;  translating  each  FORTRAN  arithmetic 
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statement  into  a series  of  subroutine  calls,  each  of  which 
performs  one  of  the  statements 's  computations.  Part  two  of 
the  design  develops  the  subroutines  which  would  perform  the 
n-bit  wordlength  effects. 


Preprocessor 

The  purpose  of  the  preprocessor  is  to  perform  a lexical 
scan  amd  amalysis  of  each  statement  in  the  FORTRAN  algorithm, 
then  convert  the  arithmetic  statements  into  a series  of  n-bit 
simulation  function  subroutine  calls. 


The  preprocessor  essentially  performs  two  major  func- 


tions ; 


1)  It  performs  a syntactic  analysis  for  arithmetic 
statements  and  converts  them  into  their  Reverse 
Polish  Notation  equivalent  representation.  (Fur- 
ther detail  concerning  the  syntax  analysis  is  con- 
tained in  Appendix  P.). 

2)  The  preprocessor  disassembles  the  Reverse  Polish 
Notation  representation  into  function  subroutine 
calls,  tailored  according  to  the  variable  data 
types,  and  arithmetic  operations  involved.  Each 
function  reference  replaces  an  arithmetic  opera- 
tion from  the  original  statement.  If  temporary 
variables  are  required  to  store  intermediate  re- 
sults, a variable  matching  the  dominant  data  type 
involved  is  created.  In  addition,  for  each  func- 
tion reference,  the  preprocessor  supplies  the  rou- 
tine name  and  line  number  of  the  FORTRAN  statement 

from  which  the  statement  was  derived. 
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The  syntax  description  for  an  arithmetic  statement  was 
shown  by  the  BNF  of  Figure  9.  Each  statement  fitting  this 
syntax  description  would  be  transformed  into  one  or  more  func- 
tion subroutine  references.  Each  function  reference  would 
replace  a single  arithmetic  operation.  Any  statement  which 
did  not  fit  the  syntax  description  would  not  be  affected. 
Figure  18  contains  an  example  of  a preprocessed  program. 

In  Figure  18,  the  temporary  variables  were  created  to 
hold  the  intermediate  computed  values.  Floating  point  var- 
iables are  stored  in  variables  preceeded  by  RTTTT  and  fixed 
point  variables  by  ITTTT.  The  last  letter  is  incremented 
from  A to  B to  C and  so  on  for  each  succeeding  intermediate 
variable  within  each  equation.  These  names  will  be  restrict- 
ed from  the  user,  but  are  named  so  that  they  are  not  likely 
to  conflict  with  any  variable  name  a user  might  create. 

The  function  reference  names  indicate  the  type  of  op- 
eration that  is  being  performed  from  the  original  arithmet- 
ic statement.  For  example,  statement  8 of  the  preprocessed 
program  performs  the  multiplication  of  15  and  B from  the 
original  program  statement  8.  The  last  three  letters  of 
the  function  name  indicate  the  type  of  operation  (e.g.  ADD 
for  addition,  and  MPY  for  multiply).  The  first  two  letters 
indicate  the  data  type  of  the  operands  being  passed  as  ar- 
guments. II  indicates  both  operands  are  fixed  point  vari- 
able data  types.  RR  indicates  both  operands  are  floating 
point.  R1  and  R2  indicate  one  of  the  operands  being  passed 
is  fixed  point  and  the  other  floating  point.  The  number  1 
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Statement 

# 

Original  FORTRAN  Algorithm 

1 

PROGRAM  SA^PLE  (INPUT,  OUTPUT) 

2 

DIMENSION  A(10),  K(30),  KEy(7) 

3 

CALL  SETNBIT  (24, 1 , 1 , 1 , 17,6,0,KEy) 

4 

B = 25.3 

5 

READ  *,  LI 

6 

LI  « LI 

7 

CALL  SUB1  (A,25,E) 

8 

NEXT  = A(1)+15*B 

9 

K(L1)  = NEXT/50*(A(2)+B) 

10 

STOP  "END  OF  SAMPLE" 

11 

END 

Original 

Statement 

# 

Preprocessed  Algorithm  Version 

1 

PROGRAN  SAMPLE  (INPUT, OUTPUT) 

2 

DIMENSION  A(10),  K(30),  KLT(7) 

3 

CALL  SETNBIT  ( 24 , 1 , 1 , 1 , 1 7 , 6 ,0,KEy ) 

4 

B = RASGN(25.3,  7HSAMPLE  ,6,1, KEY) 

5 

READ  *,  LI 

6 • 

LI  » 1ASGN(L1,  7HSAMP'LE  ,6,1,KE'Y) 

7 

CALL  SUB1(A,25,B) 

8 

RTTTTAA  = R1MPY(15,B,  7HSAMPLE  ,8,0, KEY) 

NEXT  = RRADD(A(1),RTTTTAA,7HSAM'LE  ,8,1, KEY) 

9 

ITTTTAA  = IIDVD(NEXT,50,7HSAMPLE  ,9,0, KEY) 
RTTTTAA  = RRADD(A(2) ,B, 7HSAMILE  ,9,0, KEY) 

K(L1)  = R1MPY(ITTTTAA,RTTTTAA,7HSAMPLE  ,9,1, KEY] 

10 

STOP  "END  OF  SAMILE" 

11 

L - 

END 

Figure  18  Example  Of  A N-bit  Simulated  Algorithm 
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in  R1  indicates  the  fixed  point  variable  is  the  first  argu- 
ment of  the  function,  2 indicates  it  is  the  second.  Table 
III  contains  a list  of  all  function  subroutines.  It  also 
indicates  the  arithmetic  operation  performed  and  the  operaind 


data  types  involved. 


Name  of  Function 

Operation 

Opereind  1 
Data  Type 

Operand  2 
Data  Type 

1 1 ADD 

+ 

I 

I 

RIADD 

+ 

I 

R 

R2ADD 

+ 

R 

I 

RRADD 

+ 

R 

R 

IIFNS 

- 

I 

I 

R1KNS 

- 

I 

R 

R2MNS 

- 

R 

I 

RRMNS 

- 

R 

R 

II^TY 

* 

I 

I 

RIKPy 

* 

I 

R 

R2MPY 

* 

R 

I 

RRFPY 

* 

R 

R 

IIDVD 

/ 

/ 

I 

I 

R1DVD 

/ 

I 

R 

R2DVD 

/ 

R 

I 

RRDVD 

/ 

R 

R 

IIEXP 

I 

I 

R1EXP 

I 

R 

R2EXP 

«« 

R 

I 

RREXP 

** 

R 

R 

lASGN 

I 

- 

RASGN 

S 

R 

- 

I - 

Integer  (Fixed 

Point) 

R - 

Real  (Floating  Point) 

Table  III  N-bit  Simulation  Function  Subroutine  Names 


There  are  two  forms  of  n-bit  simulation  functions:  sin- 
gle operation  replacements  and  simple  assignment  statement 
replacements  (e.g.  lASGN  and  RASGN),  The  general  format  of 
these  two  function  types  are: 

FUNCTION  (OPERANDI  ,0PERANI)2, ROUTINE, STMT#, FINALFLG, KEY) 
FUNCTION  ( OPERAND , ROUTINE , STMT#, F INALFLG , KEY ) 

The  OPERANDS  are  variables  which  will  be  n-bit  simulated 
by  the  function.  The  ROUTINE  is  the  name  of  the  routine  con- 
taining the  statement  and  STMT#  is  the  line  number  of  the 
routine  in  which  the  statement  occurs. 

FINALFIG  is  used  to  distinguish  intermediate  value  as- 
signments from  final  assignments  of  an  arithmetic  assignment 
statement.  Zero  indicates  the  result  is  an  intermediate  as- 
signment, one  indicates  a final  assignment.  Intermediate 
assignment  results  will  maintain  the  arithmetic  precision 
specified  by  the  user;  either  single,  double,  or  triple  pre- 
cision. -Final  assignments  indicate  the  final  computation  of 
the  arithmetic  statement.  The  resulting  value  from  a final 
computation  always  has  single  precision  accuracy. 

The  last  argument  of  the  function  is  KEY.  KEY  is  a 
seven  item  array  containing  values  computed  by  the  SETNBIT 
subroutine  and  is  used  by  the  n-bit  simulation  subroutines 
to  simulate  n-bit  wordlength  effects. 

Each  arithmetic  statement  in  its  original  form  appears 
in  a COMMENT  statement  (a  non-executable  statement  type). 

This  COMMENT  is  printed  just  prior  to  the  series  of  function 
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calls  derived  from  it  (Figure  19).  This  procedure  helps  im- 
prove readability  of  the  preprocessed  prograim  code  for  the 


user. 


Original  Algorithm  Statements: 

A = B+FUNC(A,5*K,K1)  where  FUNG  is  a FUNCTION 
IF(A.EQ.C)  LAST  = FINAL/R*S*25 .97 

Preprocessed  Algorithm  Statements: 

A = B+FUNC(A,5*K,K1)  < 

RTTTTAA  = RASGN  (A,  7HPR0GRI<1 , 1 , 1 , KEY) 

ITTTTAA  = IIKPy(5,K,7HIROGRM1,1,1,KEY) 

ITTTTAB  = IASGN(K1 ,7HPR0GRK1 ,1 , 1 ,KEY) 

A = RRADD(B,FUNC(RTTTTAA, ITTTTAA, ITTTTAB) , 

* 7HPR0GRM1,1,1,KEY) 

IF(A.EQ.C)  LAST  = FINAL/R*S*25 .97  < 

CONTINUE 

IF(A,EQ.C)  GO  TO  7901 
GO  TO  7902 

RTTTTAA  = RRDVD(FINAL, R,7HIR0GRM1 , 2, 0, KEY ) 
RTTTTAB  * R.RMY(RTTTTAA,S,7HPR0GRM1 ,2,0,KEY) 
.LAST  = RRMPY(RTTTTAB,25.97,7HrROGRM1,2,1,KEY) 
CONTINUE 


Figure  19  Preprocessor  Output  For  Two  Complex 
Arithmetic  Statements 

A brief  description  of  the  preprocessor's  operation  has 
now  been  given.  A detailed  description  of  the  preprocessor 
is  given  in  Appendix  C.  Chapter  5 contains  a description  of 
the  preprocessor  design  and  summarizes  the  functions  of  the 
modules  within  the  preprocessor.  The  following  section  de- 
scribes the  n-bit  simulation  subroutines  which  are  called  by 

the  n-bit  simulated  algorithm  built  by  the  preprocessor. 
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Stmt.  # 

1 

2 51 


C 


1 


C 

51 

2 

7901 


7902 
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N-bit  Simulation  Subroutines 


The  n-bit  simulation  function  subroutines  are  designed 
to  perform  the  n-bit  wordlength  effects  upon  the  operands  in- 
volved in  a single  arithmetic  operation.  A separate  fiinction 
subroutine  was  developed  for  each  operation  and  data  type 
combination  as  shown  in  Table  III . Each  routine  simulates 
n-bit  wordlength  effects  according  to  the  options  of  the  us- 
er. The  user  states  the  options  desired  by  a call  to  a SET- 
NBIT  routine  to  be  subsequently  discussed  in  detail.  The 
various  user  options  are  stated  in  Figure  20.  Options  2,  3, 
and  4 apply  only  to  floating  point  computations.  The  fol- 
lowing subsections  will  show  how  each  option  is  carried  out. 

User  specifications: 

1)  Total  number  of  bits  per  word  | 

2)  Number  of  bits  in  exponent  auid  number  of  bits  for  mantis- 
sa, in  addition  to  the  position  of  the  implied  decimal 
point  with  respect  to  the  mantissa.  The  sum  of  the  man- 
tissa and  exponent  bits  must  be  one  less  than  the  total 
number  of  bits  specified  in  1).  The  extra  bit  is  for  the 
sign. 

3)  Rounding  or  truncation  effects  for  computations. 

4)  Single,  double,  or  triple  precision  arithmetic  accuracy. 

5)  Overflow  messages  print  or  suppression. 

Figure  20  User  Specifications  (Options) 

Option  One . The  first  option  indicates  the  total  num- 
ber of  bits  in  the  wordlength  to  be  simulated.  The  SENBIT 
routine  computes  the  maximum  value  that  a fixed  point  word 
of  n-bit  wordlength  could  contain.  The  absolute  value  of 
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each  fixed  point  result  is  compared  to  this  maximum.  If  the 
maximum  is  exceeded,  overflow  has  occurred.  The  value  that 
overflowed  is  then  replaced  by  the  maximum  positive  value 
possible  if  the  positive  maximum  was  exceeded,  or  the  maxi- 
mum negative  value  possible  if  the  negative  maximum  was  ex- 
ceeded. 

Option  Two . The  second  user  option  defines  the  data 
format  to  be  used  for  floating  point  values.  The  different 
specifications  include  the  number  of  mantissa  and  exponent 
bits  along  with  implied  decimal  point  with  respect  ^ the 
mantissa  (either  to  its  left  or  to  its  right).  They  are  used 
by  SETNBIT  to  compute  the  maximum  and  minimum  positive  float- 
ing point  values  representable  for  the  n-bit  word.  If  the  : 
absolute  value  of  a floating  point  result  exceeds  these  lim- 
its, an  overflow  or  underflow  condition  has  occurred.  Un- 
derflow means  the  result  represented  was  a smaller  value 
than  could  be  stored  in  an  n-bit  word  with  the  given  number 
of  exponent  bits.  Upon  detection  of  underflow,  the  value 
that  underflowed  is  replaced  with  zero.  Overflow  detection 
is  treated  similar  to  fixed  point  overflow,  except  the  max- 
imum floating  point  values  are  used. 

Also,  for  floating  point  values,  the  proper  number  of 
significant  mantissa  bits  must  be  maintained.  This  is  a- 
chieved  by  a right-shift  (sign  extended)  followed  by  a cir- 
cular left-shift  of  the  floating  point  variable.  For  exam- 
ple, if  18  bits  were  specified  for  the  maintissa,  the  right- 
most 30  significant  bits  would  be  truncated,  as  shown  in 
Figure  21 . 
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(In  Octal) 


Original  value 

Right-shifted  30 

Circular  left-shifted  30, 


exponent 


48  bits  of  sig- 
nificance 


0000000000^^576037. 

exponent  18  bit  mantissa 
1 7465760370000000000 
18  bit  mantissa 


Figure  21  N-bit  Simulation  Of  Mantissa 


Option  Three.  The  third  user  option  is  rounding  or 
truncation  effects.-  If  truncation  is  specified  no  special 
action  has  to  be  taken,  because  the  GDC  6600/CYBER  74  com- 
puters truncate  by  default.  To  round,  however,  a value  e- 
qual  to  one  half  of  the  least  significant  bit  (LSB)  is  added 
to  the  absolute  value  of  the  word.  This  is  accomplished  by 
right-shifting  the  number  of  mantissa  bits  specified;  so 
that  the  least  significant  bit  is  in  the  second  bit  position 
from  the  right  of  the  CDC  60-bit  word  (the  59th  bit).  Then 
the  value  one  is  added  (using  fixed  point  addition)  to  the 
quantity  and  the  resulting  quantity  is  right-shifted  one  bit. 
If  the  LSB-1  is  a one,  then  the  result  would  be  rounded  up; 
if  it  is  not,  it  would  have  no  effect.  Figure  22  shows  a 
case  of  a 15-bit  meintissa  in  which  the  result  is  rounded  up. 
After  the  one  bit  right-shift,  the  floating  point  value  is 
left-shifted  back  to  its  original  position. 

A special  condition  may  result  from  this  technique  of 

rounding.  This  condition  occurs  when  the  addition  of  one  to 
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Original  value; 

00111 11J)01 1 01 0^1 1111 101 01jM  01^1 1000110101100011010001000000000 
After  right-shift  of  32  bits: 

0000000000000000000000000000000001 11110011010111111010111011 
Plus 

000000000000000000000000000000000000000000000000000000000001 
Sum  equals: 

000000000000000000000000000000000 111110011010111111010111 100 
After  right-shift  of  1 bit: 

00000000000000000000000000000000001 1111001101011111101011110 

Circular  left-shift  of  33  bits  (with  rounded  result) ; 
0011111001101011111101011 11000000000000000000000000000000000 


Figure  22  Example  Of  Rounding 

the  right-shifted  floating  point  value  results  in  an  incre- 
ment of  the  exponent.  The  n-bit  mantissa  contains  its  maxi- 
mum representation  in  these  cases.  Figure  23  shows  a round- 
ed value  for  a floating  point  value  with  a 15  bit  mantissa. 
The  rounded  result,  shown  in  (b)  has  no  significauit  mantis- 
sa bits  and  would  be  treated  from  there  on  by  the  computer 
the  same  as  a floating  point  zero.  However,  in  a n-bit 


Original  value 


(a)  17257777750031152437 


15  bit  mantissa 

Rounded  value  (incorrect)  (b)  17260000000000000000 


Correct  rounded  value 


(c)  17264000000000000000 


Figure  23  Special  Case  Of  Rounding 


I 


T 


'1 


computer  the  round  up  of  the  value  in  (a)  would  result  in  an 
increment  of  the  exponent  and  the  mantissa  would  have  a one 
in  its  most  significemt  position,  as  shown  in  (c). 

Option  Four.  The  fourth  user  option  is  to  specify  sin- 
gle, double,  or  triple  precision  arithmetic  accuracy.  This 
option  is  used  to  simulate  the  effects  of  a whole  complex 
equation  being  computed  with  multi-word  precision,  then  stor- 
ing the  result  as  a single  precision  value.  The  preproces- 
sor passes  a flag  which  indicates  if  the  n-bit  subroutine 
computation  will  result  in  an  intermeidate  value  or  a final 
assignment.  The  SETNBIT  routine  computes  the  bit  shifting 
needed  to  be  done  for  the  intermediate  values. 

A floating  point  15-bit  word  with  nine  mantissa  bits 
and  the  triple  precision  option  would  be  handled  as  shown  in 
Figure  24.  In  a triple  word  accumulator  of  a 15-bit  comput- 
er, there  would  be  59  significant  mantissa  bits  (i.e.  9 bits 
plus  two  additional  15-bit  words  used  to  extend  the  mantissa 
significance) . The  final  value  assigned  would  be  handled 
just  as  any  15-bit  single  precision  floating  number,  similar 
to  Figure  21 . 


Original  value 

17465760371015777615 

Right-shift  9 bits 

00017465760371015777 

Left-shift  9 bits 

1 746^760^1^577^00 

39  bit  mantissa 

Figure  24  Multiple  Precision  Accuracy 
Option  Five . The  fifth  option  can  be  used  to  print  or 
suppress  overflow  messages  when  the  n-bit  wordlength  is 
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exceeded.  A test  is  made  in  the  overflow  handling  routine. 
If  the  message  flag  is  off  (set  to  zero),  no  message  is 
printed.  If  the  message  flag  is  on  (set  to  one),  a mes- 
sage is  printed  indicating  the  line  number  and  routine  in 
which  the  overflow  resulted  from  and  the  operation  being 
performed  at  the  time  the  condition  was  detected. 

N-bit  Simulation  Sequence 

The  general  sequence  of  events  involved  in  simulating 
n-bit  wordlength  effects  for  each  n-bit  simulation  routine 
is  listed  below; 

1)  Perform  the  arithmetic  operation  using  the  two  in- 
put operands  and  storing  the  result.  Input  oper- 

' ands  are  assumed  to  have  n-bit  wordlength  signifi- 

cance already. 

2)  If  the  rounding  option  (which  applies  only  to 
floating  point  computations)  is  used,  rounding 
effects  are  performed  on  the  result, 

3)  Check  the  result  for  overflow.  If  overflow  con- 
ditions exist,  branch  to  overflow  handler  to  sub- 
stitute a value  for  the  overflowed  result  and  print 
an  overflow  message  (if  the  message  print  option  is 
on) . 

4)  For  floating  point  quantities,  truncate  the  mantis- 
sa of  the  result  to  the  number  of  bits  specified 
for  the  n-bit  floating  point  mantissa.  If  this  is 
an  intermediate  result  (indicated  by  the  Final  Flag- 

input  argument),  the  result's  mantissa  is  truncated 
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to  the  precision  specified  (single,  double,  or  tri- 
ple) . 

5)  Return  to  the  calling  routine. 

SETNBIT . The  subroutine  SETNBIT  computes  boundary  val- 
ues based  on  the  user  options  selected  (Figure  20).  The 
SETNBIT  subroutine  is  required  to  be  the  first  executable 
statement  in  each  prograim  unit.  The  user  specifies  the  op- 
tions desired  as  arguments  to  the  SETNBIT  subroutine.  Dif- 
ferent routines  may  execute  with  different  options.  Each 
routine’s  options  are  translated  into  seven  key  values  which 
are  localized  to  each  routine  in  aji  array  called  KEY.  The 
array  KEY  is  used  by  n-bit  simulation  function  subroutines 
to  carry  out  the  n-bit  wordlength  effects  specified  by  the 
user.  A description  of  each  of  the  seven  values  of  the  KEY 
array  and  how  each  is  computed  is  as  follows. 

KEY(1)  represents  the  maximum  fixed  point  value  repre- 
sentable with  n-bits.  The  maximum  value  is  computed  by  set- 
ting the  sign  bit  (left-most  bit  in  a CDC  60-bit  word),  then 
circular  left-shifting  n,  the  number  of  bits  in  the  word- 
length,  and  subtracting  one,  as  shown  in  Figure  25  (where  n 
equals  15). 

(In  Octal) 

1.  Set  sign  bit:  400000000000000000000 

2.  Circular  left- shift  15  bits;  000000000000000040000 

3.  And  subtract  one:  000000000000000037777 

Figure  25  Computation  Of  Kauximum  Fixed  loint  Values 
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KEY (2)  represents  the  maiximum  floating  point  value  pos- 
sible for  an  n-bit  floating  point  word  with  the  specified 
floating  point  word  characteristics  (exponent/mantissa  bit 
lengths  and  the  implied  decimal  point  position).  This  maxi- 
mum value  is  computed  by  determining  the  maximum  range  for 
the  n-bit  word  exponent.  This  exponent  maximum  is  computed 
exactly  in  the  same  manner  as  KEY ( 1 ) . The  exponent  range 
would  be  equal  to  one-half  the  maximum  exponent  represent- 
able. For  a five-bit  exponent  and  a ten-bit  mantissa,  the 
maximum  representable  would  be  (in  binary)  11111  or 
The  minimum  exponent  would  be  00000  or  2~^  (assuming  a bi- 
as of  10000).  If  the  implied  decimal  point  is  to  the  left 
of  the  maintissa  (making  the  whole  mantissa  a fractional  val- 
ue), the  maximum  exponent  is  decreased  by  the  number  of  man- 
tissa bits  (ten  for  this  example)  in  the  n-bit  word.  If  the 
implied  decimal  point  is  to  the  right,  the  maximum  exponent 
is  not  affected. 

KEY(3)  represents  the  least  significant  floating  point 
value  representable  with  the  specified  mantissa,  exponent, 
and  implied  decimal  point  position.  Figure  26  shows  this 
value  for  a iC-bit  word  with  a five-bit  exponent  eind  ten-bit 
mantissa  with  the  implied  decimal  point  to  the  right.  To 
compute  this  value,  the  fixed  point  representation  of  the 
mantissa,  shown  in  Figure  26  is  multiplied  by  two  to  the 
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Figure  26  Example  Of  Least  Significant  Floating  Point  Value 
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exponent  mantissa 


maximum  negative  31  exponent  (-31  for  a five-bit  exponent). 

If  the  mantissa  is  a fractional  representation  for  this  ex- 
ample, the  number  of  mantissa  bits  (ten)  is  subtracted  from 
the  maximum  negative  exponent  (-31).  So  the  least  signifi- 
cant value  in  this  case  would  be  the  fixed  point  mantissa 

value  times  2^”^^  “ or  This  last  value  would  be 

-31 

equivalent  to  .5  times  2 ^ . Any  value  computed  in  the  n- 
bit  simulated  algorithm  with  less  significance  would  con- 
stitute an  underflow  for  the  n-bit  word. 

KEY(4)  and  KEY(5)  simply  contained  the  options  indicated 
for  rounding/truncation  and  message  printing/suppression, 
respectively. 

KEY(6)  indicates  the  number  of  bits  a mantissa  must  be 
right-shifted  then  left-shifted  to  retain  the  proper  amount 
of  mantissa  significeince,  which  would  be  48  minus  the  num- 
ber specified.  Forty-eight  is  the  number  of  bits  in  the 
mantissa  of  CDC's  60-bit  word. 

KEY (7)  indicates  the  number  of  bits  a mantissa  must  be 
shifted  to  retain  single,  double,  or  triple  precision  ac- 
curacy, This  is  computed  by  adding  the  number  of  mantissa 
bits  to  another  full-wordlength  (n-bits),  for  double  preci- 
sion, or  to  two  wordlengths  for  triple  precision.  Then  that 
result  is  subtracted  from  48.  For  a 16-bit  word  with  a nine- 
bit  mantissa,  the  resulting  quantity  for  single  precision 
would  be  39  or  (48  - 9).  For  double  precision,  it  would  be 
23  or  (48-(9+16)),  For  triple  precision,  it  would  be  7 or 
(48-(9+16+16)). 
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Besides  precomputing  the  values  for  the  KEY  array,  SET- 
NBIT  performs  validation  checks  on  the  options  to  verify  they 
are  consistent  and  can  he  simulated  using  the  CDC  60-bit 
wordlength.  Failure  of  a validation  check  results  in  pro- 
gram termination. 

The  Overall  N-bit  Simulation  Tool 

To  use  the  n-bit  simulation  tool,  a user  must  first 
comply  with  the  restrictions  which  are  described  in  Appen- 
dix B.  Two  major  restrictions  are  to  debug  the  FORTRAN  al- 
gorithm from  all  syntax  errors  prior  to  execution,  and  to 
assign  all  external  values  of  the  program  (such  as,  through 
a read)  to  themselves.  For  the  variable  L5,  the  statement 
L5  » L5  would  reduce  the  value  in  L5  to  its  n-bit  wordlength 
significance.  One  additional  restriction  is  that  the  first 
executable  statement  of  each  program  routine  must  be  a call 
to  the  SETNBIT  subroutine  containing  the  user's  options, 
amd  the  array  variable  KEY  must  be  dimensioned  for  seven 
items.  The  SETNBIT  subroutine  call  allows  the  user  to  spec- 
ify different  n-bit  wordlength  characteristics  for  different 
subprograms  within  the  same  set  of  programs. 

The  preprocessor  will  read  the  FORTRAN  algorithm  and 
transform  each  arithmetic  statement  into  a series  of  func- 
tion subroutine  calls.  The  preprocessed  program  would  then 
be  compiled  and  loaded  with  the  n-bit  simulation  function 
subroutines  and  executed.  During  program  execution,  the  n- 
bit  function  subroutines  would  perform  the  arithmetic 
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operations  and  assignments  while  simulating  the  effects  of 
n-hit  wordlength  and  other  user  options. 


Appendix  B contains  the  User's  Guide  to  the  n-bit  simu- 
lation tool  amd  Appendix  C contains  a Detailed  Description 
of  the  preprocessor  and  the  n-bit  simulation  subroutines. 

The  next  chapter  analyzes  the  design  of  the  preprocessor  in 
further  detail  and  describes  the  modules  contained  within  it. 
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IV,  structured  Design  Of  The  Ireprocessor 

This  chapter  describes  the  final  design  of  the  n-bit 
simulation  preprocessor  which  resulted  from  the  application 
of  structured  design  techniques.  It  also  discusses  how  ap- 
plicable structured  design  techniques  were  found  to  be  for 
a small  system  design  (i.e.  the  preprocessor)  of  approxi- 
mately 1000  source  statements, 

! I 

Structured  design  is  "the  art  of  designing  the  compo-  ' 

nents  of  a system  and  the  interrelationship  between  those 
components  in  the  best  possible  way"  (REF  9:  7).  Basically, 
it  is  a form  of  top-down  design.  One  structured  design  meth- 
od is  transform  einalysis.  Transform  analysis  is  a strategy 
that  derives  initial  structure  designs  which  are  usually 
quite  good  and  require  only  modest  restructuring  to  arrive 
at  the  final  design.  The  quality  of  a design  can  be  meas- 
ured in  terms  of  its  reliability  (the  number  of  "bugs"  en- 
countered in  the  progrsun) , maintainability  (the  amount  of 
effort  required  to  fix  "bugs"  in  the  program),  and  modifi- 
ability (the  cost  of  chamging  or  extending  the  program). 

Data  Flow  Description 

Using  the  transform  analysis  approach,  a data  flow  dia- 
gram of  the  preprocessor  was  developed,  as  shown  in  Figure 
27.  The  preprocessor  inputs  a FORTEAK  statement,  then  scans 
the  statement  for  a sequence  of  characters  which  could  rep- 
resent an  object  of  am  arithmetic  statement.  An  object  re- 
fers to  the  variable  to  which  the  result  of  arithmetic 
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equation  is  assigned;  where  an  arithmetic  statement  is  of 
the  general  form  "object  = arithmetic  expression".  Any 
statement  not  having  a legal  object  is  regarded  as  a non- 
expression and  is  output  in  its  original  form  to  the  n-bit 
simulated  program  file. 

Upon  detection  of  a legal  object,  the  scsm  continues 
to  the  right  side  of  the  in  the  arithmetic  statement  to 
evaluate  the  expression  and  convert  the  statement  into  its 
equivalent  Polish  Notation  representation. 

The  right  half  of  the  data  flow  diagram  describes  the 
process  in  which  the  Polish  Notation  string  is  scanned  for 
operators.  Each  operator  found  will  result  in  the  build- 
ing of  a function  call.  To  accomplish  this,  the  operands 
which  are  being  operated  upon  are  found  in  the  Polish 
string,  then  merged  with  the  function  routine  name  and 
other  data  required  to  build  the  function  call  (such  as 
routine  name,  line  number,  and  the  reference  to  the  KEY 
array) . 

The  afferent  and  efferent  data  elements  are  defined 
in  the  data  flow  design  also.  "Afferent  data  elements  are 
those  high-level  elements  of  the  data  which  are  furthest 
removed  from  the  physical  input,  yet  still  constitute  in- 
puts to  the  system.  Efferent  data  elements  are  those  fur- 
thest removed  from  the  physical  outputs  which  may  still  be 
regarded  as  outgoing.”  (REF  9:  262).  In  the  case  of  the 
preprocessor,  the  lolish  Notation  string  representation  of 
each  arithmetic  statement  constitutes  the  afferent  as  well 
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as  the  efferent  data  elements.  All  processes  to  the  left 
of  the  afferent  data  element  are  dedicated  towards  build- 
ing the  Polish  string  from  the  input  FORTRAN  statement. 

All  those  to  the  right  are  intent  upon  converting  the  Pol- 
ish string  into  a series  of  subroutine  calls. 

In  developing  a program  structure  from  the  data  flow 
diagram,  transform  analysis  states  that  major  afferent  and 
efferent  data  elements  should  be  handled  by  separate  mod- 
ules and  that  these  modules  should  be  immediately  subordi- 
nate to  the  main  module  of  the  preprocessor.  The  basic 
module  structure  which  should  be  developed  from  this  data 
flow  diagram  is  shown  in  Figure  28.  The  fully  factored 
structure  which  resulted  for  the  final  design  is  shown  in 
Figure  29.  The  module  definitions  for  the  resulting  pre- 
processor design  and  the  module  interfaces  are  contained 
in  Appendix  C.  The  following  section  explains  the  differ- 
ences between  the  fully  factored  structure  of  Figure  29  and 
the  structure  that  would  be  expected  to  follow  from  Figure 
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Departures 

Ideally,  with  the  transform  analysis  strategy,  every 
module  of  the  afferent  branch  should  have  only  afferent  or 
transform  subordinates,  where  afferent  characteristics  im- 
ply that  data  is  floating  inward  (actually  upward)  with 
respect  to  the  system.  Transform  modules  transform  the  in- 
put data  closer  to  its  most  abstract  input  representation 
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Figure  28  Initial  Program  Structure  For  The  Preprocessor 

Data  Flow  Diagram 
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for  the  system.  Efferent  (data  flowing  outward)  modules  are 
expected  to  have  only  efferent  and  tramsform  subordinates, 
where  in  this  case  the  transform  modules  tramsform  data 
closer  to  its  output  form.  Anything  contrary  to  a consist- 
ent data  flow  (either  inward  or  outward)  is  termed  a depar- 
ture. The  departures  of  the  resulting  structure  of  Figure 
29  are  discussed  in  the  following  paragraphs. 

Two  of  the  departures  are  for  two  minor  efferent  mod- 
ules which  appear  inside  the  afferent  branch  of  the  pre- 
processor structure.  These  are  the  PRINT-STATEMENT  amd 
PFELIMINARY-OUTFUT  modules.  The  PRINT-STATEMENT  module  is 
called  when  any  module  within  the  afferent  branch  detects 
a non-arithmetic  statement.  Only  arithmetic  statements 
are  converted  into  Polish  Notation  and  subsequently  trans- 
formed into  subroutine  calls.  All  other  statements  are 
passed  along  to  the  n-bit  simulated  prograim  file  unchamged. 

The  PRELIMINARY-OUTiUT  module  is  performed  once  an 
arithmetic  statement  has  been  converted  into  its  Polish  No- 
tation equivalent  representation.  This  routine  sets  up 
prints  to  the  output  file  for  a comment  containing  the  o- 
riginal  arithmetic  statement  and  the  preliminary  output 
that  is  required  for  labels  and  IF  statements.  Perhaps 
this  module  should  have  been  included  under  the  efferent 
branch  module  structure,  but  it  requires  the  original  in- 
put statement.  No  module  in  the  efferent  branch  requires 
the  original  statement,  while  almost  every  afferent  branch 
module  does.  In  order  to  prevent  the  original  input  of  the 
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system  from  being  passed  unnecessarily  into  the  efferent 
bremch  side,  the  PRELIMINARY-OUTPUT  module  was  included  as 
the  last  step  of  the  afferent  branch. 

The  INITIAL-READ  module  which  is  directly  subordinate 
to  the  executive  module  represents  amother  departure.  This 
module  is  used  to  read  the  first  program  card  image.  All 
succeeding  reads  are  performed  by  the  GET-STATEMENT  module, 
which  operates  on  a 'look-ahead'  basis,  in  order  to  detect 
statements  that  are  continued  on  other  cards.  The  INITIAL- 
READ  module  does  the  initial  read  so  the  GET-STATEMENT  mod- 
ule does  not  have  to  repeatedly  test  for  the  initial  read 
on  all  cards  in  the  program. 

The  other  departures  of  the  final  structure  in  Figure 
29  appear  in  the  efferent  branch.  The  module  called  FIND- 
ING-THE-OPERATION  in  Figure  28  was  such  a short  process  that 
it  was  merged  into  the  SET-UP-CALLS  module.  Also,  the  GET- 
PARTS  module  was  merged  with  its  subordinate  (BUILD-SUB- 
CALLS) because  they  were  so  strongly  data  coupled.  Data 
coupling  measures  the  amount  of  data  passed  as  arguments  be- 
tween modules  as  subroutine  pareuneters.  In  this  case,  the 
GET-PAPTS  module  was  very  simple  amd  the  data  coupling  be- 
tween it  and  BUILD-SUB-CALLS  was  high  enough  that  a separate 
module  was  not  justified. 

The  Effects  Of  Structured  Design 

One  of  the  goals  of  this  thesis  was  to  determine  how 
applicable  structured  design  concepts  were  to  small  programs. 
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The  importance  of  structured  design  in  large  programs  has 
received  considerable  publicity,  but  the  programs  of  less 
than  5,000  source  statements  are  more  common  place  in  the 
programing  world.  The  n-bit  simulation  preprocessor  is  an 
example  of  a small  program  design. 

Part  of  the  evaluation  of  the  resulting  structured  de- 
sign is  based  on  its  comparison  to  a design  developed  for 
the  preprocessor  before  the  structured  design  approach  was 
begun.  The  previous,  or  first,  prograun  design  did  not  taJce 
a formalized  approach  to  structuring  the  design.  It  was 
based  upon  a flowchart  of  the  preprocessor  concept.  Many 
of  the  modules  resulted  from  blocking  off  portions  of  the 
flowchart.  The  first  design's  progrsun  modules  were  fairly 
intricate  and  hard  to  understand.  The  structured  design 
approach  emphasizes  that  the  structure  of  the  prograun  should 
resemble  the  structure  of  the  problem,  which  the  first  de- 
sign certainly  did  not. 

The  first  design  consisted  of  16  modules  that  were  in- 
terdependent on  each  other.  Modules  were  not  designed  to 
be  independent.  Modules  split  up  functional  tasks.  When 
a module  called  another,  certain  assumptions  were  being 
made.  As  a result,  a prograun  modification  to  one  module 
may  have  indirectly  affected  another.  Prograun  modifica- 
tions could  not  be  performed  without  analyzing  many  modules 
to  evaluate  the  impact  of  the  change.  The  whole  design  was 
essentially  tailored  to  solving  this  problem,  without  much 
consideration  to  how  future  changes  could  be  incorporated. 
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It  is  probable  that  as  the  program  was  modified  it's  complex- 
ity would  go  up  amd  it's  maintainability  down.  Though  the 
first  design  was  never  carried  out  to  completion,  it  was 
clearly  inferior  to  the  final  structured  design. 

The  structured  design  consisted  of  30  modules.  It 
seemed  to  be  simpler,  more  straight-forward,  and  more  gen- 
eral approach  to  the  problem.  These  modules  were  general- 
ly much  more  functional  in  nature.  In  fact,  ten  of  the  30 
modules  were  so  generally  defined  that  each  was  used  by 
three  to  eight  calling  modules.  A new  maintenance  program- 
er  could  more  readily  understand  the  design  since  the  struc- 
tured design  was  more  functionally  broken  down. 

Another  aspect  of  the  structured  design  was  that  var- 
iables were  initialized  at  the  lowest  possible  level.  In 
the  first  design,  nearly  all  initializations  were  performed 
by  one  module.  Initialization  in  the  first  module  where 
the  variable  is  needed  is  more  desirable,  because  when  that 
module  is  modified  it  is  possible  an  initialization  may 
have  to  change.  Since  the  modules  are  so  far  removed,  the 
change  might  easily  be  overlooked. 

Data  coupling  of  modules  was  used  for  the  preprocessor 
interface  in  all  but  four  cases.  In  three  of  the  cases, 
COMMON  coupling  was  chosen  over  data  coupling  because  in 
each  case  only  two  modules  needed  to  share  the  data  and 
these  modules  were  separated  by  several  module  levels.  COM- 
MON coupling  of  the  two  modules  was  believed  better  them 
passing  the  data  along  as  arguments  through  many  modules 


that  had  no  need  for  the  information.  By  COMMON  coupling 
the  two  modules,  it  prevented  other  modules  from  having 
needless  access  to  the  data. 


The  other  exception  to  data  coupling  involved  the  pass- 
ing of  a control  flag  as  a subroutine  argument.  This  flag 
was  passed  down  to  the  BU!^LD-SUBROUTINE-CALLS  module.  It 
indicated  the  subroutine  call  being  built  was  the  final  as- 
signment for  the  arithmetic  statement,  eind  no  temporary  var- 
iable would  need  to  be  created  for  storing  the  result. 

No  other  control  flags  were  used  specifically.  Alter- 
nate returns  from  subroutines  were  used  in  place  of  control 
flags.  These  alternate  returns  were  used  for  end-of-file 
conditions  and  for  signaling  the  detection  of  an  illegal 
arithmetic  statement.  The  alternate  returns  for  an  ille- 
gal arithmetic  statement  were  used  to  avoid  further  syntax 
analysis  of  the  statement.  Non-arithmetic  statements  do 
not  need  to  be  converted  to  Polish  Notation  form.  Myers 
(REF  6;  51),  however,  does  not  regard  this  usage  of  alter- 
nate returns  as  control  flags  since  the  return  indicates 
to  the  superordinate  module  that  the  function  failed. 

Myers  regards  this  type  of  return  code  as  data  since  it 
does  not  tell  the  calling  module  what  to  do,  but  says  'I've 
failed'  and  the  calling  module  can  respond  however  it  wants. 

One  other  design  improvement  over  the  first  design,  re- 
sulting from  the  structured  design  approach,  was  better 
scope  of  effect/scope  of  control.  Scope  of  effect/scope  of 
• r essentially  rationalizes  that  all  modules  effected 
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by  a decision  should  be  subordinate  to  (beneath)  the  module 
making  the  decision,  A consequence  of  disregarding  the  scope 
of  effect/scope  of  control  rule  is  that  decisions  are  repeat- 
ed as  control  flags  are  passed  back  to  the  superordinate  to 
be  retested. 

The  preprocessor  design  kept  the  scope  of  effect  for  a 
decision  within  the  scope  of  control  of  the  module  making  the 
decision.  Under  the  first  design  of  the  preprocessor  pro- 
grams, the  detection  of  an  illegal  arithmetic  statement  re- 
sulted in  passing  a control  flag  back  to  the  superordinate 
module  which  output  the  statement  for  that  subordinate's  de- 
cision and  several  other  subordinate  modules  which  detected 
syntax  errors.  In  the  structured  design,  the  statement  was 
output  by  the  module  detecting  the  illegal  syntax. 

Testing  The  Quality  Of  The  Design 

An  important  test  for  rating  the  quality  of  a design  is 
the  ease  with  which  modifications  can  be  made  to  the  exist- 
ing program.  For  one  proposed  change  to  the  program,  the  de- 
sign was  found  quite  adaptable.  This  change  is  described  in 
the  following  subsection.  The  change  is  explained,  then  the 
modules  and  interfaces  affected  are  listed.  In  addition, 
the  impact  of  the  change  upon  the  quality  of  the  present  pro- 
gram is  evaluated. 

One  change  was  to  replace  all  subroutine  calls  built  by 
the  efferent  branch  with  in-line  code  that  would  do  the  n- 
bit  simulation  effects.  The  in-line  code  would  be  tailored 
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to  the  user  options  in  effect  for  that  particular  routine. 
This  would  mean  fewer  decisions  would  be  made  (such  as  test- 
ing if  the  rounding  effects  option  is  on) . Also  subroutine 
linkage  time  would  be  saved  since  the  n-bit  simulation  code 
was  no  longer  performed  by  a subroutine.  This  would  result 
in  improved  execution  time  for  the  n-bit  simulation  of  an 
algorithm.  This  proposed  change  is  described  completely  in 
Appendix  D,  as  a future  enhancement  to  the  system. 

Basically,  this  change  would  require  modifications  to 
three  of  the  current  modules:  BUILD-SUBROUTINE-CALLS,  OUT- 
PUT-S INGLE-ASS IGNKENT,  and  NON-EXPRESS ION-HANDLER . Also  two 
new  modules  would  need  to  be  developed.  The  modifications 
to  BUILD-SUBROUTINE-CALLS  and  OUTPUT-SINGLE- ASSIGNMENT  would 
simply  be  the  replacement  of  one  block  of  statements  within 
each  one.  This  block  of  statements  presently  merges  portions 
of  the  function  subroutine  being  built.  The  new  code  intro- 
duced within  these  routines  would  merge  these  same  parts  al- 
ready present  in  a different  way  and  call  a new  module  which 
would  insert  these  parts  into  standard  positions  to  build 
the  appropriate  series  of  in-line  statements  for  n-bit  simu- 
lation. 

Pour  additional  interfaces  would  be  introduced  into  the 
present  code;  three  of  them  to  join  the  new  modules  to  their 
superordinate  modules  and  the  other  would  interface  the  two 
new  modules  through  the  use  of  COMMON  coupling.  No  present 
interfaces  would  be  changed  and  the  strength  of  each  module 
J modified  would  not  be  decreased.  The  changed  module, 
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NON-EXPRESSION-HANDLER, presently  detects  key  statements  and 
calls  the  appropriate  handling  routine.  The  CALL  SETNBIT 
subroutine  reference  could  be  detected  using  the  seime  tech- 
nique, then  a new  module  to  handle  the  statement  could  be 
called.  This  new  module  would  analyze  the  argiaments  of  the 
SETNBIT  subroutine  and  compute  the  key  values  normally  com- 
puted by  the  SETNBIT  subroutine  of  the  present  system,  then 
pass  the  keys  as  data  through  an  interface  to  the  other  new 
module  created  for  this  enhancement. 

All  changes  for  this  program  enhancement  could  be  ac- 
complished quickly,  easily  with  firm  confidence  that  the 
system  would  be  just  as  functional  and  easy  to  maintain  as 
before . 

Conclusion 

In  conclusion,  structured  design  techniques  were  found 
to  be  valuable  in  the  creation  of  a better  design.  The  fi- 
nal design  had  many  significaint  advantages  over  the  previ- 
ous design  which  relied  on  no  actual  methodology.  The 
structured  design  resulted  in  a simpler  system  since  its 
structure  resembled  the  problem  structure  and  the  modules 
were  functional  in  nature.  Simplicity  is  am  important  de- 
sign objective.  A program  that  is  simple  and  easier  to 
understand  has  a more  positive  effect  on  the  future  main- 
tenance and  modification  costs  of  the  prograun. 

Even  though  the  structured  design  analysis  took  con- 
siderably more  time,  it  was  felt  the  time  would  be  regained 
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in  the  succeeding  stages  of  the  preprocessor  system  life- 
cycle,  The  integration  of  the  preprocessor  modules  and  the 
debugging  went  quickly  and  smoothly  once  the  design  was 
coded.  The  testing  and  the  results  obtained  are  discussed 
in  the  next  chapter. 

Based  on  the  preprocessor  program,  the  structured  de- 
sign techniques  were  found  to  be  very  applicable  to  this 
small  progreun  design.  Perhaps  the  cost  of  poor  designing 
is  not  as  great  for  a small  program  as  for  a large  program, 
but  small  programs  comprise  a large  percentage  of  the  pro- 
graming systems  around  the  world.  Therefore,  the  overall 
impact  of  structured  design  could  be  even  greater  than  ap- 
plying it  to  only  large  systems.  The  next  chapter  will  re- 
view the  result  from  tests  performed  upon  the  preprocessor 
developed  from  the  structured  design  and  the  other  parts 
of  the  n-bit  simulation  tool. 
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V,  Testing  And  Results 

This  chapter  reviews  the  tests  performed  on  the  two  por- 
tions of  the  n-bit  simulation  tool,  the  preprocessor  sind  the 
n-bit  simulation  subroutine.  It  then  examines  some  of  the 
results  obtained  from  two  examples  of  FORTRAN  algorithms  to 
, demonstrate  some  of  the  effects  obtainable  through  the  n-bit 

simulation  tool. 

Preprocessor 

The  objective  of  the  preprocessor  is  to  translate  a giv- 
en FORTRAN  algorithm  into  its  equivalent  n-bit  simulation 
form.  This  process  converts  all  syntactically  correct  ar- 
ithmetic statements  into  a sequence  of  one  or  more  function 
subroutine  calls. 

A test  program  including  various  forms  of  each  FORTRAN 
statement  type  was  developed  to  verify  that  non-arithmetic 
statements  would  pass  through  the  preprocessor  unaffected 
and  that  virtually  any  form  of  an  arithmetic  statement  would 
be  correctly  translated  into  the  intended  sequence  of  func- 
tion subroutine  calls.  This  ssimple  program  is  shown  in  Fig- 
ure 30.  The  n-bit  simulated  version  of  this  program  pro- 
duced by  the  preprocessor  is  shown  in  Figure  31 . 

All  statement  types  were  processed  correctly  by  the 
preprocessor.  The  sequence  of  function  calls  for  each  ar- 
ithmetic statement  reflected  the  proper  precedence  order  of 
evaluation  for  the  operators  and  opersinds  involved.  (A  com- 
ment statement  containing  the  original  statement  preceded 
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the  n-bit  simulation  function  calls  generated).  For  example, 
the  arithmetic  statement  in  line  420  of  Figure  30  translates 
into  the  statements  shown  in  lines  590  through  630  of  Figure 
31.  In  this  example,  statements  that  would  not  normally  ap- 
pear if  n-bit  simulation  were  not  performed  are  designated 
with  a pointer  in  Figure  30. 

Certain  extremes  were  introduced  into  the  original  pro- 
gram to  demonstrate  the  versatility  and  range  of  statements 
that  could  be  handled  by  the  preprocessor.  A few  of  these 
extreme  cases  shown  in  Figure  31  were  continuation  lines 
(line  830),  function  references  with  greater  than  six  argu- 
ments (line  2270),  long  array  argument  lists  (line  2560), 
long  integer  and  real  numbers  (line  2490),  and  multiple  IF 
statements  within  a given  routine  (lines  920-1490).  The 
latter  case  was  used  to  verify  that  unique  labels  were  be- 
ing generated  from  a series  of  modified  IF  statements. 

N-bit  Simulation  Subroutines 

The  preprocessor  produced  an  n-bit  simulated  algorithm 
which  referenced  up  to  22  different  n-bit  simulation  func- 
tion subroutines.  Testing  these  n-bit  simulation  subroutines 
consisted  of  two  phases.  The  first  phase  validated  the  key 
values  which  were  computed  by  the  SETNBIT  subroutine.  The 
second  phase  involved  testing  the  22  n-bit  simulation  func- 
•^ion  subroutines  to  verify  the  n-bit  wordlength  effects  and 
other  user  options  were  being  properly  performed. 


SETNBIT  Key  Values  Verification 

A critical  portion  of  testing  the  n-bit  simulation  tool 
involved  verifying  the  correctness  of  the  key  values  computed 
by  the  SETNBIT  subroutine.  This  is  important  because  the 
performance  of  all  22  n-bit  simulation  subroutines  are  based 
upon  the  computed  key  values.  The  seven  key  values  produced 
by  the  SETNBIT  subroutine,  based  on  the  user  options,  are 
listed  below: 

1)  Maximum  positive  fixed  point  value  possible 

2)  Maximum  positive  floating  point  value  possible 

3)  Least  significant  positive  floating  point  value  pos- 
sible 

4)  Rounding  or  truncation  effects  flag 

5)  Overflow  message  printing  or  suppression  flag 

6)  Number  of  least  significant  bits  (LSB)  to  be  trun- 
cated from  floating  point  values  for  final  assign- 
ments 

7)  Number  of  LSB  to  be  truncated  from  floating  point 
values  for  intermediate  assignments 

Figure  32  shows  key  values  computed  for  various  user  op- 
tions. Figure  32(a)  has  specifications  for  a 24-bit  float- 
ing point  word  with  17  mantissa  bits,  six  exponent  bits,  and 
an  implied  decimal  point  to  the  left  of  the  mantissa.  In  a 
24-bit  word,  the  maximum  floating  point  value  representable 
would  look  like: 
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(b) 

r«-BIT  2FEC.IFICRTlQr4':  PPE: 
rHjr?BEF-  GF  IITii-  WCrD  =£4 

TPUrCPTE  ' 0 CF  FQijNC  • 1>  =0 

PRECISIOri  = 1 

OVEPFLDUi  MES;SFi6ESa>  NC  MESSFiGES  OO;-  =1 
« MflMISiF  FITS  =17 
I « EXPDhEHT  BITS  =6 
DECiMFiL  FT  LEFT  nD>  DP  RIGHT  a>  = 1 


KEY  VALUES  IH  CCTAL 
I KEY  < I > = 000  0 0 0 0 0 0 0 0 037777777 
5 KEY  <e>  = £ 0 0 077777G  0 0 0 0 0 0 0 0 0 0 
^ KEY  <3>  = 17014  0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
KEY  <4  > = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
r KEY  = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 

.KEY<6>  = 0 0 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 37  I 

;KEY<7>=  00000000000000000037 

i <IN  DECIMAL.) 

^ KEY<1>=  MAXIMUM  INTEGER  = S3SSG07 

; KEY<2>=  MAXIMUM  FLGATING  FDINT  VALUE  PGSSIELE  = £. S147£3£9££7E-*-14 
; KEY<3>=  LEAST  SIGNIFICANT  FLOATING  FT  POSSIBLE  = . 00003051757S1ES 
; KEY<4>=  ROUNDING  >.  1 • ^ TRUNCATION  -;:0>  EFFECTS  = 0 
, KEY<5>=  OVERFLOW  MESSAGE  PPir4T  • i;>  .'  SUFPRSSN - O:)  = 1 
! KEY<G>=  « BITS  TRUNCATED  FROM  RIGHT  OF  FLTG  PT  FOR 
; FINAL  assignment  =31 

KEY<7>»  « BITS  TPUNCATED  FROM  RIGHT  OF  FLTG  PT  FOP 
INTERMEDIATE  ASSIGNMENT  = 31 


Figure  32  (Continued) 
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(c)  . 

M-BIT  SPECIFICFtTICN::  RFE: 

NUMBER  GF  EITS-  WCFIi  =£4 

truncate  -0*  OF  PDUhr  <1>  =1 

PRECISION  = £ 

OVERFLOW  MESSAGES' 1>  "NO  MESSAGES ' O:'  =1 
ts  MANTISSA  BITS  =1? 

» EXPONENT  I; ITS  =G 
DECIMAL  PT  LEFT'  O)  OP  RIGHT  a > = 0 


KEY  VALUES  IN  OCTAL 
KEY<1>=  000000  I 0 0 0 0 057777777 
key  (£^  — 1 1'S'C'i  1 ( I I I 1 I • I 7 1' G 0 

' KEY  <3>  * 1 GG  0*1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
‘ KEY  <4>  = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 
KEY  <5  > = 0 0 0 0 0 0 0 0 0 0 0 'j  0 0 0 0 0 0 0 1 
key  <6  ) = 00 1?  0 0 0 0 0 0 0 C 0 0 0 0 0 0 0 3 7 
KEY<7>=  00000000000000000007 

<IN  DECIMAL- 

KEY<1>=  maximum  INTEGER  = S3S8G07 

KEY<2>»  MAXIMUM  FLCATING  POINT  VALUE  POSSIBLE  = £.  1474:?5G479-?GE+'? 

KEY<3.)=  LEAST  3IGr< IF ICANT  FLOATING  PT  POSSIBLE  = £. 3£S30G45G539E-1 0 
KEY<4'>=  POUNDING  - i:- ' TRUNCATION  -.O;-  EFFECTS  = I 
KEY<5>=  OVERFLOW  MESSAGE  PRINT • !>  . SUPPPSSN ■ 0>  = I 
KEY<6)=  « BITS  TRUNCATED  FROM  FIGHT  OF  FLTG  PT  FOP 
FINAL  ASSIGNMENT  * 31 

KEY<7>*  « BITS  TRUNCATED  FROM  FIGHT  OF  FLTG  PT  FQR  ' 

INTERMEDIATE  ASSIGNMENT  = 7 

Figure  32  (Continued) 


The  maximum  exponent  possible  for  six  bits  is  l31 , as- 
suming an  exponent  bias  of  lOOOOOg.  Therefore,  the  value 
shown  would  be  equal  to  the  contents  of  the  mantissa  times 
2*^^  , The  mantissa  above  represents  the  fractional  value  of 
(2^*^  - , When  multiplied  out,  this  fraction  times  2*^^ 

equals  2147467264.,  which  is  consistent  with  the  value  com- 
puted by  the  SETNBIT  routine  for  KEY(2), 

The  least  significant  floating  point  value  possible 
would  look  like 

00000001 0000000000000000 

The  mantissa  in  this  case  is  equal  to  .5  and  the  expo- 
nent is  equal  to  -31  (with  respect  to  the  bias  IOOOOO2). 

So  the  least  significant  number  possible  is  .5  times  2~^^ 
or  2.328306436539  * 10~^*^,  This  is  consistent  with  the  val- 
ue computed  for  KEY(3)  by  SETNBIT. 

The  verification  of  the  fixed  point  maximum  can  be  ac- 
complished by  looking  at  the  octal  representation  of  KEY(1) 
in  Figure  32(a).  The  least  significant  23  bits  are  ones; 
this  trsmslates  to  8388607  for  its  equivalent  decimal  rep- 
resentation. In  a 24-bit  fixed  point  word,  the  largest  val- 
ue would  have  1's  in  the  least  significant  13  bit  positions. 

Figure  32(b)  has  the  same  specifications  as  (a)  with 
the  exception  that  the  implied  decimal  point  is  to  the  right 
of  the  mantissa.  The  maximum  floating  point  number  possible 
for  this  representation  is  the  maximum  mantissa,  which  has 
a fixed  point  value  of  2^^  - 1 or  131071,  times  2*^^, 


yielding  2,81472829227  times  10^^,  This  value  is  also  ex- 
actly the  same  as  the  computed  value  for  KEY(2)  of  Figure 
32(b).  The  least  significant  floating  point  value  also  com- 
putes to  the  exact  value  shown  in  KEY (3)  of  (b).  The  maxi- 
mum fixed  point  value  does  not  change  with  the  moving  of  the 
decimal  point  since  only  floating  point  values  are  affected. 

In  Figure  32(c),  double  precision  arithmetic  accuracy 
was  specified.  KEY(7)  contains  the  number  7,  which  is  the 
number  of  bit  positions  that  should  be  truncated  off  the 
CDC's  48-bit  mantissa.  This  means  41  significant  mantissa 
bits  are  retained.  The  single  precision  mantissa  of  17  bits 
(shown  in  KEY(6))  plus  extra  wordlength  of  24  bits  equals  41. 

Testing  Of  The  N-bit  Simulation  Subroutines 

A unique  function  subroutine  was  developed  for  each  ar- 
ithmetic operation  and  each  possible  combination  of  input 
variable  data  types.  Twenty  of  the  22  function  subroutines 
developed  perform  arithmetic  operations  and  two  perform  sim- 
ple assignments. 

The  20  function  subroutines  can  be  divided  into  two 
general  classes:  those  which  produce  floating  point  values 
and  those  which  produce  fixed  point  values.  The  subroutines 
of  a given  class  basically  differ  only  in  the  first  state- 
ment they  execute.  That  statement  performs  the  arithmetic 
operation  being  simulated.  After  which,  they  are  constructed 
the  same.  For  instance,  for  the  operations  A*E  and  C+D  the 
first  statement  in  each  handling  subroutine  would  perform 
the  required  operation  upon  the  operands  and  store  the  result 
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in  the  variable  RESULT,  From  there  on,  exactly  the  same  n- 
bit  simulation  effects  are  performed  upon  RESULT  for  both 
subroutines.  Therefore,  verification  of  the  results  pro- 
duced by  each  combination  of  user  options  for  one  subrou- 
tine of  each  class  should  sufficiently  prove  the  other  sub- 
routines of  that  class  perform  correctly  also. 

Figure  33  shows  the  results  from  two  sets  of  fixed 
point  and  floating  point  additions.  For  each  case  within 
a set,  the  octal  representation  of  the  two  values  being 
added  are  shown  by  INFUT1  and  INFUT2.  The  NORMAL  OUTPUT 
represents  the  result  of  the  computation  using  the  full  60- 
bit  precision  of  the  GDC  computer.  The  N-BIT  OUTPUT  rep- 
resents the  value  resulting  from  the  n-bit  simulation  of 
the  two  inputs  being  added  together.  If  an  overflow  condi- 
tion occurs,  a message  is  printed  (e.g.  case  1C),  and  the 
maocimum  value  possible  replaces  the  old  value.  The  key  val- 
ues precede  each  case.  Case  one  is  executed  for  the  options: 
1 6-bits,  rounding,  single  precision  accuracy,  overflow  mes- 
sages, nine  mantissa  bits,  and  six  exponent  bits  with  the 
implied  decimal  point  to  the  left  of  the  mantissa.  Gases  1A 
through  1H  were  performed  with  fixed  point  addition.  Cases 
II  through  IP  were  performed  with  floating  point  addition. 
Case  1C  shows  an  example  where  underflow  for  the  16-bit  word 
was  detected.  Case  IN  shows  that  9 bits  of  significance  were 
maintained  in  the  mantissa  and  also  shows  a case  of  round  up. 

In  order  to  verify  correct  detection  of  overflow,  many 
borderline  cases  were  examined  (e.g.  cases  1A  through  ID). 


90 


r:EVp:o  = 

f.s'.r.  r;r‘f  ■*  • 

' •‘•rr-*'rrrr 

'•c  .Tfin  FI7'£'’p 

ieY(£)  = 

ir'-T’.rrri."' 

•; . • I';  ■'  i;  • . M m 

I'.F  ' FLCF7*;  .•  FCII'iT 

rEV'*'.'  = 

1 . a (.’‘.•r-  ■ * 

^ r, 

■ LE;-';pT  t IC-!  ■ ►ICFii'’  flT'7' 

KEV'Jp= 

0 0 *.  ' c v • 1 

• *’1 

F.7i..ii£  :r.p"  . 1 7,.  TCI..;  CFT 

kEY'^. = 

»'•  0 >■; ))  •;  * - 1.1 1;- 1.'  i;  ■. 

i V • 1 

C‘  CFFlOIi  MEirFOe  CFTICn 

kEY'|'..'  = 

L- 1.'  *.  1 ’ '.i  f'  • 1 

«;  • r 1 a *■ 

V E l-'  • 7 • - 

0 ^ 

i * .»  4 7 

iriT£C’cF  irTFi^t'F  F!.:.:TiCr-; 

r'-i  -’•‘t'-t- 

tiCf'f’F',  cv'Fi,' ■ I ; • f.' r t t'’iy 
M-  BIT  O'J ' F U T= C C ' 0 F C 


IMFUTl'=t'-i'''..'Ci''.'0'.'‘''o'’  '"OyrTrTF  = O'' i I'l.iOi'"!"'' 

r*0Fr'‘-t.  r.'j'r'U"  = Mi'  1 ' ' v'lii'"’ ( r7?r'i 
f<“Bi T ijuTF'ur=oi.ii.;c  I. 0 i ' n 0 1.  i.n.iuC'rrrr':- 


*c: 

= I Oi"'  r'r  r.i'  '•7777't-  If'Fi.  ~C 

'f  I..  : = '■  p,.  I '■  ' t I p'l  i.-'". : I",  p i:  i; 

i:-.'t:.FFl  fl.  .f-  Lll't  " 

ri-piT  0'.iTpi.i''=ij:! t'P.'C.rpc I'l  c I'ooi'CprTr?" 


r-  ! f ' F I.  ~c  ■=  C 1"  i1  '■  >;  0 1 j i" p''  0 0 ■'  0 1‘  00  c. 

p. ; Pi.'.  1 ■ i;  i; 

Lll-t  s:  1 FGP  FCriTICr  ♦ 


irUrt+i  + -i 

iH«''<.'Tl=777777T7r77?r770pX.r..:-.  :r)F'<JT£= 
ItCF'fp'FL  C'.''Fi.'’'  = ;’7,T777  7"777:'700p..'’.1 
M-BIT  D'.iri^'Jf=r77777777777777i:ii-i;''  1 


;'h"  * i'-'C'i.p 

Jp.Ci.p'psp  • '.PI  ' ' : ■' I..  ■"7: — 

'■'iF'-IOt.  L'.  I • ■ _ 'I  ■.|■:P 

0’  CFFL.L''  iri  : t O'l- F Llf'f.  ’• 
M-fiT  C'.'“Fi.'T=i.'i:.'  '•  ..i.v', O'.’i.  p'0i‘'7'777 


: FCF  Fr'I'IIICh 


Ifip.jT  1x777777' 
nCrfHl.  Cl.'F^I. 

C''.£FFLr.i 
'i-BIT  CU'P'J 


riFi.'  -r  - rln-p, 

'•p'Ft^'il.  I‘='.' 

H-BIT  C'-iTF'.- 


. fpt.  ..  1 • ► -p 

ri  r...  P I 

'.7.Fi>'h.  Cl. '^p. 


■fii;i':i:'  f I r-F  1.1 '£=777 777777 777~''~T 
1 FCP  Finn  TIC 


i.i  (..i.r*-  I ! I.- 
f 1 r.  1 j J.1  r 0 1.. 


1 

: --'T  FT' 


Figure  33  Examples  Of  N-bit  Simulation  Additions 

91 


SET  2 


( EY ■ 1 > » 

► EV  '.'.i'  = 
^EY<;j■  = 
‘■eY'4.= 
t ev  • 5;>  * 

► .EY'6j  = 
KEY • 7 J = 


ooiii,iic.’jo':-r  i.iiC  fii.rrrrr 
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Also  computations  which  result  in  values  well  within  the 
range  of  the  16-bit  word  are  shown  in  Set  1. 

Set  2 has  the  same  options  as  Set  1 , with  the  exception 
of  triple  precision  arithmetic  accuracy.  Cases  2F  amd  2H 
show  examples  where  the  resulting  value  was  rounded  up.  Al- 
so these  cases  show  that  41  significant  bits  of  the  mantissa 
are  retained  for  the  double  precision  arithmetic  accuracy. 


N-bit  Simulation  Tool  Effects 

After  verifying  the  two  major  parts  of  the  system  (i.e. 
preprocessor  and  n-bit  simulation  subroutines),  the  n-bit 
simulation  tool  was  tested  as  a whole.  The  sample  FORTFLAK 
algorithm  to  be  n-bit  simulated  is  shown  in  Figure  34.  It 
performs  very  simple  computations  designed  to  demonstrate 
losses  of  significance  for  shorter  wordlengths.  The  n-bit 
simulated  version  of  the  program  (built  by  the  preprocessor) 
is  shown  in  Figure  35. 

The  n-bit  simulated  prograun  was  then  executed  with  var- 
ious options.  The  results  from  these  various  executions  are 
shown  in  Figure  36  (b) , (c),  (d),  and  (e).  The  results  from 
the  execution  of  the  original  FORTRAN  algorithm  without  n- 
bit  simulation  effects  is  shown  in  Figure  36(a). 

Figure  36(b)  used  a 24  bit  wordlength  with  17  mantissa 
bits  (with  the  implied  decimal  point  to  its  left)  and  6 bits 
for  the  exponent.  It  was  performed  with  single  precision 

I 

arithmetic  accuracy  and  with  rounding  effects.  The  result- 
ing output  for  TCTAL  was  376784.  while  the  original  program 
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1 

oROGRflN  TES7  4 ( INPUT, OUTPUT, TAPE1= OUTPUT, TAPE2  = INPUT) 

000100 

nUENSION  K(l"i),KEY(7) 

000110 

Bl>55210. 

000160 

J 92*7.327 

000170 

'■  no  8 1=1,50 

000180 

; TOTftL=Pl/02*TOT AL 

' 000190 

CONTINUE 

000200 

: PRINT  T0T1L=", TOTAL 

000210 

1 HALF*TOrAL-TOTAL/2. 

000220 

: PRINT  HSL='  TOTAL  = “,HALF- 

C00230 

Xl=. 5791123 

000240 

X2=.  5790  995 

000250 

; oiff=xi-x2 

000260 

WT1IF==ni'^F»lJ00. 

000270 

I PRINT  OIFF=",DIFF,“  WTDIFF=“,WTOIFF 

000280 

1 XT*. 97265  ... 

000290 

i Xs.*.  9857 134  . : *■ 

* 

' 000300 

. XT.0IFF=X3-X4  , 

000310 

i RATI0=WT0IP'/X340IFF»1CC 

000320 

! PRINT  X-«L0IFP=:”,x34riFF,“  RATIO  =”, RATIO 

000330 

! FUNC»FUNC1(23..56.  34,5*i.‘ . ,1,0,RATIO>*1.25 

000340 

; PRINT  FONC=~,FUNC 

000350 

CALL  SU31 (PA T:0, X 34niF F-. 453, OlFF, TOTAL) 

000360 

PRINT  SU9  RATIO*", RATIO 

‘ 000370 

STOP  " TEST4  COMPLETE" 

000390 

FIJNCTIO'J  '•■Jt;Cl(a,5,T,J,K) 

0T1EN3I3N  KEY (7) 

(-2) +3** (-2) 

FIINCI»P.1/300  0 . 

IFIFUNC1.LT.  .:3025)  FI'NC1=. 0 J 00 25 

FETURN 
, ENO 

I S»«ROUTrNE  sun  (Rl,P2,R3,R4)  ' 

niMENSION  KEY(Y) 

P1»R1»R2yF3*R<. 

Rl»Rl**(-2> 

RETURN 

EN!) 


000400 

000410 

000420 

000450 

000460 

000470 

000480 

000490 

000500 

000510 

000540 

000550 

000560 

000570 


Figure  34  Example  Of  An  Algorithm 
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9f 0 .KEY) 
f 9,1, KEY) 


, 11,0,KFY) 
, 12,0, KEY) 
, 13,1, KEY) 


TESTt(IMOUT,0UTPUT,TftPEl=0UTPUT,T4PE2=IMPUT) 
IHENSION  K<10),KEY<7) 

^1*55210  . 

Pl*'?ftSSN{r>521j.,’'HTESTi»  , 3,0, KEY) 

•?*7.327 

92s^AS5NCr.32’,7MTE?T‘,  , 4,0, KEY) 

j no  a t*l,50 

! TnTftLsBi/PZ^rnTaL 

j PTTTT4a  = ^P0tfn(31,'i2,rHTFST4  , 6,0, KEY) 

I TorAL=»^4in{=!TTTTAA,TOTf  L,7HTEST4  , 6,1, KEY) 

i \ C04TIUUE 

; BP'NT  T0TAL=", total 

“ILFsTOTAL-TOTAL/’ 

I PTTTTaA=P2DYD(TOT5L,3, 7HT£ST^  , 

i HA.P*P»'fYS  (TDTaL,-7TT-^T  AA,7HTEST4 

; P’lNT  MALF=",HaLF 

] XI*. 5791123 

X1SRAS3N ( .57911 23, 7HTEST4 
! X’s. 5790995 

j X?s^asiN(  .5793995, 7HTEST4 

I nirF=Xl-X2 

1 l5tCF=R^MMF(Yl,X2,7HTEST4 

i ■ WTOIFFsOI''P»  1 3 0 9 . 

i ir9IFF=’=MOY (OIPC, 1333 .,7MTEST4  , 14,1,XEY) 

i P'^ItJT  OIFF  = ",-)irF,"  WTOIPF=“,MTDIFF 

i X’». 97265 

' X3s^AS3N (.97253, rHTE9T4  , 16, 3, KEY) 

'■  X4». 9557134 

, X4«RA95N(. 9557134, 7HTEST4  , 17,0, KEY) 

i XX4PIFe=X7-X4 

I X3«OIF5*2P4mS (X3, X4,7HTEST4  , 18,1, KEY) 

j PATIO=<<TOI?'^/X34nirc,i(;r 

' P’'»’TTAl*9P9\/n(WTnir'-,X3‘.0IFr,7HTE5T4  , 19,0, KEY) 

j paTIOs92*1DY(PTTT-AA,i;C.7HTE3T4  , 19,1, KEY) 

j PPINT  r349I=’=-  = “,X34''I5’E,“  ^ATIO  = ", RATIO 

j P')*IC=5'JSCl<2  3-5fi.T4,5*4r  . ,1,  A,  RATIO)*  1.25 

RTTTTAAaRi-oY(5 .^l. ,7HTpsT4  , 21,0, KEY) 
pTTTT43*R£33M(2T45s.3(,,7htE5T4  , 21,0, KEY) 

ITTTTAA  = IA<-.‘}M(1,'4T>^5T4  , 21,0, KEY) 

ITTTTA9=IAFGN'()  ,7HTE5T4  , 2l,C!,i<EY) 
».TTTTA2=RA3r.t4(RA'c:0,74Tc^-r4  , 21,0, KEY) 

97rTTA9E’..'''  N(-'3*!ri  ( = TTTTAn,:^TTTTAA,i:TTTAA,:TTTTA3,RTTTTAC> 
• ,7UT£5T-  , 21,0, KEY) 

r.;-|C=RRKPY(RTTTTA9,1.25,7HT£ST4  , 21,1, KEY) 

PRINT 

call  SU91  (RATH  ,X'hPIff-.4  = 3,9IFF, TOTAL) 

PPINT  SU”  RATIO*", ?£TIC 
?TOP  ~ TES74  C3*iPLETE" 

FN1 

EIMCTIO'J  FlINCK  A,9,I,J,C) 
nr^ENEION  <£Y(^) 

®laA**  (-2)  ♦=>**  ( -2) 

RTTTTAA*R2*<o(1  ,-2,’4Fl.l'  Cl 
PTTTTA9aR2>YO(9,.2,-'w5U'  Cl 
“IsRRAOOCRTTTTAA,  ■’TTTTa'',7HrU‘ICl 
P'HCl*31/3330. 

p''i:iaRR3V(Rl,  33j3.,7“ci)';Cl  , 4,i,kE7) 

T'<FU*I':i.LT.  .33C.2I)  Filf  c:».30JC25 
!®(5un:i.lt..:3:25) 
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FU’i  1 . uu  uu  , ^mi-iinu  i 

t;o-lTINUE 

RETURN 

END 

subroutine  SUBKRl, '’2,93,54) 
nUENSION  t<EV(7) 

®l=Rl>R?  + R3+‘’‘* 

OTTTT4ft=R=ftTn(Rl,='2,''HSl’5l  , 
R'-TTT4BsRRiDn(R'’’TTT4a,  R7,7HS'J.R1 
R1=RR409(RTTTTAS,R4,7HSI'B1  , 


»,U,KtT) 


3.0, KET) 

, 3,0, KEY) 

3.1,  KEY) 


Rl=Rl»*(-2) 

R1=R2EX»(R1,-2,7HSUR1  , 4,1, KEY) 

RETURN 

ENO  
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obtained  a result  of  376757.1993995,  a difference  of  approx- 
imately +24.  When  truncation  effects  were  applied  (shown  in 
(c))  to  the  n-bit  simulated  program,  the  results  were  slight- 
ly less  accurate,  TOTAL  equalled  376688,  a difference  of  -69. 
Employing  double  precision  arithmetic  accuracy  with  single 
precision  storage  (shown  in  (d))  did  not  impact  the  accura- 
cy. However,  double  precision  24-bit  storage  variables, 
shown  in  (?),  increased  the  accuracy  of  this  algorithm's  com- 
putation to  almost  the  exact  accuracy  accomplished  by  the  o- 
riginal  algorithm  in  (a). 

Double  precision  storage  can  be  simulated  by  keeping 
the  exponent  the  same  length  as  in  the  single  precision  24- 
bit  word,  but  increasing  the  mantissa  and  the  overall  n-bit 
length  by  24  (another  full  word  of  accuracy) . The  only  cau- 
tion a user  must  exercise  in  simulating  double  precision 
variables  in  this  manner,  is  that  the  fixed  point  maximum 
values  computed  for  detecting  fixed  point  overflow  will  not 
be  correct.  Therefore,  fixed  point  values  that  would  have 
resulted  in  overflow  in  an  n-bit  machine  may  not  be  detected 
by  the  n-bit  simulation  tool. 

The  results  shown  in  (g)  of  Figure  36  are  for  the  full 
60-bits  of  the  GDC  word.  The  n-bit  simulated  60-bit  results 
were  exactly  the  same  as  those  produced  from  the  original  al- 
gorithm and  the  CDC's  full  60-bit  accuracy. 

For  another  example  of  n-bit  wordlength  effects,  the 
small  FORTRAN  algorithm  shown  in  Figure  37  was  n-bit  simu- 
lated for  various  combinations  of  user  options.  The  original 
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Figure  37  Example  Of  A Second  Algorithm 


algorithm  resulted  in  a SUM  of  219.693088999.  Numerous  oth- 
er bit  lengths  were  tried  and  their  results  are  shown  in  Fig- 
ure 38.  Particularly  interesting  is  the  varying  effect 
rounding  has  upon  the  results  of  this  algorithm  at  differ- 
ent bit  lengths,  A 12-bit  word  with  rounding  resulted  in  a 
SUM  equal  to  11,  while  without  rounding  SUM  equalled  128  as 
shown  in  (a)  and  (d)  respectively  of  Figure  38.  Figure  38 
(b)  and  (e)  show  a case  where  rounding  Increased  the  accuracy 
of  the  SUM  from  204.5  to  213.5  for  a l6-bit  wordlength. 

Cases  (f),  (g),  (h) , and  (i)  show  the  results  from  various 
wordlength  and  exponent,  mantissa  combinations. 

Cases  (j)  and  (k)  show  user  options  which  were  incorrect 
and  flagged  as  erroneous  by  the  SETNBIT  subroutine.  Case 
(j)  had  too  many  mantissa  and  exponent  bits  specified  for  a 
16-bit  wordlength  and  case  (k)  specified  64-bits,  but  60-bits 
is  the  maximum  wordlength  this  n-bit  simulation  tool  cam  sim- 
ulate on  the  CDC  machine. 

Cases  (l)  and  (m)  contrast  the  results  obtainable  from 
this  algorithm  with  a single  word  storage  24-bit  wordlength 
machine  and  a 24-bit  machine  with  double  word  storage  capa- 
bilities, both  involving  single  precision  arithmetic  accura- 
cy. 

It  would  be  up  to  the  user  to  determine  the  best  choice 
of  bit  length  and  characteristics  to  provide  sufficient  accu- 
tacy  for  the  algorithms  being  considered. 
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Core  And  Execution  Costs  Of  The  N-bit  Simulation  Tool 


Although  not  regarded  originally  as  a limiting  factor 
in  the  design  of  the  n-bit  simulation  tool,  the  additional 
core  and  execution  time  requirements  for  n-bit  simulating  an 
algorithm  may  inhibit  some  users  from  readily  using  the  tool. 

The  resulting  n-bit  simulated  version  of  a FORTRAN  algorithm 
requires  substeintially  more  execution  time  and  core  space 
than  the  original  algorithm. 

For  the  algorithms  shown  in  Figures  30,  34,  and  37  the 
core  space  required  for  the  n-bit  simulated  version  was  four 
to  five  times  that  used  by  the  original  algorithms.  More- 
over, the  n-bit  simulation  subroutines  added  an  extra  2000 
words  to  the  core  size  requirements.  However,  for  the  types  || 

of  algorithms  for  which  this  tool  is  intended,  the  original 
algorithms  usually  do  not  require  more  than  5000  words.  An 
expansion  of  these  kinds  of  algorithm*!  by  a factor  of  four 
or  five  would  probably  not  impact  the  user  very  much.  Those 
algorithms  that  did  exceed  5000  probably  have  much  of  their 
core  space  reserved  for  dimensioned  arrays.  The  n-bit  simu- 
lation tool  would  increase  these  algorithm's  core  space  re- 
quirements to  a lesser  degree  because  the  core  expansion  fac- 
tor is  dependent  upon  the  number  of  arithmetic  computations 
within  an  algorithm. 

Execution  time  for  an  n-bit  simulated  algorithm  may  be 
a.  more  significant  problem.  The  programs  in  Figure  34  and 
35  were  modified  so  that  the  statements  within  the  program 
were  executed  500  times.  The  original  algorithm  required 
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,277  CFU  seconds  for  execution  while  the  n-bit  simulated  ver- 
sion required  2,4  CIU  seconds,  an  increase  of  almost  a fac- 
tor of  ten.  Other  tests  found  cases  in  which  the  CFU  execu- 
tion time  was  increased  by  more  than  30  times. 

This  tremendous  increase  in  execution  time  is  due  to  the 
fact  that  for  each  arithmetic  computation  within  a FORTRAN 
arithmetic  statement,  six  to  ten  FORTRAN  instructions  must 
be  executed  to  perform  the  necessary  n-bit  wordlength  effects. 
In  addition,  there  is  the  subroutine  linkage  time  required 
to  call  n-bit  simulation  subroutine.  This  introduction  of 
new  subroutines  into  a FORTRAN  algorithm  causes  the  FORTRAN 
compiler  to  save  all  register  values  somewhere  in  storage  be- 
fore ijumping  to  a subroutine,  then  restore  those  values  after 
returning.  An  increase  in  execution  time  is  not  surprising 
considering  these  factors. 

This  additional  execution  time  may  greatly  impact  a 
user's  readiness  to  take  advantage  of  the  n-bit  simulation 
tool  if  the  user's  original  algorithm  required  more  than 
120  seconds  for  execution.  The  n-bit  simulation  version 
of  that  algorithm  could  require  anywhere  from  1200  to  3600 
seconds.  Now  the  user  would  be  dealing  in  terms  of  hours 
instead  of  minutes  and  that  is  only  CFU  time,  not  total  run 
time,  which  would  be  even  more.  One  method  to  lessen  the 
execution  time  explosion  factor  is  discussed  in  the  follow- 
ing chapter. 
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Conclusions 


In  conclusion,  the  n-bit  simulation  tool  was  found  to 
be  an  accurate  and  fairly  versatile  tool.  Through  the  em- 
ployment of  various  options,  a user  could  simulate  the  ef- 
fects of  various  wordlengths  with  rounding  or  truncation  ef- 
fects, double  or  triple  precision  arithmetic  accuracy,  var- 
ious mantissa  and  exponent  bit  lengths,  and  simulate  double 
wordlength  storage  variables  (for  floating  point  values) . 

The  tool  does,  however,  cost  the  user  additional  core  space 
and  execution  time.  Nevertheless,  it  is  a very  useful  tool 
for  those  applications  where  these  limiting  factors  do  not 
come  into  play  to  a substantial  degree. 


VI.  Conclusions  And  Recommendations 


This  chapter  provides  the  reader  with  conclusions  aiid 
recommendations  for  the  n-bit  simulation  tool. 

Conclusions 

The  n-bit  simulation  tool  provides  the  user  with  a means 
to  evaluate  the  effects  of  varying  wordlength  upon  the  numer- 
ical accuracy  of  an  algorithm.  The  results  can  then  be  eval- 
uated against  the  performance  specification  to  determine  the 
wordlength  and  precision  requirements. 

The  n-bit  simulation  tool  requires  that  an  algorithm  be- 
ing n-bit  simulated  is  executed  on  the  CDC  6600  or  CYBER  74 
computer  systems  and  that  the  algorithm  is  expressed  in  a 
slightly  restricted  version  of  FORTRAN  IV,  as  described  in 
Appendix  B. 

Through  the  various  user  options  available,  the  n-bit 

simulation  tool  can  simulate  the  numerical  effects  for  com- 

\ 

putations  performed  on  a fairly  wide  range  of  computer 
types.  By  exercising  various  options,  the  user  can  speci- 
fy a computer  type  which  has  various  wordlengths,  varying 
exponent  and  mantissa  lengths  (where  the  mantissa  repre- 
sents either  a fractional  or  fixed  point  value),  rounding 
or  truncation  effects  for  computations,  single,  double,  or 
triple  wordlength  accumulators,  and  doubleword  storage  for 
wordlengths  of  up  to  24-bits. 

In  addition,  the  user  may  specify  different  n-bit  word- 
length  effects  for  different  portions  of  an  algorithm. 
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This  feature  would  be  especially  useful  in  simulating  the 
algorithms  of  an  integrated  system;  where  different  portions 
of  the  total  integrated  system  (such  as  the  Kalman  Filter, 
navigation,  flight  control,  suid  weapon  delivery  systems)  re- 
quire differing  amounts  of  precision. 

Pecommendations 

Two  areas  for  immediate/near  future  development  are 
recommended : 

1 , Include  am  option  for  the  n-bit  simulation  tool 
which  would  provide  the  user  the  exact  results  ob- 
tainable from  an  n-bit  wordlength  machine;  specif- 
ically giving  the  values  that  overflow  the  exact 
representation  which  would  result  from  am  overflow 
of  the  n-bit  wordlength  machine. 

2.  Modify  the  n-bit  simulation  preprocessor  to  sub- 
stitute in-line  code,  which  would  perform  the  n- 
bit  wordlength  effects,  instead  of  substitution  of 
subroutine  calls.  This  change  would  improve  the 
execution  time  of  the  n-bit  simulated  prograira. 

Case  One.  The  result  of  an  overflow  varies  for  differ- 
ent computer  types.  With  the  present  n-bit  simulation  tool, 
maximum  values  are  substituted  for  values  that  overflow.  In 
the  real  world,  it  may  be  that  the  most  significamt  bits  are 

simply  truncated  off  for  fixed  point  addition  overflows  while 

• 

multiplication  overflows  result  in  the  correct  mamtissa  and 
the  maximum  exponent  possible. 
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An  option  to  the  n-bit  simulator  which  would  produce 
the  exact  results  from  an  overflow,  would  enable  the  user  to 
analyze  the  exact  results  obtainable  from  an  algorithm  that 
is  processed  on  an  n-bit  wordlength  machine.  This  new  op- 
tion would  require  an  analysis  of  the  ways  in  which  differ- 
ent computer  types  overflow  and  incorporate  simulation  tech- 
niques into  the  n-bit  simulation  tool  accordingly. 

Case  Two , Presently,  n-bit  simulation  effects  are  per- 
formed by  substitution  of  subroutine  calls  for  each  arith- 
metic operation  of  an  algorithm.  Upon  execution  of  the  n- 

bit  simulation  version  of  an  algorithm,  the  subroutines  ref-  i 

j 

erenced  perform  the  n-bit  wordlength  effects.  However,  the  ; 

I 
i 

additional  time  required  to  simulate  these  effects  for  all 
operations  in  an  algorithm  is  a very  high  ratio  (15  to  40 
times)  to  the  algorithms  normal  execution  time. 

It  is  anticipated  that  the  introduction  of  in-line  code 
to  perform  the  n-bit  wordlength  effects  would  reduce  the  ad- 
ditional execution  time  by  at  least  one- third.  This  proposed 
enhancement  would  allow  the  tailoring  of  the  code,  performed 
for  the  n-bit  simulation  task,  to  the  user  options  in  effect 
at  the  time.  For  example,  the  present  system  tests  within 
each  subroutine  for  the  rounding  option.  The  subroutine  ex- 
ecutes differently  depending  upon  the  option.  By  tailoring 

the  in-line  code  to  the  user  option,  decision-making  time 

« 

would  be  saved.  In  addition,  subroutine  linkage  time  (sav- 
ing and  restoring  of  CCr'PASS  registers)  would  be  saved  with 
the  in-line  code  approach.  The  total  aunount  of  time  saved 
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could  be  substantial  depending  on  the  number  of  arithmetic 
operations  performed  within  an  algorithm. 

Although  this  enhancement  would  improve  execution  time, 
it  would  cost  the  user  twice  the  amount  of  core  space  re- 
quired by  the  present  n-bit  simulation  tool.  However,  this 
additional  core  space  may  be  more  acceptable  to  the  user 
than  the  additional  execution  time  of  the  present  system. 

The  changes  to  the  preprocessor  required  for  this  modifica- 
tion are  described  in  detail  in  Appendix  D. 
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Approaches  To  Finding  Extra  COMPASS  Registers 

The  solution  to  the  problem  of  finding  extra  registers 
would  be  to  find  some  standard  value  that  is  always  stored 
in  a known  register  or  to  devise  some  method  of  saving  and 
restoring  registers.  Several  possibilities  were  investigated: 

1)  The  AO  register  contents  may  be  predictable.  How- 
ever, no  consistencies  were  found  and  destroying 
its  contents  caused  irregular  results.  AO  is  used 
as  the  relative  starting  address  in  a block  copy 
operation. 

2)  The  CDC  manuals  stated  it  was  a standard  practice 
of  CDC  to  use  the  B1/B7  registers  to  store  the  val- 
ue 1 . If  this  were  the  case,  that  register  could 
be  used  as  a spare  register  at  anytime  as  long  as 
it  was  restored  to  the  value  1 when  it  was  no  long- 
er needed.  However,  several  cases  of  exception 
were  found  to  this  standard  practice-so  it  was  of 
no  use. 

3)  The  Exchange  Jump  instruction  was  known  to  save  the 
contents  of  all  registers,  but  its  use  is  restricted 
to  the  operating  system  for  switching  between  two 
central  programs,  leaving  the  first  program  in  a 
useable  state  for  later  re-entry. 

4)  If  CDC  had  a standard  register  saving/restoring  con- 
vention similar  to  lE.V's,  the  extra  register  problem 
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would  be  solved.  However,  CDC  does  not  have  one 
per  se;  instead  before  jumping  to  amy  external  sub- 
routine the  register  values  which  will  be  required 
later  (within  that  program  unit)  are  stored,  then 
restored  when  needed. 

5)  Another  theory  was  that  since  the  A and  X registers 
were  associated  one-for-one,  the  address  in  regis- 
ter (where  n is  1-5)  could  be  used  to  recover 
the  corresponding  X register.  This  would  enable 
the  current  value  of  the  X register  to  be  written 
over,  since  it  could  be  restored  later.  Once  one 
spare  60-bit  register  was  found,  all  other  registers 
could  be  saved  using  a series  of  register  moves  and 
write-outs  to  a known  save  area.  However,  the  val- 
ue of  X register  does  not  always  correspond  to  the 
associated  A register.  A technique  was  devised 
which  saved  the  contents  of  an  X whenever  it  was 
changed.  In  that  mainner  it  could  be  restored  later. 
However,  there  was  still  a need  to  have  an  associ- 
ated pair  of  A6/X6  or  A7/X6  registers  to  make  this 
possible.  All  this  would  require  a lot  of  analysis 
of  the  COMPASS  code  which  was  felt  to  be  beyond  the 
scope  of  this  thesis. 
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Appendix  B 


User*  s Kemual  For  The  N-bit  Simulation  Tool 

The  n-bit  simulation  tool  enables  the  user  a means  to 
evaluate  the  numerical  effects  that  varying  wordlengths  have 
upon  a FORTRAK-programmed  algorithm.  When  applied  to  a FOR- 
TRAN-programmed  algortihm,  this  tool  will  produce  the  near 
exact  arithmetic  accuracy  as  a computer  of  the  specified  n- 
bit  wordlength.  This  n-bit  simulation  tool  is  useable  on 
either  the  CDC  6600/CYEER  74  computer  systems  and  requires 
that  algorithms  to  be  n-bit  simulated  be  expressed  in  a some- 
what restricted  version  of  the  FORTRAN  IV  (EXTENDED)  pro- 
gramming language. 

This  User's  Manual  will  describe  the  seven  options  a- 
vailable  to  the  user  in  specifying  different  n-bit  wordlengths, 
the  using  programmer's  restrictions  will  be  described,  and 
the  required  control  cards  will  be  shown. 

User  Options 

Seven  options  were  considered  necessary  to  provide  the 
user  with  a versatile  n-bit  simulation  tool  which  could  be 
used  to  simulate  the  characteristics  of  various  computer 
types.  Different  options  could  be  applied  to  every  routine 
within  a program.  These  options  are  expressed  by  a call  to 
the  SETNBIT  subroutine  as  the  first  executable  statement  of 
each  routine  (i.e.  following  the  DIi;ENSION,  COMMON,  DATA 
statements,  and  so  on).  The  general  form  of  this  subroutine 
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call,  where  each  # represents  a user  option  and  KEY  is  a var- 
iable name  which  must  be  included  as  the  eighth  argument,  is 
shown  below: 

CALL  SETNBIT  (#1 , #2, #3, #4, #5, #6, #7, KEY) 

The  variable  naune  KEY  must  be  dimensioned  to  seven  (e.g. 
DIMENSION  KEY(7))  within  each  program  routine.  The  KEY  ar- 
ray is  used  to  store  values  which  are  computed  by  the  SETNBIT 
subroutine  based  upon  the  user  options  specified.  These  val- 
ues are  used  to  simulate  the  effects  of  n-bit  wordlength  for 
each  arithmetic  operation  of  an  arithmetic  statement. 

Option  #1  specifies  the  number  of  bits  per  wordlength 
(from  8 to  60  bits) . 

Option  #2  specifies  whether  arithmetic  operations  will 
be  performed  with  rounding  or  truncation  effects.  The  val- 
ue 1 would  indicate  rounding  effects  and  0 would  indicate 
truncation  effects. 

Option  #3  allows  the  user  to  specify  whether  arithmetic 
calculations  are  performed  in  single,  double,  or  triple  n- 
bit  wordlength  precision,  with  the  results  being  stored  as 
single  precision  quantities.  This  option  essentially  sim- 
ulates the  effects  of  a computer  having  a single,  double, 
or  triple  wordlength  accumulator  but  only  single  precision 
storage  variables.  Single  is  indicated  by  the  value  1,  dou- 
ble by  2,  and  triple  precision  by  3. 

Option  #£  gives  the  user  the  option  to  have  overflow 
messages  printed  when  they  are  detected  (i.e.  overflows  of 
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the  n-bit  wordlength) . The  value  1 would  indicate  overflow 
messages  will  be  printed  when  overflows  occur  and  0 would 
indicate  the  messages  would  not  be  printed.  An  overflow  mes- 
sage would  tell  the  user  which  routine  and  the  line  nvimber 
within  that  routine  where  the  n-bit  wordlength  overflow  was 
detected.  The  line  number  would  reference  the  line  number 
from  the  original  FORTRAN  program  (before  the  n-bit  simula- 
tion modifications  were  made).  In  addition,  the  arithmetic 
operation  being  performed  at  the  time  is  indicated. 

Option  #5  indicates  the  number  of  bits  which  are  to  be 
used  for  the  mantissa  of  floating  point  values. 

Option  #6  indicates  the  number  of  bits  which  are  to  be 
used  for  the  exponent  of  floating  point  values.  The  total 
number  of  bits  specified  for  the  mantissa  and  exponent  must 
be  one  less  than  the  total  number  of  bits  specified  in  op- 
tion #1  (this  other  bit  is  used  as  the  sign  bit). 

Option  #7  specifies  the  position  of  the  implied  decimal 
point  with  respect  to  the  mantissa,  either  to  its  left  (mak- 
ing the  mantissa  a fractional  representation)  or  to  its 
right  (making  the  mantissa  a fixed  point  representation). 

The  value  1 would  position  the  decimal  point  on  the  right 
aind  0 would  position  the  decimal  point  on  the  left  end  of 
the  mantissa. 

An  example  of  a typical  SETNBIT  subroutine  call  would 
be 

CALL  SETNBIT  ( 1 6 , 1 , 2,0, 9, 6 ,0, KEY) 
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The  options  shown  as  arguments  of  SETNBIT  specify  a 16- 
bit  wordlength  with  9 bits  in  the  mantissa  and  6 bits  in  the 

I 

exponent.  The  implied  decimal  point  would  be  positioned  to 
the  left  of  the  mantissa  (making  it  a fractional  representa-  1 

tion) . Rounding  effects  are  specified  for  arithmetic  calcu- 
lations, double  precision  arithmetic  accuracy  is  specified, 
and  overflow  messages  are  supposed  to  be  suppressed. 

Double  precision  storage  can  be  simulated  by  adding  the 
original  wordlength  to  the  mantissa  specification  (Option  #5) 
and  doubling  the  bit  length  specified  in  option  #1 . All  oth- 
er user  options  should  be  maintained  as  they  were  for  the 
single  precision  storage  word.  For  instance,  for  a 16-bit 
wordlength  with  9 bits  in  the  mantissa  and  6 bits  in  the  ex- 
ponent, the  double  precision  storage  could  be  simulated  by 
specifying  a 32-bit  wordlength,  a 25  bit  mantissa,  and  a 6 
bit  exponent. 

Double  precision  arithmetic  accuracy  should  be  distin- 
guished from  double  precision  storage,  with  double  preci- 
sion arithmetic  accuracy  (which  can  be  specified  under  op- 
tion #3)  a FORTRAN  expressed  arithmetic  statement  will  be 
computed  with  double  precision  accuracy,  but  after  the  final 
computation  of  the  statement  is  performed  the  result  stored 
will  contain  only  single  precision  accuracy.  In  contrast, 
double  precision  storage  would  have  double  precision  accura- 
cy, in  the  computation  as  well  as  double  precision  storage 
of  the  result.  Since  the  n-bit  simulation  tool  does  not 
presently  have  a separate  option  to  designate  double  precision 
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storage,  the  double  precision  can  be  accomplished  by  doubling 
the  original  wordlength  and  increasing  the  mantissa  by  one 
wordlength  (discussed  previously) , 

User  Progreunming  Restrictions 

The  user  programming  restrictions  which  must  be  adhered 
to  for  programs  to  be  n-bit  simulated  are  listed  in  this  ap- 
pendix. Restrictions  are  listed  in  the  order  in  which  they 
are  most  likely  to  impact  the  user.  Each  restriction  will 
be  briefly  justified.  Other  things  a user  should  be  aware 
of  are  listed  in  a section  following  the  restrictions. 

1 . Every  routine  must  have  the  variable  KEY  dimen- 
sioned to  seven  (DIMENSION  KEY  (7)).  Also,  the  first  ex- 
ecutable statement  within  each  routine  must  be  CALL  SETNBIT 
(seven  options  and  KEY),  where  the  user  specifies  the  seven 
options  to  be  applied  to  that  particular  routine.  These  ad- 
ditions are  required  so  each  routine  can  perform  the  n-bit 
simulation  effects. 

2.  Single  assignment  statements  are  required  for  all 
values  coming  from  a source  of  different  bit  length  speci- 
fications. Thus  all  values  read  in  from  cards  or  vdlues 
passed  as  parameters  from  a routine  of  different  n-bit  spec- 
ifications must  be  set  equal  to  themselves,  in  order  for 
them  to  take  on  their  n-bit  significance  (since  all  n-bit 
simulated  subroutines  assume  the  values  being  operated  upon 
have  their  n-bit  significance  already).  Therefore,  if  the. 
variable  A was  read  into  an  n-bit  routine  then  A should  be 


set  equal  to  itself  (A  » A)  before  it  is  used  anywhere  else 
in  the  program, 

3.  The  entire  program  must  be  free  of  syntax  errors, 
before  it  is  run  through  the  n-bit  simulation  preprocessor. 
Debugging  errors  would  be  easier  with  the  original  program 
and  the  preprocessor  has  not  been  verified  to  be  able  to 
handle  all  syntax  errors. 

4.  TAPE1  must  be  an  output  unit  for  each  n-bit  simu- 
lated program  (overflow  messages  are  printed  to  the  unit 
number) . 

5.  CONTINUE  statements  must  be  used  as  the  last  state- 
ment for  DO  Loops  due  to  the  way  labels  are  handled  by  the 
preprocessor  for  arithmetic  statements. 

6.  Variable  data  types  must  use  default  data  types 
(Reals  (A-H)  and  (0-Z),  Integers  (I-N)),  Declaration  of  da- 
ta types  using  INTEGER  and  REAL  statements  are  not  allowed, 
nor  are  IM LICIT  INTEGER  statements  allowed.  The  preproces- 
sor selects  handling  routines  based  on  default  data  types. 

7.  All  arrays  must  be  declared  by  DIKENSION  statements 
(not  by  COMMON,  INTEGER,  or  REAL  declarations).  This  is  so 
the  preprocessor  can  detect  functions  from  arrays  (it  finds 
all  arrays  in  the  DII'ENSICN  statement). 

8.  Logical,  double  precision,  and  alphanumeric  data 
handling  is  not  allowed.  The  n-bit  simulated  program  ex- 
pects all  numeric  values  in  assignment  amd  arithmetic  state- 
ments. 
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9.  These  are  certain  reserved  subroutine  names  (shown 
in  Figure  B-1)  and  reserved  variable  names  which  must  not 

be  used  by  user  (except  for  KEY  & SETNBIT).  These  names  are 
used  by  the  preprocessor. 

10.  No  statement  functions  are  allowed. 

11.  User  should  not  use  shift  operations  to  perform 
multiplications  or  divisions  (by  powers  of  two).  These  op- 
erations would  not  be  detected  as  arithmetic  operations  and 
would  not  be  n-bit  simulated  properly. 

12.  Data  statements  should  not  be  used  for  real  val- 
ues or  values  that  might  exceed  an  n-bit  simulated  word- 
length.  These  Data  statements  are  not  analyzed  by  the  pre- 
processor. 

13.  Reserved  labels  are  79110  - 79199. 

14.  Only  one  statement  per  line  maximum  (i.e.  no  A=B= 

0=0) . 

i 


Figure  B-1  Reserved  Variable  Names  For 
N-bit  Simulated  Irograms 

Other  User  Considerations ; 

1.  Intrinsic  functions  (such  as  SIN,CGS,ATAN) , other 

1 

FORTRAN  library  routines,  and  all  other  routines  not  read 
into  the  n-bit  simulation  preprocessor  are  not  n-bit  simu- 
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lated  completely  since  they  are  coded  in  COMPASS  and  the  pre- 
processor can  only  handle  FORTRAN  programs.  However,  the  ar- 
guments are  n-bit  simulated  and  the  resulting  value  from  the 
function  is  n-bit  simulated. 

2.  Array  arguments  are  not  n-bit  simulated  . 

3.  Only  assignment  and  arithmetic  statements  are  affect- 
ed by  the  n-bit  simulation  tool.  Expressions  occurring  else- 
where will  not  be  affected  (such  as  inside  IP  statements). 

Also  data  initialized  within  DATA  statements  will  not  be  af- 
fected . 

4.  For  exponentiation,  n-bit  simulation  affects  are 
only  applied  to  the  result  of  the  exponentiation  operation, 
so  the  results  do  not  reflect  on  an  n-bit  wordlength  machine 
exactly  in  some  cases. 

5.  Equivalence  statements  are  not  recommended,  but  can 
be  used  if  the  progrsunmer  fully  understands  the  impact  n-bit 
simulation  may  have  upon  the  equivalenced  values. 

6.  Any  user  should  be  aware  that  any  arithmetic  error 
which  would  result  in  abnormal  termination  in  a CDC  FORTRAN 
program  will  terminate  the  n-bit  simulated  program  also. 

Even  though  no  abnormal  termination  resulted  from  execution 
in  a normal  CDC  execution,  the  effects  of  n-bit  wordlength 
may  cause  a value  to  be  generated  which  could  terminate  it 
(for  instance,  .000001  might  be  truncated  to  0.  for  a 16-bit 
word  and  would  result  in  underflow  termination  when  used  as 
a divisor) . 

11 
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The  pointers  shown  on  the  left  of  Figure  B-2  indicate 
those  statements  which  would  not  have  appeared  in  the  pro- 


gram if  the  program  was  not  to  be  n-bit  simulated. 


Control  Cards 

The  following  card  sequence  could  be  used  to  n-bit 

simulate  a FOPTRAK  prograunmed  algorithm: 

FTK.L^DUMI'T. 

LGO, ,,NBIT. 

REWIND, KBIT. 

FTN  ,L=DUFiKY2  ,B=SUBS . 

REWIND, SUBS. 

EDITLIB. 

LIBRARY (NS UBS) . 

FTN , I=KBIT , B=NBITGC , L=DUKr-.Y3 . 

REWIND, NBITGC . 

NBITGO. 

7/8/9  Card 

***  Preprocessor  source  program  *** 

7/8/9  Card 

***  USER'S  SOURCE  PROGRAM  *** 

7/8/9  Card 

***  N-bit  Simulation  Subroutines  (Source)  *** 

7/8/9  Card 
LIBRARY (KSUBS, NEW) 

ADD (*, SUES) 

FIMSE  . 

EKDRUN . 

7/8/9  Card 

6/7/8/9  End-of-Job  Card 

To  n-bit  simulate  the  USER'S  SOURCE  ihOGRAF,  the  program- 
ming restrictions  previously  described  must  be  strictly  fol- 
lowed. Then,  the  USER'S  SOURCE  PROGP.AI'!  is  read  by  the  n-bit 
simulation  preprocessor  program,  which  transforms  arithmetic 
statements  of  the  original  program  into  subroutine  calls, 
which  are  designed  to  accomplish  the  n-bit  simulation  word- 
length  effects.  The  preprocessor  outputs  the  n-bit  simulated 
version  of  the  USE7-'S  SOUFCE  FRCGFAF.  Afterwhich  the  n-bit 
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simulation  subroutines  which  perform  the  n-bit  wordlength  ef- 
fects are  compiled  and  libraried.  Then,  the  n-bit  simulated 
user  program  is  compiled  and  executed.  If  the  program  uses 
additional  input/output  files,  these  files  would  be  put  in 
their  proper  position  along  with  the  NBITGO  card. 

If  the  object  forms  of  the  preprocessor  eind  n-bit  sim- 
ulation subroutines  are  used  instead  of  the  source  versions, 
the  control  card  sequence  would  be: 

COPYBR, ,SUBS,1 . 

INPUT, ,, KBIT. 

REWIND,  SUBS. 

EDITLIB. 

LIBRARY (N SUBS) . 

FTN , I=NBIT , B=NBITGO , L=DUMMY3 . 

RE WIND, NBITGO. 

NBITGC . 

7/8/9  Card 

***  N-bit  Simulation  Subroutines  (Object)  *** 
7/8/9  Card 

*•»*  Prepro'cessor  object  deck  *** 

7/8/9  Card 

***  USER’S  SOURCE  PROGRAM  **♦ 

7/8/9  Card 
LIBRARY(NSUES,NEW) 

ADD(*,SUBS) 

FINISH. 

ENDRUN . 

7/8/9  Card 

6/7/8/9  End  -of-Job  Card 

If  the  user  wishes  to  see  the  n-bit  simulated  version  of 
the  USER'S  SOURCE  PROGRAM,  the  L=DUMMY3  parameter  on  the  last 
FTK  card  should  be  removed.  Normally  this  print  would  pro- 
bably be  of  little  interest  to  the  user. 

The  user  should  be  aware  that  the  n-bit  simulation  of 
• a program  will  probably  increase  the  core  requirements  by 
four  to  five  times  and  the  execution  time  may  increase  up  to 
40  times  the  prograim's  normal  execution  time  (without  n-bit 
simulation  effects) . 
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With  the  present  n-bit  simulation  tool,  values  which 
overflow  are  replaced  by  the  maximum  significant  value  that 
the  n-bit  wordlength  can  represent.  For  example,  if  the 
resulting  value  of  a computation  exceeded  the  maximum  neg- 
ative number  representable  by  the  n-bit  word,  the  value  would 
be  replaced  by  that  maximum  negative  number.  In  addition,  an 
overflow  message  would  be  printed  if  the  user  option  was  on. 
Similarly,  an  overflow  of  the  positive  maximum  woiold  result 
in  the  replacement  of  the  value  that  overflowed  with  the 
positive  maximum  number.  Underflow  results  in  the  replacement 
of  zero  for  the  value  that  underflowed.  Same  results  of  the 
n-bit  simulation  tool  are  shown  in  Chapter  5 of  the  main  text. 
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Appendix  C 


Maintenance  Manual  For  The 
N-bit  STinulation  Tool 

The  n-bit  simulation  tool  consists  of  two  parts:  the 
preprocessor  program  suid  the  n-bit  simulation  subroutines. 

This  manual  includes  a description  of  the  overall  n-bit 
simulation  tool  system,  general  descriptions  of  the  pre- 
processor and  n-bit  simulation  subroutines,  structure  charts, 
interface  charts,  and  module  definitions  for  the  preprocessor. 
In  addition,  program  listings  of  the  preprocessor  and  the 
n-bit  simulation  subroutines  are  included. 

Overall  H-bit  Simulation  Tool 

The  preprocessor  reads  a FORTRAN  program  (to  be  n-bit 
simulated)  and  converts  each  arithmetic  statement  into  one  or 
more  subroutine  calls.  A subroutine  call  is  generated  for 
each  arithmetic  operation  of  an  arithmetic  statement.  The 
subroutine  calls, generated  by  the  preprocessor, reference  the 
n-bit  simulation  subroutines,  wfhen  the  n-bit  simulated  ver- 
sion of  the  original  algorithm  is  executed,  the  n-bit  simu- 
lation subroutines  perform  the  n-bit  wordlength  effects  upon 
each  arithmetic  statement  operation. 

The  n-bit  wordlength  effects  are  specified  by  the  user 
within  each  routine  of  the  original  FORTRAN  prograun.  The 
results  from  the  execution  of  the  n-bit  simulated  version 
are  nearly  the  exact  result  obtainable  from  a computer  of 
n-bit  wordlength. 
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Preprocessor 

The  preprocessor  basically  reads  and  syntactically  ana- 
lyzes each  statement  of  the  user's  program.  Each  arithmetic 
statement  is  transformed  into  a Polish  Notation  string,  then 
decomposed  one  operation  at  a time.  For  each  operation,  a 
subroutine  call  is  generated.  The  subroutine  referenced  will 
perform  the  arithmetic  computation  indicated  by  the  subroutine 
name.  Each  combination  of  variable  data  types  is  also  handled 
by  a separate  subroutine.  The  operauids  which  the  operation 
is  being  performed  upon  may  be  either  integer  or  real.  A dif- 
ferent subroutine  reference  is  generated  for  each  possible 
combination  of  operand  data  types  and  arithmetic  operation. 

The  module  structure  of  the  preprocessor  is  shown  in 
Figure  C-1.  The  module  definitions  follow,  along  with  a 
chart  describing  the  preprocessor's  module  interfaces.  The 
program  listing  of  the  preprocessor  follow  Figure  C-1 . 
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Module  Definitions  And  Module  Interfaces 
For  The  Preprocessor 

PPEPROCESSOR-KXEC  is  the  main  routine  of  the  preproc- 
essor program.  It  calls  the  afferent  and  efferent  branches 
(EVALUATE-STATEMENT  and  SET-UP-CALLS)  which  carry  out  the 
translation  of  a FORTRAN  algorithm  to  the  n-bit  simulation 
equivalent  algorithm  in  addition  to  the  INITIAL-READ  module. 

INITIAl-PEAD  does  the  initial  read  of  the  FORTRAN  al- 
gorithm. It  reads  the  contents  of  the  first  program  card 
into  a temporary  array  which  is  used  by  the  GET-STATEMENT 
module,  which  functions  on  a 'look-ahead*  basis. 

EVALUATE-STATEMENT  supervises  the  evaluation  of  (legal) 
arithmetic  statements  into  their  Polish  Notation  equivalent 
representation . 

GET-OBJECT  supervises  the  check  of  a statement  for  a 
legal  object,  where  an  arithmetic  statement  structure  con- 
sists of  object  = arithmetic  expression  . 

GET-STATEMENT  reads  input  program's  cards,  passing  a 
'statement*  a time  to  its  superordinates,  GET-OBJECT.  A 
statement  may  stretch  over  several  cards  through  the  use  of 
continuation  cards,  so  a look-ahead  read  is  performed. 

FIND-OBJECT  determines  if  a statement  has  a legal  ob- 
ject. 

• NON-EXPRESS  ION-HANDLER  calls  for  the  print  of  a non- 
expression  statement  and  determines  if  the  non-expression 
is  a key  statement  type,  for  dimensioning  arrays  or  the 
beginning  of  a new  routine. 
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NEW- ROUTINE-HANDLER  extracts  the  program  unit's  name 
and  initializes  the  line  counter  for  that  new  unit  to  one. 
(This  information  is  used  for  each  arithmetic  function  sub- 
routine call  built  by  the  efferent  branch) . Also  it  clears 
the  list  of  variable  names  declaured  as  arrays  for  the  pre- 
vious program  unit. 

DIMENS ICN-STATENENT-HANDLER  stores  all  variable  names 
which  have  been  declared  arrays. 

PRINT-STATEMENT  prints  the  statement  which  may  span 
several  cards  to  the  n-bit  simulated  program  output  file. 

EVALUATE- EXPRESS I ON  is  a recursive  module  which  deter- 
mines (through  a syntax  analysis)  if  the  expression  is  a 
legal  arithmetic  expression,  aind  converts  the  statement  in- 
to its  Polish  Notation  equivalent  representation.  It  is 
recursive  because  parentheses  may  be  used  to  enclose  another 
complete  expression  and  parentheses  may  occur  throughout  an 
expression . 

ANALYZE-CPEP.AND  determines  if  the  syntax  of  a string 
of  characteristics  constitiutes  a legal  variable  name  (a 
variable  may  also  be  an  array  or  function  reference). 

ANALYZE-KUMBER  determines  if  the  syntax  of  a string 
of  characters  constitutes  a legal  number. 

ANALYZE-P CSSIELE-APRAY  determines  if  a variable  name 
has  been  declared  an  array. 

CLASS IFY -CHARACTER  types  input  character  as  a letter, 
digit,  operator,  parenthesis,  decimal  point, or  comma. 


L 
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STORE-CHAPACTER-IN-I-OLISH-STPING  stores  a variable  name 
or  operator  in  Polish  string  according  to  Reverse  Polish  No- 
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tation  smd  the  operator  precedence. 

PUT-VARIABLE-NAME  stores  a string  of  chaoracters  in  the 
NAME  array  and  passes  back  to  superordinate  a pointer  to  the 
position  of  the  string  in  the  NAME  array. 

PPEIIMINAPY-OUTIUT  outputs  a COMMENT  statement  contain- 
ing the  input  statement  aind  outputs  preliminary  output  for 
labels  and  IF  statements. 

SET-UP-CALLS  scans  Polish  string  left  to  right  and  calls 
a handling  routine  for  each  operator  found. 

EUILD-SUE-CALLS  assembles  and  merges  parts  of  a subrou- 
tine call. 

PRINT-SUE-CALL  prints  80  character  card  imaiges  to  the 
n-bit  simulated  program  output  file. 

BUILD-FUNC-AFGUMENTS  merges  function  arguments  together, 
generating  rails  for  single  arguments  if  necessary. 

FINl.II.- FUNCTION  merges  the  function  name  with  its  argu- 
ment list  (oailt  by  BUILD-FUNC-AFGUMENTS)  and  generates  a 
single  assignment  call  for  the  function. 

SET-UP-SINGLE- ASSIGNMENT  sets  up  parts  for  a single  as- 
signment . 

HANDLE- UNAIY  places  a unary  sign  in  front  of  the  var- 
iable name  through  a call  to  the  PUT-VARIAEIE-NAME  module. 

OUTPUT-SINGLE- ASSIGNMENT  merges  parts  required  for  a 
single  assignment  function  call. 

I 
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GET-VARIABLE-KAME  gets  a string  of  characters  from  the 
name  array  according  to  the  pointer  given  and  passes  to  the 
caller  the  string,  its  data  type,  and  its  length. 

MERGE  merges  strings  of  characters  together. 
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SHAT  -HTc;  n^t'^ATOf  IN  THE  STACK 

C1312C 

35 

Tf  (•|ST<.'.".S*K".AX)  SO  TO  500 

or  3130 

Cu314C 

S-Ar:<  (NST<,1)  =0->S(2,I) 

CCTISO 

S-^CK  (NSTK,?)=CHAR 

003150 

no  TO  no 

00317(1 

r 

cor  US'"!  »4?(:n 

CC31S0 

<«G 

!'■  (NST<.GS.S*<‘'4X)  SO  tq  51.3 

003190 

»r*K»NsrK+t 

023200 

"*A?K  (SSTK, 1) rL^TNK 

OC3210 

C 

fJ’NK  SETS  LETT  I’PEN  SANK  PRECEPENOE  INTO 

STACK 

003213 

S"lCK(NSrK,?) =SM4< 

031220 

r,o  TO  13  3 

003230 

c 

GO'  s.ISNT  "AP-N 

303240 

5: 

r-  (NST<.E'’.l)  GO  TO  CCC 

207250 

TtiSTACKIMSTK.?)  . -O.LCfi«'Eti)  30  TO  54 

003250 

c 

'•■'ECT  ‘■•in  3P  A F^'ll''TIO^'  TOO 

203270 

7“  (MPQL  .r.T.  ool-'A  <)  SO  TO  30  3 

333200 

•nALsP“0!.»l 

003290 

PO'.tS’H  INPQL.  1)  =-S*t<'K(MST'‘,2) 

001320 

P'l'.rS'MuPOL,  2>  =3 

2C331C 

c - 

rr  TT  HJ5  4 !-J‘/''TION  f OTCATO"’  “ETIIRN,  INPICATOP 

A 10^0  TO  “O'.ISH 

OC  3320 

TT(STAC<(PST<.-»)  .SO  .Four)  SO  TO  Su 

2C3330 

NSTKaNSr<-l 

;C37*.(J 

c 

TP  IT  was  a ‘■'!N''TION  INPICATOR  RETUTN , IfiO ICt TO“  ADSEO  TT  ■’OUISM 

lC335u 

sn  n =5 

CC3T50 

c 

►in  i lENOtfE  LEST  PA 'EM  OR  ElINCTIOM  FROM  STACK 

001770 

MPTKs'!STK-1 

0 01,130 

sn  TO  133 

C':339C 

C 

PAS  PRDCFSST  |-, — CLPAP  OPE^ATO®  STACK 

33143C 

Tc; 

ITCiTACKAMS’’'*'.!)  .Fn,->oS)  no  -'0  90 

CC3410 

TP  (STAIK  ('1ST<  f2)  . '■''.LTASEN)  GO  TO  5C3 

237420 

*P  (NOPU.GP.POL  lAX)  GO  '0  =03 

j j3-3o 

j:3-4C 

OAi.  ISHCJPSL.l'^-S'ArKlMTTic,?) 

0C.145Q 

p-ii.I3>4(HPriL,?)  = -! 

.C3<*5t 

"S'-<  = ‘ISr<-l 

0C31.7C 

-•1  TO  'S 

.27430 

1.' 

•'PV.=:iP3.*l 

C22-9C 

"i. 'p  )L.  3T.  •’''L  i AT  ) r.r  TO  S j 3 

^ c 3 " . 

P ’l.  r 39  ('nOL  , ! t = ‘ A s 

:c3sii 

PP'.ISA  ('(PPL,  2)  = ■ 

307520 

no 

CAI7I*.IU£ 

•2  01570 

S-'JIPN  , T 

027550 

S'  . 

■"'•IT  EP<'A->  IN  colISH  MSTK=“,  MSTK  rPOL  = 

•*,NP0U 

uC  3-570 

ro  10  3 

30X510 

c*c^ 

2 0 V5.90 

s'ii.A3'.i-i'i--  <:'t,Lrn,TYPr:,N')oo<-) 

U 1 G j 0 

”'’!.TCIT  T‘l  T =■  AP' ( ' -■' ) 

CCTGIO 

r 

'IT-iA'-*!  'TO;:;  •,  0®  .-.V'TFP.S  t.;  tuf 

N.l 

'^r  «^54v 

'CXS12 

r 

'.•13  PASS'S  'AC'  TO  sirc-.xCACNtTF  {.  t'ATNTT? 

'0  ' 

T‘<  i 

: :xM3 

r • 

POSim*!  T 1 !"*ou  7hf  s— ’US  TTOf-n  I'l 

-M£ 

. . JS14 

si ; 13  i/Ni-  's/  1 ’ '••  ( 1 ; ; , ,Lsr  i 

ATi-  ICITI  r-!(it 

;:;S3D 

•1  • 

: :7j'A 

r 

'O-'OIITS  />  3S  -iri-AS  -.EClIPr')  To  STOPS  T^-  STRING 

; ; 3-, SI. 

M r • '•;=  C.Fr-i)  f 7 

■ i,  ‘ -.GC 

r 

STO’E  LP'IG'-i  r.F  fT  '.I'  G 'TST'JE  THE  li'’  O >0 

OF 

TA-:  srEMr. 

'C  Xsss 

ST'i  »; , ->1 

. :7  3'2 

• ••,.-/._  Tr-I*l  , A1  I -ysir 

•3  IMS 

r 

- . .1  : 1- 

AA  s 111,  , 09' 0 3 

TOIGAO 

1."  I=:.''T  I*. 

3C  i'’ja 
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I 


n noo  oonoo  no 

a»  rs*  04  Ml  04 


''•’’4  iOS/^’f’')'/ .'ON'/IP*  / 0042M 

P^?T»l  0C42S0 

L**T»2’  0 0<*26C 

o’TNT  (l,r)  (TIT(I)  ,I  = F‘'ST,LST)  0G<*2''0 

00'*2a0 

(.51*1)  r.o  TO  3ff  Ow<»29C 

CO'S  IN  COLl"")  7^  OP.  lU?  217, ETC  :C429«* 

^»i.-»»0UT  IM  "IO'ITZM'IATIOM  NAoK  034100 

rt’ir  <UST*5)  =C0‘IT  004110 

L‘:*»LST*72  004320 

po^TsPPST^Tr  304330 

<■•0  TO  3 CC4340 

“PTUO,(  C04350 

p’P  0 0 4360 

5'npO'ITINE  '■.'■T'?T«T(<:TNT)  ,=cTUPflS(Nl)  3 04370 

ALT  ’iTUP*'  II  15  ll5Pf)  whpu  A‘(  IQP  W#,  are*!  CETECT-3  004330 

GEIST-IT  pPiTS  IMf'jr  PoCG’AMb  C A’PS  , ’OSS  I A ' "irAfl'IENT'  SC*332 

AT  8 TT4-  TO  IT5  5IJ<- c;^0oniijt  T E ( GE  TOO  J)  . 1 OTITEHIMT  0C-*334 

TIT  5t>8N  SCO-aM.  CAPtiS  TH-’O'JGH  THE  USI  CONTINUim'l  CAPOS,  004135 
SO  A LO.TT-AMaio  opflr  is  PEQUIPEO.  00^'»36 

T-PCICIT  IMTEor=»  (A-7)  0C439C 

0T3£N5t3'l  S’’'''(l)  004430 


CH'IOM/CMTS/lICMT 
Cl<H0«l/r"'‘^P''/TP')?(7'») 
oo  t40‘l/MA)*-0/  U"-  1100,3)  ,L5T‘( 

'•AT A LSTCOl./*:  / ,P0«;/3PErS/,7EP0/l<»C/,  3I.)."</1P  /,ENnPtLF/3PE0P/, 

* 4AXCP(.5/i‘.  ./,I.  IC‘IT/1/ 

MCI'JG  '■13  r'lr’TAL  =FA3  , JUST  A noPO  OF  “LANKS 
LASTaLSTCCt 

’•IT'IALIZE  L5’’  I FTJ  NAPES 
usTNaa 

C-lpflX  POP  3 3 ■•'"'IIT  CApP 

*t(TP40(  1) . "C. iPC.OP.TEyP »i) .PO.lP*.OP.'’EH? (1) .EO.l’i)  GO  TO  70 

'^'8  2 Isl,L5T'''>,_ 

5”  !T<I)sTP  |P(:) 

CO  tTfl'IE 

■'■'0  I'lO'IT  CAPdCOL*^  t-72)  IMT)  TEMO 

'’"10(2,3)  (’•"  l'(  T)  ,T  = 1,U5TC0L) 

P8’*‘17  (7  pa*.) 

""(PO"  (’)  .•!•:.  ;)  CO  ■'C  5.) 

LT '1T  = LM"“T*: 

' -ca  'Pit  ■•IT  r.!'-'’"'-------- 

■PlTP^aC.)  .F P(l)  .PI. IP*. 0-..T"  •’(!)  .50. IP’)  CO  PJ  30 
TC  T ;nvj-;-i  in  l,f)P  rc  tj  1. 

TPC  TPpa(  a)  .PO  .07.  Tc  )'  = ( ^)  .EO.lDiO  .VO  TO  JO 

COT  :o*r TM'i\pron  lint 

'VxCDLs  "PO’  Tjrcaf:*  p jVpijQ'rs 50'  15  l-.*.'(2  CAP'’S) 

^LL  SriTP-  'P  I’P  OONCfo  TmA:J  fMIS  HILL  "0’  OE  AMALT'iO 
pp(UA5T.',p.'.'.t-;-ilC)  CO  TJ  fjO 

•P^aiOa’L'K 
'’FI  It,  :,i,l5':.)l 
•.T*L17P*l 
1*  "TiT(|.»5T)  = Tr  IMT) 

CO  TO  u 


C a IT  P05  •'A'F-'’  ft  LAST  CCLIH'l*!  OF  OUTCOftG  pT8TP;:i|* 


JC**420 
3C4-,3C 
3G>*-,4C 
CCi,4P0 
30*460 
JC4470 
aG4-*t30 
0 G <,490 
GG4uG0 
GC<*532 
.C451G 
334520 
OC4510 
CC45S0 

00. .553 
0045^0 
C C4.>70 
00-C3C 
0 C aOU 

'.I '*5 10 

1.  J * U t J 
j « 46  30 
004640 
Z L 4630 

0 L 4632 

1 CitobO 
C C 4670 
G3463G 
CC**60& 
3G*7:o 
0,4710 
CO-’H 
:;*77: 
jO'.’iO 
•C*>7'.0 
'C  ,74C 
Z L ‘‘"*C 
aC4770 
;g4mo 
:0*732 

• . .733 

. 't.TO'’ 
oc4):o 


5 rrjTiNti; 

"'■’’IJRN 

F-n 

S'lT^nUTI'IF  'r.j-T'’T,pfTuof..<;(,4^) 

’1’LICrT  (4--') 

C :*IT"’T  O'tFS  TfT'I’L  r'AO  .)=■  TM"  FCnTAH  /> 

C :T  ?E4T;  rn-  roraflTS  OF  TmF  CAnO  isn 

C A TF'1’0'?4*Y  t5  = 4»  WHTCW  xs  US£0  THE  'JET-STAY-i^-iJ' 

c 'lo.ti'LSt  ««’■';«  FU'ioricN*'  on  a look-ahfao  basis. 

cn  •'<0f</T£YP-''’/T'^MO(7»» 

OATA 

'■“40(2,?*)  ,1  = 1,72) 

2B  "■O’HAT  (*?oi) 

C 'JHFCK  "O’,  f ('’■V  FILE 

Y'*(eOF  (2)  .*)-. ')  .'>FTUP)-  Ml 

P-”')^M 

c>n 

S'nBP')Tr-|E  )''J  (S*NT  ,F)  jBET'JBNSCNl) 
r A'.T  Z'i  '(SFO  O^CF  AN  EOF  HA?  NrEM  OETECTEO 

C S -'PF’Y’FEt;  JUC  rMF''<  OF  a STtT'''EMT  'OP  A LES4.  OBJECT 

C MHF9=  an  A-’IT4”FTIC  ftaTENENT  STfuCTUF.E  CO^'SISTJ  of 

C OTJE**  = AFI'uMETir  FXP?E??ION 

T'*3LICr*  IN*'-.E={l-7) 

'**1ENFI3N  5T‘T  (1)  ,IF'?T"t“(3) 

rOHN0N/'>0l.’''‘u(-;/,j?2(_  ,,|c:tik 

n*,ri  TF^T-r/i  3- ,1->F  ,15  (/,3LM'</i=  /,'»f037wt/$/,IFTY'>f/i^1/, 

* rg-j/^rr-ic/.Loa^c  i/n  ( /,  -ajoE-j/i  p > / 

C rNCriALITr  - CT.'S  f,!,;  :ur  POLISH  FTFINGjSrN'.F 

r NE  HILL  4 riiw  one  fo<<  each  FTaxPHPxjT 

S ‘|->1L=NSr  < = 1 

rV.L  F,r75THT(-''“T)  ,fE7'J'- !1?  (6J5) 

•3r|sqP0?T>-:* 

C Fyrgn  ir-jo  c;Tt«’-  TVFF,yMTOH  OO'ITAIM  4V  4BIT4NETIC 

r TATFHF-IT  TO  )F  ExECUTEO  upon  A TBI'E  CPNOITION 

“O  15  1=1,1 

FALL  :L4  3IFY  < -*•-  , J)  ( = ic,  =cc) 

*“(FT <T(  )) . I-. 7 "CTMT (I) ) GO  TO  21 

1" 

C ■.■•)*  4'1  :=  TT-. -r  plate  iri  iN  ^OLISH  SEJIN*. 

0"*C*|T=1 

C f:>.'.C.E  OJEP  p-vr-ftTr  orpTiPN  OF  IF  STITFnHNT  TO  077  TO  O^J^CT 

?0  )»t*t 

’■F(ST1*(  J)  .-o.'O'.rTv,  ^4vCHT=o;Fr..iT-l 

T-r<3TNT(  J)  .Fo  .L  o'.PE'')  •■’4oC*;r=PAP0NT»l 

•Ff  ST 'T(  J)  ,--o  . r IP,  r.O  TO  -)„C 

■)  '■,o''*.o  2: 

C ',0  O.M<  NI'.i.  ‘I^TCHIM'  ►■'ir.HT  F.i.r,;  gp  jp  c;TmT  “’EOIOATE 

0 IS  FOUNT 

OA.L  3“=)L  (rF“Y''F, 

■>=  ( 

71  '■••■.L  F.no  = j(p-.,-, ,^PT|'5.-;S(3) 

->■3  rMO^, 

-OS  Fri  (•••,  , 13 -^TUM.-ITr  pTITC^lPIlT 

'■'tt  '".'.l  ■■’nn: ( = *••*) 

*0  TO  ? 

r Tg.i-.-;  - I.,.-!  J-J,  ; I.^ST  " T J T PM  (I  - ( >.  3?  | f.  P r,  T"|  nc 

r.  rTi-;T”-i  , 

'?')“  “)'.L  =~ir('-M-i 

I--.,,.,  ..L 

"•P  , 

t"  , yg,,  TT  , ,T  ( OijT  ) 

r ■>  >r  33I--  - I-  s-'-p..-  .■-■•.V  '•Fji.i  SF';FFV.  I’.jrj,  73- 

c - ir  N-i  r ' ■ • • rr-->~  Fti  f. 

'•  •-■-r. T-.j  - re.) 


003710 
003721 
0C3730 
u3.37i,0 
003750 
0C3760 
00376A 
C037f.5 
QQ376S 
’ 003767 
303770 
003790 
OC37<}0 
003900 
0C3902 
003910  - 
G0393J 
003840 
ii0  3950 
0C3d6a 
:u39&3 
t;C3967 
0n“')63 
01  'T 
u G A M 9 C 
00399C 
003900 
0F3910 
003030 

or  3033 

0C394Q 

003950 

;;*9&iJ 

or  3090 

003992 
!u  3990 
0 0 4 J 0 y 
CChOIO 
OC4020 
O'**  3 *3 
, J , 4 3 
rC4043 
::h05o 
5 y ^ y 63 
C04C7C 
004093 
00<*a90 
0 04u9“ 
C5-»ij99 
::4ioc 
OOhIIC 
00h12C 
y C ^1«0 
1)04143 
014150 
ocixirc 
004176 
: Cm  177 
COliO 
0Cm19C 
0 C m203 
ay421C 
::42?o 

OjA??* 
:.422'* 
yT42 JC 
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OUtJU  ou  uu  u uucuo  oo 


fit  can. 

r 3'IT  OUT  ALL  -'I^THCQ  cc*  TIMUarrONj  FOT;  THar  STATEI^MT 

Li'T*Lsr:oL 

^"aod,?)  (TFO  (I>,  iJl.lFTOCL) 

Te'fEO'f  (»)  .‘K.  j)  <;0  '0  50 

: SHiCK  F=  :T4'ic‘n  statement 

TC(TtMP(  U .os.TF'P  <i)  .E0.1’>.0R.TE“P  (1>  .EO.l?!)  GO  TO 

: ''aECK  ''O®  CT'-'Mijr  $-aTEH£MT 

T-<TEMe<S)  .co.tcp^.OP.TFMFJ?.)  .=0.3LH<)  GO  TO  1, 

(TO  TO  SI 

69  ran.  p-^NTtrEM^^T 

r. o  TO  SI 

: OT  a comment — p^iNT  n ano  pead  next  ca^o 

73  CALL  ®PMT(TE‘;-') 

P"'0(?,S)  (T5MP(I)  ,I=l,t.STC0L» 

Tr(fOP(»)  .ME,  :)  GO  TO  jJ 

U^IOUTsLNCH’’*! 

GO  TO  1 
"'ll 

SM'l^ODTr'IE  - MOOPJ  ($'■'(7  ,F)  .pET'IRMSTUD 

EM’iaOJ  OE’’"7“t'lES  IF  A GTATEMEMT  HAS  A L^GAL  OIJEGT, 

ai. r  p.etj^v  'll  IS  USED  when  illegal  apithetic  stmt 

implicit  I'|TF*.-p  (^.7) 
n’lE'isiO'j  S’'MT(i) 
ri'A  EO.TYPi/l/ 

CALL  CLASIPyI  ”''*t“tTvoF)  ,?"T'.i;?NS  (3C0 ,2U  3) 

CMtC<  03J-D’’  ir  AEIThEiTIC  STATEPEUT  TO  SEE  IF  IT'S  A LSGA. 
OOr’AMO/O-;  i’..^AY 

'•A.L  OP^MPf  P'  '*  .=  tTY^E)  »7ET1I^NF(2:C,2C3»2G0» 

:r<TTPE.Mr.'T>L’Y->^)  Gf  TO  2jii 

S-'T^e  FOUAL  Gfl  IM  o-jlI'H  ST-ETnO 

>nt  GO  PEC'.Esr  .-r-T  -jf  TrlTH'lETIC  STATEMENT 

raLL  SToolTST  1- (O) ) 

'’“TUIN 

GIT  NO‘I-A=:''m 'ITT',  STATEMEMT-  go  ChECX  key  STATEME-iT  types 
213  Cl'.L  ’IC'|Fy3<  ST '*1 

ES-')J'|  Ml 

F'tl 

'•’MPTITI'M?  Mr  isyrat^TMT) 

E’XITIT  :'I*--.E'(h-7) 

lO'l'yp  '■Al.'.  GnT  thc  cint  0'’  1 NON-E<"Fr.5<;iQ^  gTATEIEMT 
AVI  OFTT,  irMS'  IF  It  E f'n.|-cyP'^ESSnf  IS  A KEY  ST'^’’ 

■ ’•y>rjFp!’  o:  T-NriCMlIG  tY>AYS  OF  "EGIiiNIHr  r|F4  nJTIMES 


TT'ENsrON  .^T-TC.) 

-lA'A  KYw I's/s,  1 -G,  1 -1,1-  -r ITS, IT  1, -‘I, 

' 13,  '.T',*  "M-  'JflP'J.l'r.l  PT •'N,  IE-', 

. '.1,  1 1 M I-M.,  , 

12,1EP  '.F,  1-.  a,  IV'.  , • J-,  laU,  1®*',  t lEI.  1 o, 


I G > i3 1 0 f 3 > 
’ » G , 3 » 


I'-f  1®!  f ! I'-’' , :E'  , r-T  , 1 ir , , 1 (C,  1 MU,  ;CM,1»;,  IM^T  , 

« 1 T ^ 1 -"-a  , 1 - >'/ 

'".'a  ’L  l*'/'. s ' 

'•ur'fc'  T-.3  rr  ;,  n’S3'JT:'lE.''F'G'’A'',oc-6L  Fn'TTON* 

I , (>'  nr.TiCU 

■'  = 1 
LM-3 

-T  'j  Isl,') 

Pn_  '■) 

z (i  , *1  »< 

• j=  :•! 


co<*aiQ 

uaYS2o 
jC4)3C 
0C4940 
}C-»6S0 
OC407C 
004075 
S3  uC43a0 

0 G 4054 
OC4590 
004420 
3Q443Q 
:Cm440 
0G445O 
OC4460 
0Gm970 
0ti49dC 
4C499U 
OwGOOO 
005010 
0C5020 
005023 
CG5'J30 
C 05031 
305036 
005040 
CC5050 
J05C6C 
005070 
005076 
005077 
005050 
305u9C 
Ca513C 
305105 
0C5110 
0C5120 
uC5l23 
;c >130 
OCGIyC  ' 
0Cil50 
;:5i6o 
0:5170 
GG5171 
GG5172 
005173' 
0G5174 
005179 
0351SC 
uC51'lC 
j:52o: 
1.0  J2iC 

1 i3522C 
OCG’TO 
.C'240 

0 5250 
..OGDGC 
t05270 
•i:s»ar 

,i.529i 
10530!! 
OGSTIO 
0i,53»J 
i;;?33L 
0 i,»340 
uC535C 


oo">n>  o o onnooc^oon 


1 


1 


C IMPUT  S'iTCHr.jT  /GAINST  THf'ROS  ( ruaNK*"  ^^E  IS'OP.SQ) 

IT  OiOfl 

T*’(ST»iT(’)  .FT.ILNO  GO  TO  17 
IP/STXTIO)  .')r.<VW5C'5(Jf:)  ) 50  TO  20 

IS  OTiriN'j" 

c t3i  iNor:aT-s  rr  mas  a dimension  statemeht 

T«’H.SD.l)  50  TO  ?ao 
50  TO  40  a 
20  CriTTN'JE 

C IT  OT"*  NOT  NATSM  AMY  SEYMOROS 
50  TO  5? 3 

330  CALL  0I'1-NS'l(ST'<T,O) 

50  TO  533 

C COT  A NEW  PO'JTTflE-SO  SAVE  NAME  ANO  INITIALI7E  VIA  NW?JTIN 
400  CAUL  NWaiJ'TN(S’‘'T,P) 

530  PALL  ORNKST^T) 

return 

EN3 

S'HROIITINE  MW^'JTrM(?TMT,P) 

I-*»l.ICIT  lNrE5ER(a-7) 

NM‘*UTrN  EXTPaOTS  the  P»0G'’A“  UNITS  NAME  AMO  iMITiatlTES  THE 
LIME  COJMT^'R  POR  that  MFH  UNIT  TO  O'E.CTHIS  iMFOPMAriON  IF 
IS  USED  FOo  EAFH  a=>ITH'-ETI>D  FIIMOTirN  SURROUTIME  r.\..  3UTLT 
3Y  TUE  EPFEP-M"  APflMCH).  ALSO  IT  TLEAPS  THE  LIST  3?  YARieS'.E 
NAMES  DECLARED  AS  ARRAYS  FOR  The  PPEVIOUS  °RCGRaM  JNlT. 

TftvTS  ROUTINE  NAME  fluD  RE INITIALS £S  LINE  COUNT 

DIMENSION  STMTtt) 

PDMHON/aRRSYS/  AR^A  V((,0)  , ARRYMAX 
rOMMON/GOTON/  50'DNUM 
r'»MNON/RD"T:)A-1/  ROUTINE 
rOMMON/CNTF/  i.NFNT 

data  LPARFM/1R(/,-OS/3F''OS/,dlANXS/10H  /,D.N<71R  / 

RD JTINF=0l4NKS 
■'DNT  =3 

’’NTTIALITE  5TTD  n FOR  THE  r;EN  ROUTINE  JUST  DET'^CTEO 
TNTTTALTTE  Lit-:  "''UNT  Tr  t FJR  HEN  ROUTINE 

DL-AR  Array  ddun:  (s-nce  This  ue'i  -outine  ADESt.'T  hay* 
any  as  yet* 

51'gN')M=10 

LM>iT=l 

aR’YM.AX=3 

? Ps’*! 

L^’TT  PAR-N  INDIFATE*-  '■M'-  of  routine  name 
T'^«STHTI?1..-D  .i.dA-’EN)  FO  to  30 
•'stSTMTIP)  .f-OLNO  GO  TO  E 
’■P(S7"T(P)  .E^.EiDE)  50  TO  33 
D'DR*-'  CHfp'.F'EP  INTO  OUTINE  NAME  SHE  MCPD 
L--T=(?-p:Nr) ' 1 

pdjTI'4E=SM'<-T(S''»-(c:,  ,LrFT)  .OR.  (5MirTC1ASX(i,'.)  ,LE":i  .ANO.pQ  flNE) 
PD  tYiRCN’*! 

53  TO  5 
>3*  n-niRN 
-NO 

su’ROtJTIN-  a:  i-N'‘M(>JT''T,P) 

sricrr  » VA'-JAELE  NA 'ES  “hIDH  h\Y£  REFU  3ElLaRE0 
ARRAYS  To?  t;  T.icjiH  STATENEfiT 

■»n->r*P  H IS  rv-.— -r-  T-  'r  C-;rTH5  AT  L.''ST  L-TT'‘R  5'’  'IMENSICN  NORI 
’"LtDIT  i:-7\ 

-3l"3M/A-*'Avr/i7->„Y  (..(.I.APC.YUX 
■ ' ■t“'M~  13  .|  ? T a r ( . ) 

"■A  '.PAR. -./I  . 't  •.-'•'/'.P) /,;o'.'''TD/,.3.LNY/ip  /, 


035354 

0C536C 

0C5370 

005300 

005350 

JC5395 

005400 

005410 

305420 

005430 

005440 

C05450 

005460 

305470 

4C54S0 

305450 

005530 

305510 

CC5523 

C05530 

305531 

C05532 

305537  I 

CC5534  1 

JC5535 

C05-J36 

035537 

005540 

305549  ) 

005550  \ 

305560  j 

005570  ! 

405533 

3C5590 

0 3 5600 

3 35p13 

0C5630 

305623 

009624 

aOS625 

305626 

JC5633  i 

305643  1 

005650  j 

3 0 566C 

035665 

305670 

ccssao 

309690 
305695 
035730 
0C5710 
0C6773 
.39730 
305743 
309750 
C 05760 
C.C5T63 
; 157 75 

;C5776 
005777 
3 3 5776 
305730 
J337D0 
0C53QQ 
0 kSa  13 
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1 


k 4»1xRHMK«? 

HONTxO 
5 P*i>*l 

C '^■*EC<  -D’  =^?GINMIN<;  OF  OlHfNSION 

I'iSTlTta)  .“O.L'’a=‘EN)  HO  TO 
I'tSTHTCO)  .ci^.ijLuO  50  TO  5 
TP{STt1T<0)  .rO.^OS  J '•.0  TO  80 
LFFT»(q-'<c')r)  ‘S 

C •^TORC  eWflY  OHFH>IONFO  7AOIA3LE  MflMEdUE  CHAP4CTE1  4T  & TIME) 

NA'IxSHr- y<«;T'1T<o)  ,L?FT)  .09.  (fiaSK(r)CNT*8)  .AN0.N4M) 

►^MTaNCMrYl 
r.o  TO  5 

C '■•#$  "'cEN  DF-'=-CTrD 

C HE  MA'/E  POUNO  A ftPRAY  SC  SAVE  IT  THEN  PROCESS  ITS  A^SUIENTS 

SO  PC'IT*! 

AREYHAXxAPPY'nXtl 
A=PAY(  ARRY'IAX)  =mA»< 

C OPOCESSIMG  IIMTII.  HE  FIND  FIGHT  PAREN 
61  o=’*l 

tf(5THT( 3) . EO.Roaofn)  go  TO  67 
C IF  -PR09A3LE  FR’OP  STCO  PrOCESSSIMG 
! TF(3T1Tt»)  .E''.EOS>  GO  TO  30 

IF(SThT(P) .EO.LPApen)  fCNTx^CHT*! 
i GO  TO  fit 

< 67  rF(PC'<T.F0.1)  GO  TO  7L 

PCMTaPEVT-l 
GO  TO  SI 

C CO  IMA  GEPARATES  A?RY  OEFLARATIONS 

Tfc  as3«.l 

TF{STMT(o»  .FO.OOMMA)  CO  TO  :* 

TFfST>1T(=»  .FO.EOSJ  GO  TO  80 
GO  TO  7A 
as  "eturn 

F'n 

^•nPooTIME  'ET'/F 
ri’LICIT  I‘jr-GFP  fA-’l 

SFTJP  p-iu''?')  tTPp'G  left  TO  PIGHT  ANO  CALLS  T -U'lOLING 

P'TI’TM-  FT*  E'."^  OPE^f^OP  FOUST. 

PF  <E*(SIOS  ,T  AEf  f.y  (c  ) ,0TJECT  (F  ) 

OF'F«.|FIO  I PPS  ( ^)  , ^ULilFPl  T 1 

->»  *>5.1  ''o-j-fi) 

~P  'Aon/GOTT!/  n'OUIIM 
FO  tHO'J/CVT'/L  IC'I'' 

GO  IMOS/ROU*!'". '/  ppiriuE 
-OIMON/POL/pglISM (CD,?) y 

01*1  rr£Mr’ii/*MTT*TTAA/,»'*£yp-lC/7M?TTTTAA/,E0L/-‘rV/,  VEG/IH-/ 

OF  OPFFJ?.-*  / *'■  = ,?! 

'ATI  OOF/-?.'  ,-G’  ,-63/ 

Pin  IFTYC/lP^/.GTTPSUH/lu/.OOST/luH  GOMT,  mi>lJE  /, 

•*  COSTLi'i/l!-/ 

--T..O-J=f H 

o-tmp-JsR  'r  .o»p 

. F-IL  TS  A F'.AG  RMTGH  TOTi'ATiE  Tpr  LIST  OPPPATIFS  =f«fT?MEO 

I‘l  THE  tTr*H""T:F  STATC.IEM*  (OESCKIIfO  TY  T-ir  polish  NTTATIOS) 

c 

CML3  i 


FL-HLE'Isai 

G ’’IILP  SnsP'PP  "''!LC''FP.  FOI  'LL  SU  ' OIL'.S  TF-^jV  ER•)^|  TuF 

r POL  IF-*  G"*'  /G  ic  1 F-.-5LE  *•  ii*H'irr:c  sniFy.-sT 

F-tOOTF  (TLL  I . “LL  I-'P)  RC'IT  r''F,u(  C'lTjFSL 

I f-)3  .;t  / I , I • ,1  ■*. , T ; , ; I*. , : 1,  /i',  X'  y ) ) 

2 

>pFj\ro?F.i  F :ocFi  T(  <F  L'ss  thas  i is  polis«  gtjing 

G "(  P;l.  if  ■*'■,•)  . ..  T,'-.)  c ■ TT  1. 

■"  (Py).::-'  ( . '.  I . f.'.  .’rr'Fj)  Gl  TO  f ; 


ilOSSZO 

dcssao 

006840 

aossso 

Qosass 

(j&saso 

0S5370 

005380 

005390 

005894 

005900 

OC5910 

005930 

005932 

0 G394(i 
005950 
005960 
S 05970 
005990 
006000 
0 06(110 
006023 
GC603G 
0C604O 
006350 
006060 
OGbOTO 
006030 
006130 
006110 
036120 
006130 
00614C 
006150 
006160 
3 05 170 
JCS130 

oosiai 

006133 

acsiTo 

306230 

0:5210 
3:S2?0 
CC623Q 
3062*3 
OOS250 
3C62fi0 
OF j»70 
006230 
0CS2P0 
uCSTOE 
OCSJIC 
336320 
31.6325 

OOfiTJfi 

0:63?7 

FC6T30 

>lCS34d 

1 ' 5 345 
J.S346 
301350 
,35360 
3*6370 
..1310 
3 053 3C 
;f.SADC 
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Os3^1 

006»10 

ro  TO  5 

a06L2& 

c 

FINO  >tFXT  3 

306430 

I*; 

N’=?»l 

306440 

17 

T"  (“OLI'JH  (‘<•’,1).  LT.  C)  GO  TO  23 

00645C 

c 

POC  ■•=••  (SINGLE  ASSIGNMfH’  INDICATION) 

30 6460 

I"  (POLISH  (N'’,l>.  TO.  2REOS)  GO  TO  4C 

306470 

SC6400 

GO  TO  IT 

006490 

2C 

I'  (POLISH  ('l-'.J).  EQ.  FQL)  FNL  = l 

006'!  1)0 

INCOOr (pLLH-EN,1,plLWEP>  pOOTINE.LNCNT.FNL 

006910 

GO  TO  70 

306S20 

C 

.....  CHOHLO  ALWfTS  SET  FNL  = 1 ANO  STO®  IT  AL. 

306930 

c 

should  or  AM  "DUAL  SIG'  IF  IT  IS  THE  ONLY  OPr=,ATOT  IN  STRIMT 

006940 

4C 

IF  (OOLISH  (p,1).  ES.  EFL)  go  TO  i.2 

306990 

POINT*,  "EP’OT-MO  EDUAL" 

306960 

PETUPN 

006970 

C 

SET  JP  A single  ASSIGNMENT 

306990 

i,2 

CALL  GPT^N''(O0LIOH,",CNr) 

006990 

CALL  GETMA"- (COLTSM (ONO, 1) ,NAHE,TYPF,LEM) 

3C&S3C 

r 

HOH  STOPE  Tur  OOJECT 

036610 

CA.L  SETpNO^FILlFHjOxjQ.rPJ) 

3C6620 

CALL  Gcr'iSir  ( ^OL I EH  ( 0?  J . 1 1 , 0 IJECT  , 0"!  Y^F , OOLEN) 

036630 

c 

— NO  MFFO  TO  STOs;  NAME  S’TNCE  IT  IS  THE  LAST  0(1^ 

306640 

CALL  SMSLAS*.  (MANE,  LEM,  TYPE,  OAJcCT,  OcL  EM,  FLLHERjFLLH.EN) 

a C 669.3 

c 

NOW  CLEAP  30LIPH 

3 06660 

PPLISH(03J,1)=) 

306670 

- POi.ISH(PMn,l)  =) 

006600 

c 

GO  SE’  MEYT  POLISH  'TPINS  OPEEATOP 

036690 

GO  TO  2LJ 

C067.30 

rc 

00  72  1=1,0 

006710 

TP(P0LIEH(’, 1) . =0 .QOS (I) ) GO  TO  ru 

006720 

72 

CO  ITIMIJE 

C0S730 

C 

SQ67LQ 

■’■•T'lT  *,"  OAP  OCfSATOE” 

3C6790 

P-’-'JPM 

CC67S0 

7l» 

YFd.LE.T)  *.0  TO  2:0 

0C5770 

TPd.EO.7)  GO  *0  117 

3C6790 

’■Pd.E''.-!)  '.0  TO  At 

036790 

■’Ffl.EO.Y)  GO  gq  175 

006903 

3 CO-dO 

r 

*oT  = FT'i:''!  F"M'r'iot 

■'io9..0 

OA.L  Fi'ijr.iG  (P)L:Gi..3,  I-F'-P;  TEMPS,  F:!L,FLL  WE',  flLM.-M) 

I L 6 a 9 C 

ro  TO 

0 :6d6C 

c 

GOT  ••,”  •’''Gi.o  P6PT  0'  A 'inCTtOM 

0i;6970 

c 

•* '.•***TO^,LL  CJMC-’LC  (P0(  ISH,P)  • » 

006090 

11 

' OO-ITIMUE 

:c6890 

'A-.L  E'I'J‘;'’L  ■)(  FOLI'M  , 0,  T’F.vO;,,  e-'m=h,flLWEP,  FLLHL'M) 

336930 

00  TO  o,: 

006913 

c 

GOT  'IMA  = Y 

'.0  6-320 

12 

r,£T“**T»(  r,  c A*r) 

i,:69-.P 

■■’..L  GF'-!.'-':  {.'h*  - u ( <';1, 1 ) ,M\''E,TY<’r,LEt.) 

JCGOEO 

*•.••5  AY  c ) rviEG 

306960 

. 

-l.'Mrl 

u «.69'*0 

«co-r  (t:  -T  -v  , TL  "i 1-iE.LEM) 

006990 

'•.L  c UT  N M 'ELY, TL'f  ,Tvpc,oT') 

C':3T9C 

'■>.  i'V  ( HM”,  1 ) r3-3 

:c73)i) 

CM  TO  ’M 

C • 7 d it 

c cTP-'Tr  ( 7 .-u  , P.  TT'“e.-.,  .;T£''FM,  F ML  , F 1.  L WE  ' , 'L  L N . I M) 

;d  )23 

r 

C »T  '0  GE"  MEV'  toc^.Tce 

J li  7d  70 

?'» 

• -'ML  .'■'.It  CO  TO  .d  - 

'■:''o-.o 

♦ “S-  1 D 

■: . 7 jgo 

r.>  > 

. '7g60 

C 

■I  |T  ITI' 

• OcTlTO 

r 

-.7  VIO 

(T.  r-v  a;  GO  TO  52E 

.M-’gJO 
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'■•I'loo"  ^^,ro•^T)  gotomui 

go  ciin^T(3H  ONTIMI'^T) 

rv.L  a-UrJTI’^  (C'I'M"  .rONT'.rS) 

52-  '•TITIM'lE 
P-T'J!?N 

C H'EL  5G3  HAmT  Tg  SF7  A FL'G  0?,  ALTFcmatE  RETtR'l 

F*n 

«!'I1<»0'ITI*<F  ft  (P^LIf  H,o,ITeMOn,(-TE'1<’H,FMl.|FLLH'^..FLl.WLcgj 

IMPLICIT  I-l’’ '^‘.IS(A-7) 

FIMFF'IC  IF'/jcg  Tur  F'lMrTIO'l  NA'fE  WITH  ITS  «»GU‘*~'ir  LI5T 
(giiTLT  '^•iN':.'’Ln)  it  r GF'4E”ATES  A 'SINGLE  ASSIG'HI'IT 
CALL  FOP  THE  F'lNCTIO*:. 

HIIEH'JIOH  ptLl  T'M  -,q  , 3)  ,Ffl?GS  (2i)  , FNA  H ( 20  ) , OnjECT  ( ?) 

O’-'ENSIDH 

0’*'4  LOA^r’I/lL  ( /,  TPflEEN/lL)  /,EWAL/1/ 

nATA  F«PGE/?C'10H  / , fhAH/ZC » 13H  / 

C O-”  F'jNCTlnj  A'^G'IMrijTg 

CA'.L  GETOW'I  ( = 'iLT':H,o,p.jr) 
rF(O0LISH(3'!G,2»  . :-O.E'/-‘.l  ) GO  TO  15 

C Gf  1 ■-‘•OJP.S’Y  rPJECT  -O’.  THE  SIMILE  AP.GU.HE'ir  

f fall  GETMSHrC’OLIFHCJNO.D.EAPGE.TYPE.FLEM) 

00'.- Ms 7 

i TS-fTYOF.co.j)  GO  TO  7 

PO  is-CT  <1)  sprE  13*1 
3*s.<o»j-^rcH3H  n.  A ) : 3 COO 
GO  TO  11 

7 oO)ECT  (1)  sttf  ON 

:7r-.(p,,siT"Ho*)*i ::  ICC  j’ 

11  FALL  EMGLSFG (- 1 "GS , FLC N , T YPE , OO JECT , 0 PL^ N , FLLWEP, PLLWL ) 
‘■A7GS(l)sn.3jc3T(l) 

4?GLfMsO^LF-l 
GO  TO  17 

i IE  "ALL  GF’'IA1F(ool:Sh(ono,i),ea^OC,  -UM«^Y,aPGLEN) 

P MOH  TO  G^T  rM'l*''T'<>l  NAHE 

I"  fa^L  GET3’;''("ALTGh,cnO,fmE) 

FT.  L GftnA"'('''<'.  tF“(F!;C.l'.F.iA‘i,TYPF,PLEM) 

"A.l  lE’GSf.--:'  !,-LE‘  1) 

"•.L  HTPO.-fr';-.  t,ri  T.';,CA  GP,'.  CL^S) 

FA'.L  •••E3Gr(r.;-.  I , r lE',  , P.' V "N  , c 1 

•‘O;"''  E ^ (-'  I'’,  1 * s , 

:f(-HL.£0.1)  go  TO  ‘.f 

o'.ESs’ 

T-/ ty’F. 70. ’)  GO  TO  3: 

00>7CT  (D-sFrr  ..T.) 

oii-^TF'-'ifl'i  2 ncpo 
FO  TO  M 

-T--.r5.|-TTr  1:  j .1  'I  •02 
i.  • "OITIMHF 

•f>  I ~ yo  ~vor 

;-!GLASG(-  I.-  I.^lEN,"  >jr/'’T,jA  JF'-t.j.ti  7'|,'-|,L''7  3,cli.-ILEM) 
r . •>.»»  t 0^_0!..  2?  :-oj:''G  .ot  * 70  T-f  t'-Tii  Fii'iC ; o-j  - • • • 

r ...  is-ir.  .r  (-‘la.s  1)  --  so  2etl'e.‘' 

"'LL  1 ’ r 'T  ,0-'L  '.PTP) 

1-1  "G  I (r*;2  , ■_  ) r ’ 

GO  TO  -2 

;•  '-..L  G -TS'i- 

. g Tf  .|a  f TO'.  * '‘ML-"' , -•)  .0  ijrcr ,ru  ' 1'' .i.ojlfm) 

; M.  >T- , ■>) : .-V .AL 

"A.L  Gm'l.-.E'.  f •••I’  • .Pi.F.-. , - V 7 ^ , , j-,T,nr  iL''!,  fi.LP  ’ ’ . ^L  L ' 'I ) 

G;  "-fPl 
" 1' 

■''■•'7r-!7  - - .••'7  (■  : , .-.I*,  *y -F, 

■■'LT'.It  2i*r.--  (A-’) 


j07ino 
007110 
C07120 
GQ713C 
007140 
007160 
00717C 
007183 
0C7190 
007192 
CC7193 
007194 
3 C 7195 
CC7199 
007200 
007210 
007220 
007230 
307240 
007250 
007260 
3C7270 
007280 
007290 
007300 
007310 
007320 
3C7330 
0C7740 
307350 
307360 
CO/T’O 
907380 
00^390 

;CT.,3t 

0C7H1C 
007420 
307430 
;..74hc 
007450 
30’ *60 

2 CT^rc 

007480 
:07>.?Q 
GC’SOC 
007510 
OG7520 
.i7F-*0 
QC7540 
j 07550 
jr7';;-.0 
•;C7'373 
.C738C 

uJTS  )0 

337630 
O'.  7610 
0 .762C 
03 7530 
- :’G  »C 
307650 
;o  ’•.‘>3 
:.767C 
..  0 ’6A0 
T JT6HG 
.’700 
.3  ./’IC 
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G-rNSHE  '■.FTS  A 0=  CH!,-<ACT£kS  F^O”  THE  ftSHE 

A'JSO’.’DtN'-  T-)  PO^NTFr-',  r,I7-H  AMO  PAS'^Ej  TO  THC  ri.Li? 

THE  STFI'ir,,  PATfi  TYPE,  ANO  ITS  LFNSTH 

n'T'^HON/HAHES/  HAH  E ( IT  C , ' ) ,LSTN 
OT^ENSIOH  njTd) 

LE»=NAHr  (HPp-;,?) 

’■Y’E=HAHE  ('JPOS.T) 

HH>ns=  (LEH-l)/ia*l 
m 3 j*if  'I'.'ip-'s 
OMT(J)  » HA  HP  CPOS.l) 

N'’1S=M’0S  + 1 
'»“Tu<?r< 

E'lO 

Sm^OUri'lE  »5I'ITIT  (OUT,  LEt.') 

PotsT';  AC  'T'-iA'’ArTEP  CAPfl  images  to  TH^  N-“IT  SlHjj.ATEO 
=^0P?AH  O'JTPUT  FILE 
•AK-^S  Cft^'  or  ''CHTiMUATIOriS 
Th’LICIT  (A-7) 

'YTHENSTOH  0JT(1) 

OA'A  COHTI'C.I/LIH  ♦ / 

riH-iOS=  (LEN-1)  /IG«-1 
J=l 
<=’ 

T”(’4V|30S.LT.r)  YzHWanc 
'’■’I'lr  OAhO  I lA'-.P  Tp  cirpiJT  FILE 

P^’^NTd,-.)  {,Tj*d),I  = J,K» 

'■OTHAT  (’AlC) 

(HWPOP.LE.T)  go  TC  Zu 

pp.int  3j;r-:otHG  sarc  ihages  as  comtinueo  statph-cht 

OO  IH  isH,  HH’)<;,  *, 

•’■K.GT.HHEn')  |<  = MW-<0E 
pot'iT  (l,E)  (OUT(I),  I = J,<) 

J 1=  J 

1.  rO'ITIH'J- 

■<  = <-*, 

'LEV-''  OH-  1 --  ?£  p::hTE0  TM  SHFE'ETIn'-.  sLrvlTE 

■;}"  ■•3PA ''-liM  0','E  CAP,' 

?Z  OT  -H  ^ LL  - JJ. '' 

'UP'I 
P'H 

<;if>’e.o,,7T,jb  pii;-;  Pl. -1  ( r y , p,  tphpm ; r jpf'-o  , FLLH.- FLLHL  EH) 

’■‘'•LT'TT  T-  T^Apo  (ft.T) 

'^'JMC-o.o  M-'-,-"  Pi'H''"'or:  •.  P'.'j-.P  jT'  TO". pp .gEM"^ d t' HO 

O'lL;  =■00  TT'I'-L'  p' M '■■p;.T  " I?  fiF  C E A "Y  . 

o-tE'lo:il  "Ll..'EOfi) 

O-'EHSIOM  r , r A"E  (O  >HEZ  t’C) 

o’.'A  ■A'/A./',/,  'o'i:/iit,/ 

0";  / , lAMF’/ZC*  1 )H  / 

"'.PHI  A'lo  PL"M?  Or'-OCS'HT  pr  T E” <"  1 T APY  HA-*E 

•'.'■r.  pTL 

O'-.L  '."ro-  ' 'L;'  o,  •;  '< 

-r  >r-.|  r;  o.. . ) Tnr- >•!  T S " PP'’AT‘p  A '^TIGL" 

-.-r '-M  " •>  ' " *.■;.  )■  -'(f  (':o  tt  r T'^ 

)H  Thp  '.'-07-  o • ->: TPP P fj  ’p ) , 
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C ftMO  C£'|-r  ^vnt.'X  IT  FILL'^  Vcv4InI‘I0  =0^TT0N  .13514 

C T*  Tmf  I. ■'me  '^T'M  BLlMk":  (<!Hr=-TIK'''.  T“£  '^'’ST  nouN  TEN  3’ICES)  C 13615 

n’-»ENSI0N  OJT(l)  u 13520 

OT'A  ■■•L&'K'/I  :h  / 01C630 

’■'■(LE'I.LE.b;)  OO  to  o?  3105^ 

T’'aEM.GT.Fj.-i-|''.Lf‘J.LT.7:)  50  TO  3 010650 

TF{MOO(.ENo3)  .L-'.ID)  GO  TO  5 010550 

GO  TO  03  010670 

5 ri'.L  NE?'5f  <D'JT,LEN,9LAKk'S,lu)  010530 

qq  O-TUON  CIOG-JO 

“^NT  010700 
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N-bit  Simulation  Subroutines 


Twenty-two  n-bit  simulation  subroutines  were  built  to 
handle  each  combination  of  data  types  and  arithmetic  operations. 
In  addition,  there  were  four  other  subroutines  developed 
(SETNBIT,  ROUNDER,  OVERFLW,  and  UNDRPLW) . 

The  SETNBIT  subroutine  was  developed  to  compute  seven 
key  values  upon  which  the  other  n-bit  subroutines  base  their 
n-bit  wordlength  simulation  effects.  The  OVERFLW  and  UNDRFLW 
subroutines  are  used  to  print  overflow  messages.  A subrou- 
tine called  ROUNDER  is  used  whenever  a special  case  results 
where  a round  up  increments  the  exponent  of  the  floating 
point  value  (due  to  the  fixed  point  add  performed  to  round 
floating  point  values) . 

The  22  n-bit  simulation  subroutines  basically  are  the 
saune.  They  can  be  divided  into  two  classes:  those  which 
produce  real  values  and  those  which  produce  integer  values. 
After  the  operation  is  performed  in  each  subroutine,  the 
rest  of  the  routine  is  eventually  the  same  as  each  of  the 
others  within  its  own  class. 

The  source  listing  of  the  n-bit  simulation  subroutines 


follow^: 


1 


\ 

! 


N-bit  Simulation  Subroutines 


Source  Listing 


UOOUU  OO  (J  O O U'C 


^'IIROUTINf  SF'NgiKNaiTJ.IRHOTRjIPRCSN.MESSGSilANTSi.tiXPNT, 

• IPTOOS.K'Y) 

THIS  ’OUTIHc  '^STA'’l.ISHS  THE  MIN/H/IX  VALUES  THAT  'lOULO  SE 
REORESENTET  IM  FIXEC  ANP  PLOATIfT,  «T  fS  OF  "^HE  SI'/EN  N-BIT 
SoECIFICATIoHS-^OI’TINE  TERHIt'ATCS  PROGRAH  FOR  ILL 

nUENSIOM  KFVC’) 
flATA  HANFL'’S/4A/ 

Iff  (H4‘jTSA*t'’XPnr*H  .nE.NSITS)  GO  TO  31 
TPtNRITS.G'T.EJ)  GO  TO  30 
▼c<I9H0TR.LT.3)  GO  TP  32 

«»**»*SOFE  Tucsi^  "A'lT  USE  REGULAR  UNITS— OEPEUOS  ON  Hit 


OCOlOO 
OOOtlG 
000120 
000130 
00 0140 
OCOISO 
3C0160 
0Q017C 
000190 
000190 
000230 
003210 
PROGRAN****000220 


MIGHT  hahT  INTEGE“  -ax  to  SE  MS1TSt2  TO  NAXE  UP  FOR  SIGN  PIT  000230 

XEV  (1)  sSHI^r  ('•A<^x  fl)  .NPITS)  -1  0 00240 

TP(IoRCSN.r.r.T.0P.l"Rr.5‘;.LT.l)  go  to  33  000250 

I«’f  lEXONT.GT.ll)  GO  TP  75  000250 

TPC-AN’'SA.''Q.NAMTLEH.ANr.iR»,rTr,,„r,p,  gg  TO  34  003270 

MAXMANTi-AMTSA*  {f’RTGfi-l)  •uniTS  000290 

IsfHAX-ANT.GT.NAMTLFN)  GO  TO  34  000290 

xe-V(4)  sIRrOTR  000300 

S^T  mESSAGc  o7T‘r/SUPPFESSI0N  FLAG  000310 

XEX(5) =NESSGS  000320 

<'^V(9)=NAM''LFN--ANTSA  000330 

x«T(7) sHAntlEN-hayhANT  000340 

•••‘••NEEO  50  SI’  iLLTWfJ'CF  IN  THIS  PME  AS  OF  17  Nnv-MRER*»*»  3C3350 

TF(IEX3NT.ro.ti.ANP.HANTSA.ET.43)  GO  TP  25  000360 

TRNGE=shTFT(mi^V( « ) ,IEXFNT)-1  000370 

PIG*!.  •»TRN'.E  303393 

TFtloTPOS.FC.O)  ''IG=2.*»(IRMGE-HAUT5A)  003390 

T'»-SY  = SHrFT(“AS  <(''ANTSA)  .MANTSA)  0 004CO 

pn1X=IpNAy»oig  C00410 

onAYsP-AX.O®.  (SHT-T  (HAGif  (mantlEH)  .hANTLEN)  » 003420 

'»MivrGuTc»(?uir'-(^M*Y,-ifEY(7)  ) ,XEY(7)  ) 3C3430 

SHTc-T^'iG  ■>41Y  '’ISA’LFS  the  HOOE  fCNVFPSION  TO  TNT'SFR  FOR  < ET  ( 2 T C 3 ‘•'•(J 
<‘T  ( 2) -Sh:c- (7-1AY,  3)  003453 

H'.YHiCyTr-  (-40Y  (1  ) , ‘'AMTEA  ) 003460 

R'>’  = H«Y-  333470 

'•HTNs'’P5»  7.  » • ( -jp  >;gF)  3C3483 

TCfiDTooS.F-'.-')  =“:N  = FPR*2.»*(-rRr'GF-HANTSA)  010490 

YCy(Y).s<;mtct(3.<X.|,C)  1 000530 

»Y  o'-turn  ■ 0C3S13 

SOETIA).  CT>  G,’  r.jT  V'ORO  HAY, ITS  P"AL  -AX  WOl.'LP  S0N9  Thf  JOS  303573 

mhEN  any  70  1PAPTS  ilfPF  -APEt  SO  ITS  hAO-  SLIGHTLY  LESS  300530 

25  PTtTTN'JF  3 C 3540 


(?) 

i?Y7G7?7  7"T7  77'7T7  779 

0:3553 

Yrv  (3, 

:C03i-*3:3'33::3-:o-0'0“ 

OCJSSC 

TO 

JT 

CC0G70 

• ’•riT  -tAuY  PITS  SCPGICTPP;  J'-IIT  HAY  IS  53  ", 

CC35.A0 

ft 

••  — "tMSiTs,"  vc’E  s=’‘':if:':'3~ 

sc:  JPG 

r-.fy  TO 

4 0 

0306:3 

’1 

-VjTjtg:  ♦ CY'C’FNT  ♦ 1 S|i:-"Lr  foijal  'N'-’tTS" 

C 3 f 6 0 

to 

-0 

330723 

■»  ^ 

’LfS'-L  ■ '"‘.''•vG/TRUNPATTP-;  0=TI0:!  VALUE  ", 

: C3670 

ft 

••  '’7JM3T  M-,  I 1/  •■S'.nCATON  * u" 

C'JS'H) 

on  Ta 

40 

C : 3 950 

eAr*^r 

• , ^LL-*AL  A •pl'HT  pt  PPPPI^'P';  SPF'>»pTro  “ 

::3663 

03:6Yo 

r.A  TO 

-1 

: 

•*•5  ^ fcj  ^ 

, ^ 3*Tc  ■5nrrT-“Tco  **gp  ( — 'IT^  * T A * p R I<  f 3^;)  *• 

•3 : OGpr, 

ft 

••  Til,.  m,vTM'Isi  I,  , *1 

161 

*1  r rt  T , a 

«>-N  • < t.  IV  I I I Ml  * • IMM  4 ■ I I IWWV***  < 1/  M I*  t M 


F 


35 

bO 


B 

loa 


c 


4 


<» 


?6 

P 

?5 


C 

31 


m 

r 


It  ’ 


■*r, 

33 


r 

■tr 


RO  TO  40 

P’TMT  FKOO'Ir'iT  TOO  lAPGF  FOR  CDC--  * 11  =3TTS~ 

ORTHT  FR'^03-  tlLFOAI  ^ETM^it  DATA  PARAIETFR  FO'J'IO" 

^Toa  ~ q«n  INOIJ-'  *0  SETfRIT  K0UTIHF“ 


FM1 

TMTER^R  'UNCTTOM  TIAOO  (11,12, ROUTINFtLMPNT, IFMt, KEY) 
OTHFNsION  KFY(7) 

TY>00=ri*I2 

Y'^dARSdlAOO)  .RT.KFYd))  00  TO  100 
PCTURN 

IP(KEV(5) .NE.O)  CALL  OVERFLH ( ROUTINE ,LNCNT ,8HAOOITION» 

Tceitino.GT.O)  TTATn=YEY<l) 

YftllAOO.LT.':)  TIAr)0  = -5EY(l> 

P.O  TO  5 
FHO 

RtHROMTlNE  0V-R'’LM(O0UTTNF,LNCNT,0PFR) 

IMPLICIT  INTFOER  (A-Z) 

•••*•  MAY  WANT  TO  PHANOE  THIS  PRIN'^  FROM  A PORMATTEO  •••♦ 


Pp:MT(1,4)  R0'ITT'|p,LMCN’,0P‘'R 

ro5MAT(-  .•••  ovEpflOM  in  ".Ae,**  LINE  • ",I5,"  FOP. 

RPTURM 

FHO 

S'nROUTIMF  UNORPLN (ROUTINE. LNCNT, OPER) 

TMOLtriT  TMTFREP  («-7) 


OPTNT  (1,4)  PO!Jtt«)E,LNPHT  ,0PER 

Pn?MAT(~  UMOERFLOH  IN  “,Aa,”  LINE  »",I5,"  FOR  ",A10,"  *••“) 

RETURN 

ENO 

PPIL  function  R01?0(R1  ,f2, routine, Lt.'CNTjIFNL, FEY) 

PHENSION  YEY(7) 

PRJOO^Rl ♦P? 

TP(i<EY(4)  .NE.3)  00  TP  30 

IP|A3S(RRA00)  .p,t.<:mIFT(kfy(3)  ,0))  GO'O  100 

IRfAOS  (PRAOO)  .L  •'.EHTFT  (FEYd)  ,0) . A.ORAOP.NE.O)  GO  TO  110 

riirri<  pgo  poUMOruO 

R-iAOOsRHTF-CRHTFT  ('’PA00,-KEY(T)  ) ,kEY(7)  ) 

TF(IFML.EO.l)  °'»aOOsSHTFT(SHIFT(RRAOO,-KEY(6)) ,<EY(s)) 

RRTUR‘1 
P.O  RO'1‘10 
TM'*l 

TPfaJin-t.LT.o.)  TNC=-1 

TFT  ( P'J'  r'  (S'^IPT  (PTfiOO,-  (vcv  (7)  ..  ) ) ) , ■<-v  (7  ) ) 

T=’'.3u:pt(<;mxRT(op*0'',12>  ,-12)  .EO.O)  CALL  OOI'NOE^  ( RR '■  DP) 

PO  TO  24 
-P(=o.\OI.GT.O.) 

TP(.3t>*,'1,O.l,T.0.) 
rucr/  FTP  .iFOF.lOES 
TPfKEY (I) .ME. j) 

P.O  TO  ?■! 

RO) 00=0 

TP(<EY  (?)  .‘IP,.-  ) 

PO  TO  33 

P‘|0 

OF-.L  PU-;PTTO.|  33a-P(7l,T2,R0UTINF,LrCf''',IFNL,YEY) 

'-..-•l.'XOu  <rv(T, 

c-lOOs  M *17 

Tr(w-Fy  ))  .-.Q  TP  tq 

r’(c 'vS  (''■’. *,po)  .p,T. 'htt' (' rY( ')  ,0) ) p-o  to  loc 

T-f  3 3-' roi!  OP)  '■'PT  (iTFy  TT)  , 7)  .0  .P2  •O"  .MF.IJ)  go  TO  HO 

•>7-.U  1-.P.4T.-*(C  JT"  ( •,0O  .-kP  M .M  ) , EE''  (7)  ) 

T’-'  TFOL.  'o.  1)  •2'0'*phIFT(ShIF  T(P2A  J';,-<EY  (4) ) .CT  ( 3 ) > 

o-'-l  J-l 


o=Ano=SHiFr CFEY(2) ,0) 

SP^OPs-PHiFT (KFY (2) ,0) 

Cf  LL  0VEPFLW(P0UTIN'^,  LNPMT,3U'.00IT  ION) 
CfLL  U‘;O?FLN("0UTINF,LHruT,,«HJOOITT0N) 


'•j'  r 

T-<  o-. ,1  * , 7 , ) IPPP-I 


ceoTio 

0C0720 

000730 

000740 

0C0750 

000760 

C30770 

0C0780 

0007RO 

COOSOO 

ccasio 

3C002C 

000830 

000840 

1)00850 

000860 

00087C 

500880 

0C08R0 

OOBROC 

OCORIO 

000420 

000430 

000440 

000450 

000960 

000470 

000400 

C00440 

001000 

001010 

0C102F 

C01C30 

001040 
001050 
0C1360 
0C1370 
001380 
001390 
GO1103 
351110 
351120 
C C1130 
0Cli4C 
3C115C 
001160 
501170 
0C1130 
001140 
3C1?t)0 
JCIEIO 
001?20 
001230 
501240 

orizsc 

151260 
: 5 1270 
. C1700 
0C124C 
Q&l’CO 
351310 
C 01 320 
5 OISTO 
'.Cl  340 
cot  350 


^62 


131 

110 


zr> 

?8 

C 

30 


113 


110 


133 


1 

103 


«%0  \ 9 $ ^ I%»  $ • 

TT(SHI'fT<«;mrTj<??ano,12)  ,-12)  .EO.O)  CSLL  ROI!NOERC<2»On) 
r.o  TO  ?5 

I='(‘>T1'?3.GT.0.1  ®2«''P=SHrPT(WEY(2)  ,C) 

▼•■(REAnO.LT.O.)  R2rOD=-S‘<TFT  (K'^Y  (2)  ,C) 

TP(<EY(5)  .Ne^.l)  CrtL  OVERFLH(«>OUTINE,LMCMT,HnniTION» 

f-n  TO  23 

P?lOO=3 


ir(KEY(5) .NE.O)  C^LL  UN0RFLH(P0UTIl!E,LNCMT,8Hi03rTI0N» 

GO  TO  28 

EHO 

RE»L  P'lHCTION  Rift  00  tl  1 ,P2  .ROUTINE, LNCNT,  IF‘IL  ,yEY) 
nTMFNEION  K'^Yl?) 

»H003ll*P2 

T"’('<EY  (4)  .MF.3)  r,0  TO  JO 

IP(ftOE<Rlflm)  =;mTpT('<’EY(2),0))  go  to  100 
I'  t&3S(Rli'”»  .L’’.SKTFT(KFY(3),C)  .fl.PllOC.NE.C)  33  TO  IIC 

oisnO=GHTrT(P'<IPT  (oii)'ip,-i<EY(7)  ),  KEY(T)  ) 

!«■(  IFNL.pp.I)  R1400=SPIFT<SHIFT  (R1A''3,-KEY(E)  I ,KEY<G1  ) 
Rptu9>| 


GO  ROUNO  IT  NON 


TSG»1 

tp(R140D.LT.O  .)  I'!C  = -1 


R1100=G*<I'”(GMTP’f  S«TFT(01flOO,-(KEY(7)-l))  ••It!'"  ,-l » , <EY  ( 7 ) ) 
I'’{SMIPT(SMlcT(ouon,i2)  ,-12)  .EQ.O)  CALL  ROUNDER  (Rl  A 00) 

GO  TO  25 


TPIRRftOO.G'^.O.) 
T^lROftOO.LT.O.) 
TP(KEY (5) .UE.1) 
GO  TO  28 
PllOO=0 

I^KEY  (5)  .NE.3) 
GO  TO  28 


PRftOO>SHIFT(KEY(2) ,0) 
roapra-sHIPT (7EY (2)  ,C  ) 

GALL  OVERFLH (ROUTINE, UNCMT ,8HA00ITI3N) 


CALL  UN0.RFLH<P0U7INE,LN0NT,8HA00ITI3N) 


PUO 


I‘I’’FGER  FUMO'^TON  ri““Y(Il,I2,R0IITINF,LNCNT,IFNL,KEY) 
OT>icnGION  <EY(7) 

77soysii»i2 

T'^dftOSl  II^PY)  .r.T.KEYd))  GO  TO  100 
PET'JRN 


GAUL  0V"'’'’L'(  (-•.UT»“IE.L”GNT,3HMULTIPLY) 
T<r(Iti'>Y  .r.T. 'I  i:«OYi<EY(l) 

’•'•(IIMT/.L''.  IT'iPYs-KEYd) 

GO  ,TO  G 

Cijll 


T'lTFr,^-?  F'IMOTIO‘1  ri«NS  Cl  ,I2,'’0UTIUF,LNGNT,  IFNL.KEY) 
"’■'■NStON  '(EY<T) 

I’-»NG=I1-T?  . 

’■'■(IftOSdIM'IG)  .GT.YFYd))  GO  TO  100 

F'^7U9«| 

■''•KEY  (5)  .‘I- .3)  '•all  OVPTFLM(ROUTIUE,L*l(;NT,6MSUTTRft!3T  ) 

TT(ir''«!G.G‘'.C)  TP'UFri^FVd) 

T»-(  jtm.iE.lt. ;>  'T'1N‘’i-K£Y  d) 

GO  TO  0 
T*n 

TM- cr.r^  cM.jEr  jp.j  T J Ov,f)  ( J j , jj  ^ j jfip  . NT  , IFML  , K'’Y ) 
'’T<E%GT'’'l  ^ty(-) 

”070=  Ji  / T7 

-■f  tftOTCrji.in,  ,r,-,  ^-vd)  ) GO  TG  IPf 

OT'IIS.) 

r rvTT|_..^(...-,,T  jf,r,Lfj<-.jT,  *HOrvTE:0-j) 
:''(if''/o.G'. '’)  T!''/G-/rY(i) 

f'cin /G.i  T,^,  riGYOs.i«’'-Yd) 

"0  TG  9 


"NT 

J.,rc-r-j  r.i.;G-jG>i  d'" 'R  (■'1  ,I'», ’CUT  INr.LNFNr  , IFMi.,KFy) 

T • 1 r ,w  1 I 

»^rvO  - Tf  ..TT 


001370 
001780 
G013R0 
001400 
C0141C 
001420 
0C1430 
001440 
001450 
001450 
001470 
001480 
001490 
001500 
001510 
0G1520 
001530 
C01540 
001550 
OC1560 
001570 
001580 
CC1590 
CC1600 
001510 
001620 
001530 
001540 
001650 
001650 
3C1670 
001580 
001690 
001700 
CC1710 
e"i72e 
001730 
301740 
3 OITGC 
001750 
0C1770 
;01780 
301790 
001300 
001810 
0C1820 
001830 
001840 
001850 
031850 
CC1J70 
001880 
30139C 
OCIHO 

oonio 

c:i9?o 

001940 

uCl')50 

031950 

001970 

OulPAO 

031990 

O'^JOJ 

0':2:i3 

•5  89  » 
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T'fdA'JSdrpyo)  .r.T.KfYd) » co  to  lod 

5 »PTUR*I 

110  !=■( <*’Y  «S)  .N^.O)  CSLL  0VE'’«rLH(P0UTIH<i,L''ICNT,8HEXP3'»E'JT> 

TPtliexP.rtT.O)  riEXP=KFY(l) 

IPdlEXO.LP.O)  TIEXP»-KEY(1) 

RO  TO  5 
P>n 

"PIL  PUMCTTON  PPIPY (P1,P2,pOOTINF,LNCMT,IF'».»PEY) 

OT'^ENSIOS  <EYC'> 

PP'4PY3Pl*o? 

TPtXEVC*)  .M- .0)  CO  TC  JO 
?1  (9P«PY)  .CT.CMIPT  (I'CYC?)  ,0)»  GOTO  100 

T'?(4ni{o»HPy)  ,LT.tmtFT(yEY(3),0)  .A.t»R“PY.ME.O)  GO  TO  110 
28  »»'<OY*S«TPT(>;-^IF'f  (RPT^Y, -K^Y  C')  ) , XEYC*)  ) 

TFdFNL.EO.l)  “OYsSHIFT (SHIFT (ORPPY, -KEY  (G))  ,KEY(51) 

RETURN 

C GO  ROUNO 

33  THGal 

f^fRRHPY.LT.O.)  T*lCs-l 

R?'1PYaSMIFT(';Hie'T  (SHIPT<p'>*10Y,-(KEY(7)-1))  ♦IMC , -1) , K- Y ( 7 > ) 
TP{SHlFT(P«IPf(''T‘"’Y,l2)  . -12)  .CO.O)  CALL  ROUNOE?  ( M'-Y) 

GO  TO  2> 

no  I''(RR“PY.rT.0.>  FP‘'PY=SHIFT(<EV<2)  ,C) 

TC(poHpv.L’'.0.)  pRmPY=-SHIct(KEY<2),0) 

C OH“C<  FO®  HFSSIC'^S 

TF(KEY (5) .MF.3>  CALL  OYEPFLM (ROUTINE, L'lCNT , BHMULTIPLY ) 

GO  TO  28 
no  P-»tPY»0 

T®(KET(5)  .HF.l)  CALL  UNObF(.W  (ROUTI NF  , LMCNT  , IHMULTI  OUY  ) 

CO  TO  28 

p>p 

"11.  FIInJTTOM  “ommS(R1,'‘2,®OUTINE,I.NCNT,IB'IL,KEY) 

OTYEMSTON  KFY(7» 
pbhnS*R1-R2 

TFfKFYd.)  .H-.-')  GO  TO  30 

2G  TP(AnS  (Romm^;)  .'-.T.'^MTcr  (KFY(2) , 0) ) GO  TO  no 

TP(AO«  (RRhmc;)  .lT.phIPT  (YETCJ)  . 0)  .a.BOHNR.NT.a)  GOTO  HO 

28  p-»...pa<;mrr  (PMt  p-  ( p - yry  (7)  ) , KFY  (7)  ) 

Tb(IP*'L.F0.1)  (S'iift  (RPMriS, -Cy  (E)  ) ,<'’Yf=i)  ) 

O'TIIRN 

C GO  RO'IMO 

TO  T-r  = t 

»^fooM.j5.LY.0.)  T‘iC  = -l 

o■»<•|«;sP^t'■T  ( phytt  ($UJCT(PP'V|$  ,-  (<  EY  (T)-l ) t yING,  -i»  . <-Y{  7)  ) 
TP(S>'tFT  (SHIFT  (=R'NF,12)  ,-12)  .PO.  0)  CALL  ■’O'JHOFR  ( RRYMS ) 

CO  TO  .28 

no  IB(»P<NS’.CT.P.)  OR.>iNP»«MTFT(<rY(?)  ,3) 

Tr(POv,.,5,LT.i: .)  RRH?  Ss-GHTFT(<FY(2)  ,0) 

r rHTK  FO’  “FOGAGr- 

r‘(KTY  (5)  .MC.  J)  till  0VFbFi.h(R01ITI‘|F,1.:,''NT,5*-SU0T9ACT) 

GO  TO  28 
no 

tt»<FY  (G)  ,*IB.-)  CALL  tJf;ORPLH(»r'UTtNF,L*l''NT,YH<;jaT  = AC*) 

GO  TO  28 
*■•10 

")L  bumCTTo*'  T-'  .-.(■.l,''2,T0'jTrMF,t.»lCflT,I''r;L,<tY) 

‘'Z  jA  TO  ! y*  (T  ) 
p-* 

TF,:<rY  (j.)  ,.;r,  <)  r-  ,3 

3-  r-t  A Os  ('”-'''). -.T, go  n n; 

TPf  AOS  (■••>''■(")  .L I T)  .*).' ..'  js'o -,3  TO  in 

2«  ps  ../’’iPu-r-c?  i;  ft  <T  ( 7)  ) 

TB,': :'.j ) * •.  -i  *'".  ' '■’i — -y;-' ( i)  1 :i ) 

P-TmH 


002130 
0 020>»0 
CC2030 
002060 
002070 
002080 
002090 
002100 
002110 
002120 
0C21J0 
002140 
CS2190 
002160 
002170 
a(.218C 
002190 
002200 
002210 
CC2220 
002210 
002240 
002290 
QC226C 
002270 
002280 
002290 
002300 
C0231C 
002320 
002330 
002340 
002390 
002390 
0P23T0 
002380 
C02393 
102s03 
0C2V10 

GC2431) 
1 C 2'*40 
uC2',G9 

1 2 2<*oC 
1 G?470 
OGPsrtO 
(i02-*90 
:r25Q0 

oc2‘;io 

002920 
C :2933 
002940 
JC2990 
.10  2960 
':C2973 
C129.n 

IJ  "2910 
3c:’'.:3 
0 "2912 
'C2920 
0:29T0 

n2:.-o 
: 0 ’If- ) 

■ i"  ."’97'' 
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Tf(RROVO.LT.n.)  TNC=-1 

SHIFT  (CMTCT  (SHIFT  (oo070.-(<EY  (7)-H)  TlNCt  -H  , <EY  ( 7 ) ) 
IF(SHIFT(SHIcT(:?-3rivn,ie)  ,-12)  .£0.0)  CAUL  ROUNOER(R?O\/0) 
on  TO  26 

r VO. GT.O.)  PP^VOsSHIFTtKEY  <2>,0) 

rF(RR':'V''.LT.C.>  PRUvr  = -SHIFT  (KEY<2),0) 


C/LU  UHORFLHCPOUTIHFiLNCNT.aHOIVISION) 


re’tpPrjvo.GT.O.)  PP^VOsSHIFTtKEY  <2>,0) 

rF(RR':'V''.LT.C.>  PPUvr  = -SHIFT  (KEY<2),0) 

CHFCK  FOR  M''ssar.-s 

IF(KFY  (5)  .N'^.a)  CALL  0 VERFLH  (ROOTINF,  LNCMT .aHOI'/ISION) 

GO  TO  28 
RPOVOaO 

IF(KFY(5) .H'.O)  C/LL  UHORFLH (POUTIHE i LNCNT, SHOIVIS ION) 

GO  TO  28 

PNO 

REAL  FUNCTION  R2HNS (91  iI2 iRPUTINE ,LNCMT , IFNL,KEY) 

OtiENSION  KFY(7) 

P?HNS=R1-I? 

CHFC^  FO®  PO'INOING 
TFKEY  (l»)  .HE. 3)  GO  TO  70 

IF(A<;S(P2«HP)  .GT.Ch:ft(>.'FY{2)  ,0))  go  to  150 

tF(aOS  (R’-HE)  .L  y.Gh-'itt  (i'EY(3)  ,0)  .a.P2YHS.He.O)  GO  TO  110 

“2HNSa  SHIFT  (SHTc-t  (o  , -KE  Y ( 7 ) ) , KE''(7)  ) 

IFfIFNL.EO.l)  P2HHS=EHIFT(SHlFT(R2MNS,-<eY(6)) ,<17(5)  ) 

fftuRN 

GO  ROUND  IT  NOH 

YNC  = 1 

rF(R2MMS.LT.0.)  INr=-l 

R»h»jS»SHIFT(  SHTF’-  (SHIFT  (P’mnS,- (KEY  (7)-l))  f INF  , -1)  , <Ey  ( 7 ) ) 
TF(SHIFT (SHIFT (03mne,12) , -12) .EQ.O)  CALL  P0UN0E9(92MNS) 

GO  TO  23 

TF(92HNS.C.T  .C  .)  92ri»'S=SHIFT(KEY  (2)  ,0) 

IF(R2HNS.LT .0.)  o2hnS=-SHIFT (KEY (2) , 0) 

TFIKEY  (5)  .►!=■.))  C/LL  OVEFFL'/l(POUTTNE,LNrNT,3HSinTRaCT) 


TF(92HNS.C.T  .C  .)  92ri»'S=SHIFT(KEY  (2)  ,0) 

IF(R2HNS.LT .0.)  o2hnS=-SHIFT (KEY (2) , 0) 

TFIKEY  (5)  .►!=■.))  C/LL  OVEFFL'/l(POUTTNE,LNrNT,3HSinTRaC1 

GO  TO  29 
REHNSsO 

T‘^(KFY  (5)  .NT.O)  call  UN0PFLM(F0UTINE,LNCHT,A.H3J3TRACT) 

GO  TO  29 
ENO 

“'•IL  fijmotion  91‘'HS(Il,F2tROUTINE,LMCNT,IrNL,KFY) 

''’■yfnSIOH  YFYC*) 

®lTHS=Tt 

OHcr<  cr)r  ^ruiniIHG  rrHON 

TF(KEY (U) .yr. ])  GO  rc  TO 

TT(  AOS  (P1-,S)  . -T.  TVI-CT  C-'CV  (2)  ,0)  ) GO  TO  irg 

TF  (AOS(PIH'IS)  .LT.ShIFT(»'^-Y(3)  .0)  .A.RlHNS.NF.n)  CD  TO  11 
■»l.''HS  = ^HIFr  ( GHT  TT  ( PI-*'-:,  - KEY  (T)  ) , '<ry  ( 7)  ) 

TFdFML.'FO.l)  PlH>;3-=''VIFT($HIFr  (R1M':S,-KEY(6)  ) ,<Fr  (i)  ) 
orrj|9„ 

GO  POHHO  IT  ?10H 
IH'  = 1 

TT(“1'<H«.LT.ti  thc  = -1 

ofHSiEHIP’-  ('■  -t'T  - ( - HIFT  (p1T';S,-(<EY  (7)-1))  ♦INC,-!)  , <ty  (/)  ) 
T=-(SuTCT('-.iirT('’l“‘ls,l/'  ,-l?)  .'0.0  CILL  ROU/.C' V ( ? 1 “HS ) 

GO  TO  23 

TT»  -)1  1.;;  , P 1 ■■'•’S':E'‘IPT  (KEY  (2)  » 3) 

’"(Rt-HS.Lr.O.)  '1  =:-E''IFT(<Tv  (.->),:) 

-r  (1/ .TY  (5  ) .'IP  . ) GAEL  O'.’FGCLH  ('•''H'’  IH-:  , LHONT  , JH<?:nT^  IGT  ) 


33  TO  110 


TT» ->1  i;  . , 7 I-, 

’"(Rt-HS.Lr.O.)  '1 

*P(*'-Y  (A)  .'IP.  ';) 

GO  TO  29 
p«*.4S-(5 

*F(  KFV  (S)  .‘IF  , 3) 

GO  TO  23 
F'lO 

-.-AE  Tl  .i-TTS”  pTlTYi 
O'  'EM-’ : )•!  /'Yf  ’i 
*1  ’ 


'.'.LL  U:i:PFLy(PO')GIHE,LH:‘iT,  ’“SJOTRACT) 


' . yvJTlME  ,LVG‘;T,  TEHL  ,‘'£Y) 


00260C 
002700 
002710 
002720 
002730 
002740 
0C2750 
002760 
002770 
002780 
002790 
G0280Q 
002810 
002020 
CC283C 
002840 
002850 
002860 
002370 
002980 
002890 
002900 
002910 
002920 
002930 
002940 
002950 
C 02960 
002970 
002930 
002990 
003050 
003019 
003920 
C C3530 
003040 
003050 
003060 
CC307C 
•3  03  09Q 
003095 
003190 
003110 
503120 
0 031.30 
003140 
053150 
QC3165 
003170 
003130 
0C319C 
3C3200 
00.3210 
033220 
OC  3230 
0 0 321.0 
)G32sp 
0C326C 
00327: 
C 33250 
003299 
j 0 33('3 

OOOTIO 

053320 

3''T3TC 
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6 

'>W  1 w w 

w «l  ^ 

I<=’(48S(^2‘<'"')  ,r,T. 

SHTFT(>'fy(2)  ,0))  GO  TO  IOC 

OGJTGB 

Ic»aRS (a?Moy) .lt. 

SHTFT{FFY(I),0).A.P2!'3Y.MF.a)  GO  TO  110 

Q0336Q 

28 

o^MOys^Hie^T  (3HIcT(o?fPv,-KFY(7)  ) , <EY<7)  ) 

003370 

T'dFNL,  P'3.1) 

P?hpy=Shtpt(SHIFT (P2MOY,-KEY(6)) ,<;Y(5) ) 

aC3360 

RPTuhh 

003390 

C 

GO  ROUN*)  t7  NOW 

003430 

33 

TN0»1 

aC341C 

TP(R2NPY.LT.P.) 

rNc=-i 

DC3420 

«>2NPYsG^tPT(  <:-4TPT(G>IIFT  fOa-PY  , - (KFY  ( 7 ) -1 ) ) ♦ INC  , -1 ) . <FY  ( 7 ) ) 

003430 

TPCSHIFTIPHIFT (P’-Pv,!?) ,-ie> ,EO.O)  CALL  ROUNOFO (PeN^Y ) 

0C3>»40 

GO  TO  25 

003480 

103 

IPfpaMPY.GT.O.) 

o?NPY=SHIFT(KEY(?) ,3) 

003480 

IC(pjm"V.LT.P.) 

P2MPYT-?HrFT(<EY(2)  ,0) 

003470 

TP»KFYt3) .NF.D) 

CALL  CVERFLW (ROUT INE.LNCNT.SHMULTr PLY) 

003480 

GO  TO  28 

0C3490 

111 

R?1PY=0 

003500 

tP(KEY(5) .MC.a) 

CALL  UNOP,FLW(POUTIH",LHCNT,3HHtJLriPLY) 

0C3810 

GO  TO  28 

303820 

FNO 

003530 

RFIL  FUN.3TT0N  RlN’Y(Il,F2,R0UTINP,LNCNT,IFNL,t(FY) 

003540 

OIMENSION 

C 03580 

Rl'iRY=Il*02 

003580 

C 

OMFCK  FOR  “OONOING 

0C35TC 

TF(YEY  (1.)  .NE.  3) 

GO  TO  30 

003580 

26 

I”14'»S  (RIHDV)  .-.T. 

GHTf-{fFy(7),0)>  go  to  ICO 

003590 

TF  (APFIPIMOY) .lt 

• SHTFTtKFY  (3)  ,0)  .A  .Rl'IPY.NE.  C)  10  TO  110 

GG3600 

28 

oiWOYs SWIFT  < FWI FT 

< PIMOY , -key (7) ) , KFY ( T) ) 

003510 

TF(IFNL.“''.li  Rlwpv  = $wlFT(SHlFT  (RIHPY,  -KEY  (F)  ) fYSY  (5)  ) 

03352! 

PFTURN 

003830 

r 

GO  RO'JNO  IT  NOW 

003640 

30 

INOsl 

003650 

IF(R1npy.LT.O.) 

INC*-1 

00 3660 

oiwpY=$NIFT(S'!TFTfGHIFT(l?lwPY,-(KFY(7)-l))  flNC.-l)  ,YrY(7)  ) 

00  3670 

TF(SWIFr(SHTFr(PlNPv,i2) ,-12) .EO.O)  CALL  ROl'NOFR ( P 1 NOy ) 

003690 

GO  TO  25 

003690 

111 

’■FlRlNPy.GT.n.) 

oiMPYsKMIFTCKFYfa) ,)) 

CC3700 

Ic(s.mpy.lt.O.) 

Pl“OVr -shift (KFY (2) ,0) 

CC3710 

TF  (KEY(5) .NF.O  ) 

CALL  CVEPFLWIoOUTINE.LNCNTjFHNlJLTrPLY) 

003720 

GT  TO  29 

0C3730 

113 

oi'«py=  a 

003740 

T"(<'y  (3) 

CALL  ''NORFLW (PONT’NC.lNCNT , ? HMuLI ! =LY ) 

: :3TG0 

O'}.  TO  79 

013760 

F‘n 

0C3770 

“■IL  fu'ICTION  ?70V0(Pl,I2,R0'JTI:!F,L):C>ir,rFNLt'''EY) 

aC379C 

"T-lFsiSION  KEY  (7) 

CC3790 

onv'3==l-/T2 

T 

003800 

C 

«IWFC,''  fq’  bo'JNOTNG 

003910 

TF(KrYc,)  .>!r,T) 

’0  Tr  10 

0C3320 

26 

TF(40!;  ("IPi/T)  .GT. 

F-hIf'- (Yrv(2)  ,0) ) GO  TO  f O 

Ct  3930 

^'•(sq  ; (P2-'\/0)  ,L  Y. 

^MTFT  (I.rv(7)  , ,1)  .4  .yy  YVY.NE.O)  jj  73  nj 

:C394C 

••8 

pyiVOyOHIFT  (S-(:fT  ( i?-'  '0,  -''EY  ( ’ 1 ) , <FY  ( ■•)  ) 

OOTASO 

’7?V0=FuTF7(rw:rT (C2  ',n,.^CY(i) 1 ,<fy(S)  ) 

r 0 9*60 

F"*UR*I 

: j 3J7C 

C’ 

Gl  or|M“'P  lY  »)''!< 

003990 

■»  1 

Y‘r=i 

OC  9 3-90 

:y(O207T.LY.3.) 

T'jr  I - 1 

C C TR'^O 

3’Y>)  J - r-  ( GMT  FT  fGxv^  T (P  ^cv,  -1  ))  t -1)  . <'Y  (')  ) 

; r » j < a 

’’■(SHIFT  (twift  ( P2GVG,  1 2)  , -12)  . t O.t)  )ALL  p )'  '3-  ■,(  ^T'TVO) 

«nn:o 

GY  TY  ’G 

:C.{99C 

n ) 

tf|R?y.,y,--.3.) 

P?GV''fSHIFT  ( YF  Y(  ’)  , i ) 

0 0 2 '-3 

--tjy-y-n.i'.' . ) 

F >Pv.'''r-3HTFT  ( •F''  ( ’)  . C) 

:c3'>Ga 

TF(<  rv(P)  )) 

fall  ovFPrLm'C'TiN-.L'iF'ir.?  ■ ^iyisionv 

0 3 9Y5C 

GO  ’3  29 

C-.'T'Y7C 

n 1 

PYTY3I3 

UCTIsC 

■■'ll  '.I'l”'  f i_  ' I » .r  - 1 • t ‘jr , I 'IT,  ■ < Y ’ / r S ’ y*-  ) 

Y 3 > )n 

-A  T n — . 

^ e t.  r\  • <8 

166 


I w * 

ENT 

pent  PUNCTION  ^invodl, !^2, ROUTINE, L^!CHT,IFNL,lfEY) 

nT‘«ENSI0N  KEY(7) 

mvosii/R? 

C ONECY  for  PO’INOING 

IFr<EY(4)  .N<^.0)  GO  YO  30 

2f>  TC(ATS('’invn)  .r.T.<:HIcT()'cY(2)  ,0))  GO  TO  100 

IP  (AOS(RlOvn) .LY .S“IFT (KEY (3) ,C) . A.RIOVO.NE. C)  30  TO  110 
28  RnVOsSNIPT  (SHIFT  (R10vn,-KEY  (7)  ) , KEY  (7)  ) 

IP(IFNL.FO.l)  >’lOV0  = SHIFT (SHIFT (RIO VO, -KEY (? ) ) , KEY (S)  ) 

RETURN 

C GO  ROUNO  IT  NOW 
30  T‘r«l 

IF(RlOVO.LT.0.)  IMC=-1 

"lOVOs SHIFT (SHIFT (SHIFT (PIOVO,- (KEY (7)-l)) TlNr,-!) ,<pY(7) > 
IF(SHIFT(SHYFT(S10VO,12) ,-12) .EQ. 0)  CALL  ROUNOER(RIOVO) 

GO  TO  25 

100  TF(R10VO.CT.O.)  siovr=SHIFT(KEY(2) ,3) 

IF(RinV0.LT.0.>  RIOVD^-SHIFT (KEY( 2) ,3) 

IFiKEY(5) .HE.Q)  CALL  PVFFFLH(POUTINE,LMrNT,8HOIVl3rON) 

GO  TO  29 
110  PlOVO=0 

TF(KEY (5)  .NE.O)  CALL  UNORFLH (ROUTINE, LKCNT , BHOIVISION) 

GO  TO  29 
FNY 

PEAL  FUNCTION  R2FXP(R1, T2, ROUTINE ,LNCNT, IFNL, KEY) 

OIhEksion  KFY(7) 
o?-XOsP.l*»I? 

IF(KFY (4) .MF. ))  GO  T"  30 
2fi  IF(AOS (R2EYP) ,GT. shift  (K£Y(2) , 0) ) GO  TO  IOC 

IF(  AOS  (RCEyd)  .lt.  shift  (i^CY(3j  ,0)  .A. P2EXP. NE.O)  GO  TO  110 

Z.t  O’Exot  SHI  ft  (SHIFT  (R?EXP,-KEY(  Y)  ) ,KFY  (7)  ) 

TF(IFnL.fo.I)  R2FXP=EHIFT(SHIFT(P.2EXP,-KFY(6).'  ,KEr(5)  ) 
PFTURN 

r GO  ROUNO  IT  NOW 

30  ING=1 

I"(R2fxp.L‘'.C.)  TNCf-1 

“?fxo=Sh:fT( shift (SHIFT (R2EXB, -(KEY ( 7 ) -1 ) ) ♦ INC , -1 ) , <f Y ( 7 ) ) 
YF(SHIFT(SwrFT(D2Fxo,ij)  ,-i?)  .£(>,0)  call  P.0l'N0Fj(^2rxP) 

GO  TO  29 

111  TF(o:,fxb  .GT  . c . ) BYc  yofSHIFT  (K  EY  (2  ) , C ) 

TcfoiFxo.LT  . p . ) airyai.^MIFT (KFy ( 2)  , 0) 

y'*(K.FY  (C)  .NY.O)  C;LL  0VFOFLW(=0'!T:*iF,LMrriT,  VHEx-'ONEnT) 

GO  TO  29 
111  PYFyosO 

TF(KEY  (5)  .NF.U  C:lL  UNORFLM(POUY:(IF,LNCMT,eHFXOONFNT) 

GO  TO  29  ' '• 

r.p 

F'fiCTION  rASG-J(  II , POUT  !•  f ,L fiS'lT  , IFNI.  , KEY  ) 

PI  'ENSION  kfy  {■') 

’AGGN*I1 

YF( I I .GT  .kfy ( 1)  1 YaSG*rYFY(l) 

"'I  I.LT  . -''Fv  (1  ) ) I.)«  G-  f-KEY  ( 1) 

Y'TUR'I 

* -*i  1 

a<<--  ('j,rc(,TT,  F,f.-f(r,IFNL,<FY) 

' 'I-,  Ff  » ( - ) 

'll  o»l»*»-rFY(Ftv(;),0) 

.1  »>  I -y  (■•),  1) 

• • * • > ' . > .i.  PR  1*1 

■ - - . • r I , »F»  I y ) I 


004010 
0C4020 
0C4030 
004040 
004050 
C04Q60 
004070 
304090 
OC4C90 
004103 
0C4110 
004120 
004130 
004140 
004150 
0C4150 
0C417O 
004160 
004190 
004200 
004210 
3C4220 
304230 
004240 
004250 
004260 
304270 
CC4290 
004290 
OG4300 
004310 
004320 
0C4330 
CC4340 
004350 
0C4360 
CC4370 
C0439C 
C 04390 
004400 
0 3'»410 
004420 
C 0 4-30 
C 04440 
J C44GC 
004460 
004470 
004480 

0C44<)3 
004=00 
004=10 
Cf)4  = 20 
.04530 
3:4041 
■:C4G50 
J04G40 
1 04=70 

.u.«gag 

0C4G90 
0:4410 
0 0-410 
y 04420 
V 430 
'<  .4*: 
. .44=: 


F -l.."*  t 


• . ^ • 


• I (.  AW  ■<  * \ 0 

{•*»  .MF.7)  r,0  Tf  3a 

T-{A1S  .'•.’.'^MTfT  (t'£Y(2)  ,0) ) '•■0  TO  130 

T^cansc^^.e^vp)  .L'^.'^htct  (k-cYc,)  ,J,  ,J,CPEXO.HE.O)  GOTO  IIC 

PPtXPaSMrrT (?wrPT (pcrxP.-<rY,7, j ,<£VC’)  ) 

TPCIFML.EO.l)  ?pexo  = S‘^IFT(SHIFT  (RPEXP,-i:eY(^)  ) ,XPY  (5)  ) 

9F’’IJRN 

TSOal 

tP(R9PX».LT.O.)  TNC=-l 

P9'•XPaS^rF’■^‘:wrPT(  SHIFT  (PPPXP.-fXEY  (7)-J  >)  ♦INH.-t ),  yeY(7)  ) 
TP(SHtPT(5MrFT(POFXP,l2) ,-12) .EO.O)  CALL  ROUNOE? ( 

GO  TO  25 

IF(n:iEXP.r,T.C  .)  b9EXP  = SHIFT  (KEY  (2)  ,1) 
xPt ROEXP.lt. a.)  oREXCr-SMIFT (YEY (2i  ,C  ) 

ro(KFY(5) .ME.O)  CALL  OVERPLH (ROUTINE, LNCNT , 8 HEXPONEUT ) 

GO  TO  25 
OP-XP=0 

TF(<EY (5)  .NE.O)  CALL  UNORFLW (ROUTINE, LNCNT , SHEXponERT ) 

GO  TO  25 
E'(0 

op)L  FUNCTION  91EXP(Ii,P2, ROUTINE, LNCMT,IFML, KEY) 

OT-ENSION  KFY(r) 

P1EXP=I1»*02 

TP(KEY  (V)  .NE.D  CO  TG  30 

IO(APS(RlEXO)  .'.T.Chtft  (i^PY(P)  , 0))  GO  TO  IOC 
IG(AeS(RlPXP)  .LT.EhIFT(i<^EY(?1  ,Q)  .a.PlEXP.NE.O)  GO  TO  llC 
Pl'’XPsSHIFT(<:HIFT(»iexp,-xfY(.')  ),KKY(7)  ) 

tp(IonL.‘^0 .1)  oiEXP=GHIFT(SHIoT(Pl£XP,-K£Y(6)) ,Kir(5) ) 

RETURN 

TNC«1 

TF'RlEXP.LT.C.)  ivc=-l 

PtTXP=G“TrT  ( FHTrr  ( t^HTCY  (Pipyo,  , (xey  (7)  -1  ) ) fine  , -1)  , XPY  ( 7)) 
t-^(Sh:pT  (CWTC- (PI -YP,  12)  , -12)  .EO.  3)  CALL  ^ OUM’'" > ( = 1 ^ X P) 

GO  TO  Ei 

-'■(RlFvo.r.r,,,  , F1lXG  = ghifT(y-Y(?>  ,3) 

T’^fPlE/P.I'.".)  sicy  ni.cuYCT(<rY(7)  ,3) 

TPC<ET(;)  .ME.:)  CALL  rv‘^RFL'.((KOUTINE,LNC‘(T,iH£xFJNTNr) 

GO  TO  25 
’ITyPaO 

T"(<EY  (3)  .ME.O)  call  UM0'’Fm(O0UTINE  , L*'CnT  , SHFXPONFnT  ) 

GO  TO  25 
G>|0 

<:'n90MTIMF  RO'IMOg?  (PMllN) 

'’ata  oMT/ar“;.c'::'’L'nT.OG?a3oo/,7E”o/ ’■’7-3  77  77'»77-F7-,--.-.7'r,/ 
'’A'A  P1777/17-Y  i--3^;;crr;,;;3Cr37 

OAT.a  RMnQ/e.0''J7’~77  77  7Y77-’777/P/ 

irfCHT cT(p Ml <-,■;). TO, <;hi'’''(?i''”,o))  pM"“T2':o:i.03cao''';‘30Cv''': 

rr  ( siiT  cr  ( P'"|“  . ) ) . T.)  .CUTTY  (Tf. ; 33  , ,•  ) ) c)M|'.'a-,rT7T77.77,'-'7?;-*.'' 

TCfCNTF’ (9M'|-’,  1)  .l,T  .“)  ^''''fP'iSO.AMA.T-RO 

’■T(GHrGT(R'.'JM,  1 ) .GT.  3)  7VJ".OR.OME 
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Appendix  D 


Proposed  In-line  Modification 

A proposed  change  to  the  existing  n-bit  simulation  tool 
would  replace  all  subroutine  calls  generated  by  the  present 
preprocessor  program  with  in-line  code  which  would  perform  the 
n-bit  wordlength  effects.  Although  this  modification  would 
result  in  large  memory  requirements  for  the  resulting  prepro- 
cessed program,  the  execution  time  should  be  reduced  by  one- 
third.  The  in-line  code  would  be  tailored  to  the  user  options 
in  effect  for  that  particular  program  module.  This  would  meam 
fewer  decisions  would  be  made  than  in  the  present  simulation 
tool,  during  execution  of  the  n-bit  simulated  version  of  the 
algorithm.  Also,  subroutine  linkage  time  would  be  save  since 
the  n-bit  wordlength  effects  were  no  longer  performed  through 
calls  to  subroutines. 

This  modification  would  require  the  detection  of  the 
CALL  SETXBIT  subroutine  for  each  program  unit.  The  NCN- 
EXPFESSICK-HANLLER  module  could  be  modified  to  detect  CALL 
SETNEIT,  then  call  a handling  routine  which  could  be  devel- 
oped to  analyze  the  user  options  expressed  as  arguments  in 
SETNBIT.  The  same  SETKEIT  subroutine  which  is  used  presently 
could  compute  the  key  values  and  pass  them  via  CCM:  ON  to  the 
subroutines  of  the  preprocessor  which  output  the  n-bit 
simulation  code  for  the  arithmetic  expressions. 
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The  in-line  code  which  could  perform  the  n-bit  simulation 
effects  for  the  various  arithmetic  operations  and  operand  data 
types  is  shown  in  Figure  D-1 . The  underlined  words  in  the 
figure  would  be  filled  in  by  the  handling  routine  with  each 
n-bit  simulation  situation. 


For  Real 
Results: 

If  Round 
Option 
is  on 


object  = varb1  operation  varb2 

object  = SHIFT(SHIFT(SHIFT(object  .-(key?  -1))+ 
(-1**(  SHIFT  (object  .A:vD.  ) , 1 ) ) , 

-D.key?  7"^ 

IF(SHIFT(SHIFT(object  .12). -12) .eq.O) 

'CAlI  I^CUND£K( object) 

IF(  AES(object) .GT.key2  .OF . (A£S(ob ject) .LT .key3 
. AND. AbS( object) .Mi .0) 

CALL  0VRFLWR(ob ject  . operation  .routine. 
line#  . key2  ) 

object  = SH1FT(SHIFT( object  .-key6/7) . key6/7) 


For  Integer 

Results:  object 


varbi  operation  varb2 


IF (IABS( object) .GT.keyl) 

CALL  CVPFLWiTobj'ect  .operation  .routine. 
line#.  keylT 


Figure  D-1  Fossible  In-line  Code  For 
Real  and  Integer  Arithmetic  Operation  Results 


To  output  the  code  shown  in  Figure  D-1,  the  GETIRTS  emd 
SNGLASC  modules  would  have  to  be  modified  slightly.  Speci- 
fically, the  merge  operations  of  each  module  would  be  replaced 
by  the  actual  arithmetic  operation  build  followed  by  a sub- 
routine call  to  a new  module.  This  new  module  would  have 
standard  outputs  (shown  by  the  non-underlined  portions  of  Fig- 
ure D-l).  Key  values  (computed  by  SETNBIT)  would  be  placed 
in  the  standard  positions.  The  calling  module  would  have  to 
pass  to  this  new  module  the  object  name,  object  data  type,  and 

In  addition,  the  new  routine  must 
70 
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be  passed  values  computed  by  the  SETNBIT  subroutine.  In  Figure 
D-1 , if  the  final  flag  was  turned  on,  key6  would  be  the  value 
positioned  in  the  last  statement.  If  this  was  not  the  final 
assignment  key?  would  be  placed  there,  so  a multi-precision 
mantissa  would  be  saved  when  the  in-line  code  is  executed. 
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